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FOREWORD 


The  collection  of  papers  in  this  publication  covers  the  presentations  made  at  the 
Third  Tropospheric  Refraction  Effects  Meeting,  sponsored  by  the  Electronic 
Systems  Air  Force  Systems  Command  of  the  Air  Force.  It  should  be  noted 
that  since  the  speakers  were  limited  in  time,  they  would  not  cover  all  of  the 
details  on  their  work  at  that  time.  The  papers  contained  in  this  document,  in 
some  cases,  extend  the  details  of  the  oral  presentations  in  that  they  much  more 
comprehensively  treat  the  topics.  Furthermore,  the  results  of  the  most  recent 
measui'ements  data  have  been  included  where  feasible. 

The  publication  of  these  proceedings  does  not  constitute  official  Air  Force  or 
MITRE  endorsement,  acceptance,  or  approval  of  the  techniques  proposed  and 
used  in  the  work,  or  of  the  results  reported.  Rather,  it  is  a  record  of  what 
work  was  covered,  and  provides  a  useful  collection  of  basic  data  and  measure¬ 
ment  results  to  stimulate  an  exchange  of  ideas  which  will  advance  the  state-of- 
the  art.  In  addition,  it  serves  as  a  ready  reference  of  the  most  up-to-date 
thinking. 

I  wish  to  express  my  appreciation  to  the  representatives  of  the  various  agencies 
for  their  excellent  cooperation,  to  all  the  participants  of  the  meeting,  and  to 
the  Range  Measurements  Division  of  the  Directorate  of  Aerospace  Instrumen¬ 
tation  for  suggestions  and  assistance.  I  thank  the  members  of  the  Technical 
Services  and  the  Radar  Systems  and  Techniques  Departments  of  The  MITRE 
Corporation  for  the  excellent  support  they  have  given  me  during  the  meeting 
and  in  the  preparation  of  this  publication.  The  success  of  the  meeting  was 
greatly  enhanced  by  their  assistance.  No  acknowledgement  would  be  complete 
without  the  mention  of  Mrs.  Dorothy  Pierce,  who  typed  this  document,  of 
Miss  Margy  Lawlor,  who  spent  many  hours  in  the  preparation  of  the  artwork, 
and  to  Mrs.  Patricia  Chatta,  who  coordinated,  planned,  edited  and  guided  this 
document  through  all  stages  of  its  preparation. 


L.  J.  Galbiati, 
Chairman 
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(in  Two  Volumes) 

VOLUME  I 

Results  and  Status  of  Tropospheric  Effects  Measurement  Fh'ograms 

ABSTRACT 


This  publication  covers  the  presentations  made  at  the  Third  Tropospheric 
Refraction  Effects  Meeting,  sponsored  by  the  Electronic  Systems  Division  of  The 
Air  Force  Systems  Command.  The  Meeting  was  held  at  The  MITRE  Corporation 
on  July  28,  29  and  30,  1964.  The  proceedings  are  being  published  in  two  volumes 
for  ease  in  handling;  the  papers  covering  the  work  on  tropospheric  refraction 
sponsored  by  the  Electronic  ^sterns  Division  in  support  of  metric  instrumentation 
on  the  national  missile  ranges,  are  contained  in  Volume  1.  While  these  papers  are 
complete  in  themselves ,  many  of  the  specific  details  on  the  work  covered  in  the 
proceeding  of  the  second  Refraction  meeting  (ESD-TDR-64-103)  have  not  been 
repeated.  The  first  volume  also  contains  papers  on  the  Findings  and  Recommen¬ 
dations  of  the  National  Academy  of  Science  "Panel  on  Tracking  Data  Analysis"  in 
the  Area  of  Tropospheric  Fh'opagation  and  Refraction,  Refracto meter  and  Humldi- 
ometer  Development  of  the  Naval  Research  Laboratory,  Spaced  Refractometer 
Program  of  Eastern  Test  Range,  and  Canadian  Tropospheric  Measurement 
Activities. 

Presentations  on  certain  refraction  activities  at  the  Pacific  Missile  Range, 
the  National  Bureau  of  Standards,  Ft.  Huachuca,  Ballistic  Research  Laboratory, 

M.  I.  T.  Roimd  Hill,  Cold  Lake  Canada,  on  new  measurement  techniques,  and  on 
Tropospheric  and  Ionospheric  Effects  on  Interferometer  Systems  have  been 
included  in  Volume  H. 
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REVIEW  AND  APPROVAL 


F>ublication  of  this  technical  documentary  report  does  not  constitute  Air  Force 
approval  of  the  reports  findings  or  conclusions.  It  is  published  only  for  the 
exchange  and  stimulation  of  ideas. 

ROY  D.  RAGSDALE 
Col,  USAF 

Director,  Aerospace  Instrumentation 
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OPENING  REMARKS;  TROPOSPHERIC  REFRACTION  MEETING 


Col.  R.  D.  Ragsdale* 

I'm  Colonel  Ragsdale,  Director  of  Aerosp>ace  Instrumentation.  I  would  like 
to  welcome  you  men  this  morning  to  the  Third  Tropospheric  Refraction  Effects 
Meeting.  One  thing  that  has  become  very  obvious  during  my  briefings  with 
General  Davis  of  the  National  Range  Division  (NRD),  and  elsewhere  around  the 
country,  is  that  everyone  is  very  enthusiastic  about  the  active  refraction  effects 
study  programs  we  have  at  the  Electronic  Systems  Division.  Everytime  we  have 
briefed  General  Davis  he  has  commented  that  this  is  one  of  the  most  impKjrtant 
things  we  can  do  for  him.  The  NRD  staff  is  very  interested  in  making  sure  that 
all  the  numbers  derived  by  you  folks  are  fed  into  their  error  models  as  soon  as 
possible.  They  want  to  make  sure  that  the  refraction  programs  we  undertake 
are  related  back  to  measurement  accuracy. 

I  would  like  to  say  you  are  suppKJrted  very  well  and  at  very  high  levels  in 
your  efforts  in  this  refraction  business;  great  things  are  expected  from  the  effort 
in  this  area.  I  will  now  turn  you  over  to  Captain  Eames,  who  is  the  Air  Force 
head  of  the  project,  and  let  you  get  on  with  your  business. 


*Director,  Aerospace  Instrumentation,  USAFSC,  ESD,  L.  G.  Hanscom  Field, 
Bedford,  Massachusetts 
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INTRODUCTION 


Capt.  E.  G.  Eames* 

This  three-day  meeting  is  unclassified,  and  has  been  organized  for  the 
primary  purpose  of  presenting  you  with  information  on  the  many  significant 
activities  in  the  area  of  tropospheric  refraction  effects.  A  second  purpose  is 
to  provide  you  with  the  opportunity  to  meet  each  other  in  order  to  discuss  the 
various  aspects  of  the  problem.  We  have  limited  the  number  of  presentations 
at  each  session  in  order  to  provide  discussion  periods.  The  Proceedings  of  the 
second  refraction  meeting,  held  last  November,  have  been  distributed  to  all  the 
attendees  of  that  meeting.  Extra  copies  will  be  made  available  in  the  lobby. 
Proceedings  of  this  meeting  will  be  distributed  to  all  the  attendees  at  these 
sessions. 

The  Directorate  of  Aerospace  Instrumentation's  interest  in  refraction  is 
in  line  with  its  mission  to  support  the  National  Ranges.  We  are  presently  funding 
a  continuing  comprehensive  integrated  program  with  the  National  Bureau  of 
Standards  (NBS),  MITRE,  and  the  Air  Force  Cambridge  Research  Laboratory 
(AFCRL). 

The  objectives  are  to  develop  techniques  to  minimize  the  refractive  effects, 
determine  the  extent  of  the  effects  on  optical  systems,  obtain  real-time  correc¬ 
tion  techniques,  and  develop  noise  predictors  for  tracking  systems.  Again,  the 
emphasis  is  on  feeding  more  of  the  data  obtained  as  direct  inputs  into  the  error 
models  in  order  to  derive  better  corrections  to  the  range  systems. 


*  HQ  ESD  ESRIM,  L.  G.  Hanscom  Field,  Bedford,  Massachusetts 
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The  program  combines  three  separate  approaches;  namely,  fixed-end  point 
measurements,  dynamic  measurements  at  the  Eastern  Test  Range  (ETR),  and  a 
real-time  correction  device.  I  will  now  turn  the  meeting  over  to  Dr.  L.  J.  Galbiati, 
who  will  give  you  a  more  comprehensive  description  of  each  of  these  approaches. 
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REVIEW  OF  ESD  REFRACTION  TASK 


L.  J.  Galbiati* 


INTRODUCTION 


The  total  effort  under  the  refraction  task  at  the  Electronic  Systems 
Division  (ESD)  has  been  divided  into  three  parts;  as  shown  in  Fig.  1.  While 
the  individual  work  agents  have  specific,  limited  objectives,  an  effort  is  made 
to  coordinate  and  direct  tham  so  as  to  achieve  composite  results  applicable  to 
the  operational  problem  of  the  National  Range  Division  (NRD). 


Fig.  1  .  Overall  ESD  Task  Approach 
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MEASUREMENTS  AT  THE  EASTERN  TEST  RANGE 


Dr.  R.  M.  Cunningham  is  directing  the  work  being  performed  for  ESD  by 
the  Cloud  Physics  Branch  of  the  Cambridge  Research  Laboratories  (CRL). 

This  effort  involves  obtaining  dynamic  meteorological  measurements  at  the 
MISTRAM  site  of  the  Eastern  Test  Range  (ETR),  as  well  as  analysis  of  the  data. 
The  data  is  obtained  by  use  of  special  instrumented  aircraft,  radiosonde,  wire- 
sonde,  and  ground-based  meteorological  instrumentation.  (Fig.  2). 


At  the  previous  meeting  (November  1963),  I  indicated  that  we  were 
evaluating  methods  of  presenting  the  data  in  a  useful,  meaningful  fashion. 
During  the  past  six  months,  we  have  concentrated  on  one  unique  method;  the 
digital  representation  of  the  two-dimensional  variations  -  in-refractivity 
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profile.  The  profiles  are  very  difficult  to  generate  in  that  they  require  exten¬ 
sive  meteorological  measurements  over  the  specific  terrain  of  interest,  and 
interpretation  of  the  data  obtained  by  an  experienced  meteorologist.  However 
profiles  for  three  different  days  have  been  generated  to  date,  and  the  way  in 
which  they  are  being  utilized  will  be  covered  by  R.  K.  Crane.  With  these 
representations,  it  is  possible  to  use  a  computer  to  investigate  such  effects  as 
the  beam-sweeping  movement  resulting  from  tracking  a  moving  target,  different 
wind  velocities,  special  meteorological  conditions,  and  different  environmental 
conditions  for  the  individual  rays  of  an  interferometer  system. 

Dr.  Dryden,  in  one  of  his  studies,  indicated  the  effect  of  a  refractivity 
gradient  along  the  legs  of  the  MISTRAM  system,  and,  from  time  to  time,  the 
validity  of  using  one  profile  for  correcting  the  entire  interferometer  system 
has  been  questioned.  Little  progress  has  been  made  resolving  questions  in  this 
area  to  date;  we  have,  therefore,  endeavored  to  obtain  data  necessary  in 
arriving  at  a  solution  by  obtaining  simultaneous  vertical  refractivity  profile 
measurements  at  MISTRAM  Central  and  MISTRAM  end-points,  some  twenty 
miles  distant,  by  using  three  instrumented  aircraft.  The  details  on  this  phase 
of  the  program  will  be  covered  by  D.  L.  Ringwalt.  Later,  Dr.  Cunningham 
will  discuss  the  scales  and  types  of  atmospheric  refractive  index  anomalies. 

J.  H.  Meyer  will  explain  the  mechanism  of  generating  the  profile,  and  R.  K 
Crane  will  show  results  of  experiments  utilizing  the  digital  atmospheric  profile. 
These  papers  summarize  the  ETR  phase  of  our  program. 

FEASIBILITY  TEST  OF  LINE  INTEGRAL  REFRACTOMETER 

The  feasibility  field  tests  of  the  MITRE  Line  Integral  Refractometer  (LIR) 
have  included  measurements  on  a  number  of  paths  at  different  elevation  angles 
between  0.  2  and  4.  9  degrees,  as  indicated  by  the  top  diagram  in  Fig.  3.  In 
November  1963,  preliminary  results  of  the  measurements  on  a  level  path,  as 
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indicated  by  the  bottom  sketch  in  Fig.  3,  were  presented.  This  technique,  as 
you  may  recall,  involves  the  determination  of  the  differential  phase  between 
two  coherent  signals  (15.  6  and  31.  2  gc)  transmitted  over  the  actual  path  of 
interest.  The  differential  phase  shift  is  used  to  infer  the  differential  transit 
time  of  the  two  signals  and,  therefore,  provides  desired  information  on  the 
line  integral  of  the  refractivity;  the  line  integral  of  refractivity  is  the  correc¬ 
tion  one  has  to  make  to  metric  tracking  system  measurements. 


The  sensitivity  of  the  measurement  instrumentation  developed  in  this 
task  has  been  very  satisfactory,  and  the  field  measurements  are  the  only  ones 
of  this  type  to  my  knowledge.  H.  M.  Richardson  will  describe  details  of  field 
tests.  He  will  interpret  the  results  in  terms  of  tracking  requirements  and 
applioability  to  the  general  problem  of  the  line-of-sight  path  corrections. 


4 


Future  (early  1965)  plans  for  this  effort  include  measurements  to  determine 
the  effect  of  the  oxygen  constituent  in  the  transmission  path.  These  will  be 
obtained  by  measuring  the  differential  phase  of  another  pair  (45  and  90  gc)  of 
coherent  signal  frequencies  bracketing  the  oxygen  absorbtion  region.  The 
complete  water  vapor-oxygen  refractometer  should  be  tested  sometime  during 
the  first-half  of  1965. 

One  of  the  basic  differences  between  the  LIR  measurement  and  the 
National  Bureau  of  Standards  (NBS)  radio  frequency  phase  measurements 
should  be  pointed  out.  The  LIR  measures  the  differential  phase  shift  between 
two  coherent  signals  at  different  frequencies  (15.  6  and  31.  2  gc),  whereas  the 
NBS  system  measures  the  phase  shift  of  a  signal  at  one  frequency  (9.  4  gc). 
This  difference  in  technique  results  in  system  where  the  end-points  of  the  LIR 
do  not  have  to  be  fixed  and,  therefore,  can  be  mounted  on  an  aircraft 
(see  Fig.  4).  No  arrangements  have  been  made  to  do  so  under  the  current 


Fig.  4.  MITRE  —  Line  Integral 

Refractometer  Technique 
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program.  The  NBS  system,  on  the  other  hand  requires  fixed  end-points  to 
separate  out  a  very  small  portion  of  the  total  effect.  The  LIR  is  the  only 
potentially  practical  technique  available  for  obtaining  meaningful  data  on  the 
actual  microwave  refractivity  characteristics  of  the  entire  atmospheric  path 
for  a  wide  range  of  elevation  angles. 

OPTICAL  AND  RADAR  FREQUENCY  ACTIVITY 

The  ESD  refraction  activities  at  NBS  includes  work  at  both  optical  and 
radio  frequencies  (Fig.  5).  We  have  consolidated  the  data  requirements  for 
different  items  of  range  instrumentation  in  which  ESD  is  involved,  and  have 
outlined  work  statements  which  provide  all  the  necessary  information,  but  keep 
the  number  of  separate  measurement  efforts  at  a  minimum.  That  is,  each 
phase  of  the  NBS  work  does  not  necessarily  reflect  the  problem  associated  with 
only  one  specific  metric  system. 


Fig.  5.  National  Bureau  of  Standards  —  Fixed  End-Point  Technique 
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The  fixed  end-point  radio  frequency  phase  measurement  program  field 
work  at  Maui,  Hawaii  has  been  completed  and  will  be  covered  by  Dr.  M.  C. 
Thompson.  Optical  depolarization  measurements  have  been  performed  and 
will  be  reported  on  by  Dr.  J.  R.  Meyer-Arendt;  this  work  was  initiated  about 
a  year  ago,  and  includes  the  literature  survey  report  which  is  being  published 
by  NBS.  The  conclusions  of  this  survey  were  published  in  the  Proceedings  of 
the  second  meeting.  *  Early  this  year,  NBS  was  requested  to  study  the  problem 
of  fading  on  the  MISTRAM  data  links.  They  have  been  studying  the  MISTRAM 
records,  and  H.  T.  Dougherty  will  report  their  findings  and  recommendations. 

SUMMARY 

The  details  covering  the  history,  objectives,  and  organization  of  the  over¬ 
all  refraction  task  are  covered  in  the  Proceedings**  of  the  second  meeting,  and 
will  not  be  repeated  in  this  publication.  Papers  on  the  results  and  progress 
achieved  during  the  past  six  months  on  each  phase  of  the  ESD  program  will  be 
presented  by  the  individuals  responsible  for  the  work. 

We  have,  in  this  short  period  of  time,  been  able  to  instrument  aircraft 
and  collect  valuable  data  at  the  Eastern  Test  Range.  We  have  generated  the 
first  two-dimensional  digital  atmospheric  profile,  applied  it  to  an  operational 
system,  and  estimated  atmospheric  errors.  The  instrumentation  for  deter¬ 
mining  differential  phase  data  between  15.  6  and  31.  2  gc  signals  has  been 
developed,  and  the  first  long  path  field  measurements  have  been  completed. 

The  feasibility  testing  of  the  line  integral  refractometer  to  date  strongly  indi¬ 
cates  its  high  value  as  a  calibration  tool. 


*See  ESD-TDR-63-103,  Vol.  I,  pages  1-131  to  1-142. 

**ESD-TDR-64-103,  Vols.  I,  H  and  HI. 
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The  comprehensive  series  of  phase  measurements  at  Maui,  Hawaii,  have 
been  completed.  These  include  the  extensive  set  of  measurements  with  refrac- 
tometers  spaced  at  considerable  distances.  Publications  related  to  the  work 
conducted  under  this  program  include: 

Effects  of  Clouds  and  Moist  Layer  Roughness  on  Range  Rate 
Errors  of  Radio  Interferometer  Tracking  and  Guidance  Systems, 
AFCRL-63-706; 

Instrumentation  to  Measure  the  Effect  of  Meteorological 
Parameters  on  the  MISTRAM  System  and  Preliminary 
Data  Reduction,  CRL  Memo  11/3; 

Line  Integral  Ref  rac  to  meter,  MITRE  Report  SR- 112; 

Optical  Scintillation,  NBS  Tech.  Note  No.  225;  and 

Proceedings  of  the  Second  Tropospheric  Refraction 
Effects  Meeting,  ESD-TDR-64-103,  Vols.  I,  II  and  III. 

The  investigators  should  be  commended  for  the  outstanding  cooperation 
they  have  given  me  toward  achievement  of  the  overall  technical  objectives  of 
this  ESD  task.  We  initiated  a  rather  extensive  measurement  program,  and 
have  accumulated  much  experimental  data  directly  related  to  the  refraction 
problem.  Hence,  a  continuous  flow  of  results  can  be  anticipated  over  the  next 
year  or  so  from  the  analysis  of  the  data  we  have  collected  under  this  program. 
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SIMULTANEOUS  REFRACTIVE  INDEX  MEASUREMENTS 


BY  THREE  AIRCRAFT 

T.  S.  Weaver  and  D.  L.  Ringwalt* 


INTRODUCTION 

On  16  and  17  June  1964,  three  aircraft,  equipped  with  refractometers , 
obtained  simultaneous  refractive  index  measurements  in  the  MISTRAM  area 
while  flying  a  tight  V-formation,  and  also  while  separated  by  about  20  miles. 
These  measurements  were  obtained  in  preparation  for  and  in  order  to  determine 
the  problems  involved  in  more  extensive  tests  planned  for  early  August.  This 
work,  carried  out  at  Patrick  Air  Force  Base  (PAFB),  is  part  of  an  atmospheric 
research  program  aimed  specifically  toward  improving  the  tracking  accuracy 
of  MISTRAM.  The  data  to  be  presented  has  not  been  completely  checked  and 
consequently,  the  results  should  be  considered  as  preliminary. 

PURPOSE 

The  primary  purpose  of  these  tests  was  to  determine  the  feasibility  of 
airborne  intercomparison  and  calibration  of  the  refractometers  aboard  three 
separate  aircraft.  This  intercomparison  and  calibration  is  necessary  if 
accurate  simultaneous  refractive  index  measurements  are  to  be  obtained  in 
the  MISTRAM  area  by  the  three  aircraft  during  the  planned  intensive  measure¬ 
ments  in  August.  The  August  measurements  will  be  used  as  inputs  to  MISTRAM 
error  model  computations. 


^Electromagnetic  Research  Corporation,  College  Park,  Maryland 
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A  secondary  purpose  was  to  obtain  soundings,  at  MISTRAM  Central  and 
at  the  outlying  100  K.  West  and  100  K.  South  sites,  of  the  type  that  will  be 
obtained  in  August,  to  learn  of  possible  operational  or  technical  difficulties. 

EQUIPMENT 

Three  aircraft  were  used  for  the  tests: 

(a)  Aircraft  No.  30133,  a  C-130,  equipped  with  a  Texas-type 
refractometer  and  extensive  meteorological  equipment; 

(b)  Aircraft  No.  37802,  a  C-131,  equipped  with  an  ASH- 14 
refractometer;  and 

(c)  Aircraft  No.  37794,  a  C-131,  equipped  with  a  Texas -type 
refractometer  and  other  meterological  equipment. 

In  addition  to  the  three  aircraft,  the  facilities  of  the  Pan  American  Airways 
(PAA)  weather  station  at  MISTRAM  Central  were  used  for  radiosonde  and  wire- 
sonde  measurements. 

EXPERIMENTS 

A  number  of  experiments  or  calibrations  were  conducted.  The  aircraft 
maneuvers  are  described  in  the  next  section.  Each  of  the  experiments  listed 
on  the  following  page,  along  with  the  results,  is  discussed  later  in  the  report. 

The  experiments  consisted  of  refractive  index  measurements  in  the  three 
aircraft  as  they 

(a)  made  simultaneous  ascents  and  descents  while  flying  in  a 
V-formation, 

(b)  made  simultaneous  ascents  over  MISTRAM  Central,  100  K. 

West  and  100  K.  South,  and 

(c)  flew  in  a  level  flight  in  V-formation. 
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FLIGHT  PLAN 


The  same  flight  plan  was  used  on  both  days.  All  three  aircraft  departed 
form  PAFB.  The  C-130  aircraft  (133)  climbed  to  25,000  feet,  while  both  C-131 
aircraft  (794  and  802)  climbed  to  15,  000  feet.  About  16  miles  ENE  of  MISTRAM 
Central,  133  descended  to  15, 000  feet  to  join  the  other  two  aircraft. 

After  approximately  5  minutes  of  level  flight  at  15, 000  feet,  802  took  the 
lead  of  a  V-formation,  with  794  to  port  and  133  to  starboard.  The  aircraft 
commenced  a  4-minute  holding  pattern,  with  standard  rate  turns,  on  a  heading 
of  340/160  degrees  and  at  a  position  about  16  miles  ENE  of  MISTRAM  Central. 
While  maintaining  this  pattern,  all  aircraft  descended  in  formation  to  500  feet, 
at  approximately  1000  feet  per  minute,  with  802  calling  out  altitude  on  the  VHF 
radio  at  1000-foot  intervals. 

All  aircraft  maintained  the  500-foot  altitude  for  5  minutes  and  then  started 
to  ascend  at  500  feet  per  minute,  maintaining  the  4-minute  holding  pattern  during 
the  ascent  until  an  altitude  of  15,  000  feet  was  reached.  The  802  continued  to 
call  out  altitude  at  1000-foot  intervals.  Maintaining  15,000  feet,  133  proceeded 
to  100  K.  West,  794  to  MISTRAM  Central,  and  802  to  100  K.  South. 

At  a  time  indicated  by  133,  each  aircraft  began  to  descend  at  1000  feet 
per  minute,  using  a  figure-eight  pattern  centered  on  the  respective  MISTRAM 
stations.  The  figure-eight  track  was  centered  on  a  340/160-degree  line,  with 
turns  beginning  2  minutes  after  the  aircraft  passed  over  the  MISTRAM  stations. 
After  reaching  an  altitude  of  500  feet  over  the  MISTRAM  sites,  all  aircraft 
proceeded  to  a  point  approximately  10  miles  downwind  from  MISTRAM  Central 
while  climbing  to  2000  feet. 

With  802  as  the  lead  aircraft  in  a  V-formation,  and  with  794  to  port  and 
133  to  starboard,  the  aircraft  flew  toward  the  wiresonde  balloon  (at  MISTRAM 
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Central)  which  was  flying  at  about  2000  feet.  The  794  and  802  turned  left  approxi¬ 
mately  one-half  mile  from  the  balloon,  while  133  passed  as  close  as  practical  to 
the  balloon.  All  aircraft  rendezvoused  at  a  point  10  miles  downwind  from 
MISTRAM  Central  to  repeat  the  run. 

MISTRAM  Central  reported  winds  at  2000  feet  as  light  and  variable  on 
both  days.  As  a  result,  the  following  tracks  were  flown: 

June  16 


Track  No. 
Track  No. 

Track  No. 
Track  No. 
June  17 

Track  No. 
Track  No. 
Track  No. 


1  —  240  degrees  to  MISTRAM  Central  at  1900  feet; 

2  —  252  to  240  degrees  to  MISTRAM  Central  at 

1900  feet; 

3  —  340  degrees  to  MISTRAM  Central  at  1900  feet; 

4  —  45  degrees  to  MISTRAM  Central  at  1900  feet. 

1  —  30  degrees  to  MISTRAM  Central  at  2000  feet; 

2  —  270  degrees  to  MISTRAM  Central  at  2000  feet; 

3  —  40  degrees  to  MISTRAM  Central  at  2000  feet. 


OPERATIONS 


No  serious  operational  problems  were  encountered  in  flying  the  three 
aircraft  in  the  tight  V-formation,  or  in  making  simultaneous  descents  with  the 
aircraft  about  20  miles  apart.  The  greatest  difficulty  was  in  the  post-flight 
correlation  of  the  data.  A  better  method,  either  manual  or  automatic,  is 
needed  to  mark  all  recordings  simultaneously  with,  for  instance,  time,  for 
future  correlation  of  events. 
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FORMATION  ASCENTS  AND  DESCENTS 


The  formation  ascents  and  descents  of  the  three  aircraft  between  15,000 
and  500  feet  were  made  at  sea  about  16  miles  ENE  of  MISTRAM  Central,  where 
atmospheric  conditions  were  relatively  stable  and  where  good  correlation  of  the 
refractive  index  data  from  the  three  aircraft  could  be  expected.  If  the  atmos¬ 
phere  is  assumed  homogeneous  over  at  least  the  volume  in  which  the  three 
aircraft  are  flying,  then  any  lack  of  correlation  between  the  values  of  refractive 
index  recorded  in  the  three  aircraft  may  be  ascribed  to  instrumentation 
differences. 

The  air  temperature  difference  between  15, 000  and  500  feet  was  about  27 
degrees  C.  By  first  descending  to  warmer  air  and  then  ascending  to  the  colder 
air,  temperature  effects  upon  the  refractometer  measuring  cavity  are  cycled. 
Assuming  the  15, 000-foot  value  of  refractive  index  does  not  vary  in  the  time 
between  the  descent  from  and  the  return  to  15,  000  feet,  then  any  difference  in 
the  indicated  value  of  refractive  index  before  the  descent  and  after  the  ascent  is 
a  measure  of  the  hysteresis  of  the  measuring  cavity. 

Aircraft  barometric  type  altimeters  "lag"  the  actual  altitude  during  ascents 
and  descents,  depending  upon  the  rate  of  ascent  or  descent.  During  an  ascent, 
the  altimeter  will  indicate  an  altitude  lower  than  actual  pressure  altitude,  with 
greater  differences  being  indicated  for  rapid  ascents.  Similarly,  the  altimeter 
will  indicate  a  higher  altitude  than  actual  pressure  altitude  during  descent,  with 
greater  differences  being  indicated  for  more  rapid  descents. 

Corrections  for  altimeter  lag  have  not  been  made  in  the  data  to  be  shown 
in  the  figures.  Since  802  announced  the  altitude  at  1000-foot  intervals,  the  data 
for  each  aircraft  has  been  plotted  to  the  same  scale.  However,  it  should  be 
noted  that,  because  of  the  altimeter  lag  problem,  any  particular  altitude  as 
plotted  will  not  only  differ  from  the  true  pressure  altitude,  but  will  also  differ 
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for  the  ascent  and  descent.  The  overall  effect  of  this  is  to  shift,  in  indicated 
altitude,  a  particular  refractive  index  perturbation  between  ascent  and  descent. 

Although  corrections  for  altimeter  lag  have  not  been  made  in  the  data 
plots,  such  corrections  are  possible  and  will  be  made  at  a  later  date  in  the 
following  manner:  Both  133  and  794  made  comparisons  of  the  altitude,  as  indi¬ 
cated  on  both  the  barometric  altimeter  and  the  radar  altimeter,  during  the 
ascents  and  descents.  Since  the  radar  altimeter  has  no  lag  characteristics, 
it  is  thus  possible  to  correct  all  altitudes  to  a  common  and  true  altitude  scale. 

The  133  and  794  have  refractometer  systems  modified  in  a  manner 
described  by  Shaw.  This  modification  removes,  almost  entirely,  drifts  in  the 
system  output,  except  those  due  to  cavity  characteristics.  The  modification 
also  includes  provisions  for  precision  calibration  in  25  N-unit  steps.  The 
calibration  of  the  802  refractometer  is  accomplished  by  noting  the  difference 
in  the  indicated  refractive  index  while  adjustments  are  made  to  a  calibrated 
micrometer  tuning  screw  on  the  reference  cavity. 

Inspection  of  the  802  records  for  17  June  immediately  after  the  flight 
indicated  that  an  obvious  and  large  change  in  calibration  had  occurred.  A  new 
scale  (N -units/unit  recorder  deflection)  was  constructed  by  comparison  with 
the  133  data.  No  claim  is  made  that  this  correction  is  proper,  only  that  the 
corrected  scale  is  "reasonable.  "  This  fact  must  be  kept  in  mind  during 
inspection  of  the  data. 

Figures  1  and  2  show  the  over-water  descents  and  ascents  of  each  air¬ 
craft.  The  following  procedure  was  used  in  plotting  these  descents  and  ascents: 
In  all  plots,  the  absolute  value  of  refractive  index  at  15, 000  feet  was  normalized 
to  the  radiosonde  value  obtained  at  MISTRAM  Central  nearest  the  time  of  the 
flight.  The  radiosonde  data  is  plotted  at  the  right  in  Figs.  1  and  2. 
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Ascents  VS.  Descents,  16  June  1964 
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Ascents  vs.  Descents,  17  June  1964 


Because  of  separation  between  MISTRAM  Central  and  the  operating  area,  no 
particular  detailed  correlation  with  the  aircraft  data  is  expected.  The  aircraft 
descents  were  plotted  and  then,  without  changing  the  scale,  the  ascents  were 
plotted.  If  it  is  assumed  that  the  atmosphere  did  not  change,  the  refractive 
index  indication  should  return  to  the  initial  value  at  15, 000  feet.  It  is  seen 
that  the  refractive  index  value  at  15,  000  feet  does  return  to  the  initial  value 
within  3  to  4  N-units  for  all  aircraft  on  both  days  of  operation. 

Inspection  of  the  profiles  for  16  June  shows  first  that  this  was  not  a 
particularly  good  day  for  the  measurements.  This  kind  of  test  requires  a 
"smooth"  atmosphere,  and  the  profiles  show  that  this  was  not  the  situation, 
particularly  below  about  7000  feet.  Considering  that  no  corrections  have  been 
made  to  the  data  for  either  cavity  or  altimeter  characteristics  and  that  the 
atmosphere  was  not  "smooth, "  the  ascents  and  descents  for  the  three  aircraft 
compare  quite  well.  Note  particularly  the  region  between  about  6000  feet  and 
12,  000  feet.  In  this  region,  all  three  aircraft  retraced  their  profiles  within 
3  to  4  N-units. 

Deviations  between  ascents  and  descents  below  7000  feet  are  often  corre¬ 
lated  among  the  three  aircraft.  The  indicated  roughness  of  the  atmosphere 
probably  is  responsible  for  some  of  the  differences  between  aircraft,  even 
though  they  were  flying  in  a  tight  V-formation.  The  profiles  obtained  on  17 
June  show  somewhat  better  agreement  between  ascents  and  descents,  since 
the  atmosphere  was  appreciably  smoother  on  that  day  than  on  16  June.  When 
altimeter  lag  corrections  are  made  to  the  data,  their  effect  will  be  to  "lift"  the 
whole  trace  for  the  ascent  and  to  "lower"  the  whole  trace  for  the  descent. 

A  better  comparison  of  the  data  from  the  three  aircraft,  is  shown  in 
Fig.  3.  The  plots  of  the  data  obtained  while  the  three  aircraft  made  ascents  and 
descents  in  formation  are  superimposed.  Figures  1  and  2  mi^t  be  considered 


17 


4 


CO 

D) 


18 


Formation  Ascents  and  Descents,  133  vs.  794,  vs.  802 


as  internal  comparisons,  while  Fig.  3  is  a  true  comparison  of  the  three  air¬ 
craft  installations.  Altimeter  lag,  due  to  the  method  of  logging  altitude  in  each 
aircraft,  will  not  affect  the  comparisons  of  Fig.  3. 

As  explained  earlier,  the  802  data  used  the  133  data  for  calibration,  and, 
for  this  reason,  better  correlation  should  be  expected  between  133  and  802  than 
between  133  and  794.  This,  however,  is  not  obvious  in  Fig.  3.  The  tendency  for 
794  to  indicate  lesser  values  of  refractivity  than  133  might  be  due  to  the  method 
of  normalizing  the  data  to  the  radiosonde  value  at  15,  000  feet  on  descents. 

The  differences  between  the  various  profiles  are  displayed  in  Fig.  4; 
some  features  alreatfy  noted  are  shown  clearly.  The  differences  between  the 
133  data  and  that  for  the  794  and  802  have  been  plotted  against  altitude  for  the 
descents  on  both  days.  The  differences  are  less  at  the  higher  altitudes  where 
the  atmosphere  is  smoother.  The  plots  suggest  that  the  794  and  802  data  are 
better  correlated  than  the  133  and  794  or  the  133  and  802  data.  This  observa¬ 
tion  is  made  by  noting  the  similarity  in  the  shapes  of  the  two  plots  for  16  June 
and  the  two  plots  for  17  June. 

It  must  be  kept  in  mind  that  this  experiment  and  these  data  reductions 
have  been  performed  to  determine  the  problems  involved  in  obtaining  several 
simultaneous  profile  measurements,  using  several  aircraft.  These  profiles 
will  then  be  used  in  MISTRAM  error  model  calculations .  The  profile  data 
shown  in  Figs.  1  to  3  should  be  compared  in  a  manner  directly  applicable  to 
this  problem.  The  error  model  will  make  use  of  the  propagation  delay  of  a 
radio  signal,  relative  to  propagation  in  a  vacuum,  through  the  atmosphere,  as 
described  by  refractivity  profiles.  Table  1  lists  the  delays  for  vertical 
propagation  to  an  altitude  of  15,  000  feet.  Delays  for  propagation  to  the  same 
altitude  at  other  elevation  angles  may  be  determined  by  multiplying  the 
indicated  values  by  the  cosecant  of  the  elevation  angle. 
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Scxjnding  Differences  133  vs.  794  and  802 


Vertical  Propagation  Delays  to  15, 000-Foot  Altitude 


Difference  (ft. ) 

794  vs. 

802 
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0.  03 

-0.  02 

0.  03 
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794 
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0.03 
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Aircraft  No. 
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0.03 
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794 
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Descent  (D^^)  and 
Ascent  (D^)  Delays 
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Assuming  homogeneity  of  the  atmosphere  in  at  least  the  volume  in  which 
the  three  aircraft  are  flying,  the  differences  in  the  three  right-hand  columns  of 
Table  1  might  be  considered  as  "errors”  due  to  instrumentation  in  the  three 
aircraft.  The  maximum  error  is  seen  to  be  0.06  foot. 

As  will  be  seen  more  clearly  later  in  this  report,  the  atmosphere  is  not 
homogeneous  in  the  volume  occupied  by  the  three  aircraft.  This  was  first 
suggested  by  inspection  of  Figs.  1  to  4;  Fig.  4,  perhaps,  indicates  this  most 
clearly.  It  will  be  remembered  that  the  differences  between  aircraft  were 
greater  at  the  lower  altitudes,  where  the  atmosphere  was  more  turbulent,  than 
at  the  higher  altitudes,  where  the  atmosphere  was  smoother.  Thus,  if  it  were 
necessary  to  use  the  three  simultaneous  profiles  in  a  particular  calculation, 
there  is  good  reason  to  use  an  average  of  the  three  profiles  as  the  "true"  or 
"representative”  profile. 

INDIVIDUAL  SOUNDINGS  OVER  MISTRAM  CENTRAL,  100  K.  WEST,  AND 
100  K.  SOUTH  SITES 

Simultaneous  measurement  of  the  refractivity  profiles  over  the  individual 
MISTRAM  stations  is  a  goal  of  the  August  tests.  These  refractivity  profiles 
will  be  used  in  MISTRAM  error  model  calculations  to  see  if  they  result  in 
improved  MISTRAM  accuracy. 

Individual  soundings  were  obtained  over  MISTRAM  Central,  100  K.  West, 
and  100  K.  South,  to  obtain  some  preliminary  data,  as  well  as  to  determine  the 
problems  of  gathering  and  reducing  this  type  of  data. 

In  Figs.  5  and  6,  the  simultaneous  descents  by  the  three  aircraft  over 
MISTRAM  Central,  100  K.  West,  and  100  K.  South,  are  first  shown  separately, 
and  are  then  superimposed  in  the  right-hand  plot  in  each  figure.  These  pro¬ 
files  were  plotted  directly  from  the  raw  data,  with  no  corrections  except  for 
the  scale  correction  to  802  as  explained  earlier. 


22 


4 


23 


Simultaneous  Soundings:  MISTRAM  Central,  100  K.  West,  and  100  K.  South 
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Simultaneous  Soundings:  MISTRAM  Central,  100  K.  West,  and  100  K.  South 


It  is  immediately  apparent  that  the  differences  in  profile  at  the  three 
MISTRAM  stations,  as  shown  on  the  right  in  Figs.  5  and  6,  are  greater  than 
the  differences  in  profile  shown  in  Figs.  1  to  3.  This  was  expected,  of  course, 
since  about  20  miles  separate  MISTRAM  Central  and  100  K.  West  and  100  K. 
South.  If  it  is  assumed  that  the  total  differences  in  profile  in  Figs.  5  and  6 
are  due  to  the  atmosphere  (  and  noi.  to  instrumentation),  then  a  calculation  of 
the  propagation  delay  of  a  radio  wave  through  the  atmosphere,  described  by 
these  refractivity  profiles,  is  information  pertinent  to  the  determination  of 
MISTRAM  bias  errors. 

Table  2  tabulates  the  delays  for  vertical  propagation  of  a  radio  wave  to 
an  altitude  of  15,  000  feet  as  calculated  from  the  refractivity  profiles  shown  in 
Figs.  5  and  6.  As  mentioned  earlier,  the  delay  to  the  same  altitude  for  other 
elevation  angles  may  be  determined  by  multiplying  the  values  in  Table  2  by  the 
cosecant  of  the  elevation  angle. 

If  the  profile  at  MISTRAM  Central  (Fig.  5  or  6)  had  been  assumed  to  be 
present  over  all  of  MISTRAM,  and  if  it  had  further  been  assumed  that  the 
beams  from  each  station  operate  principally  in  the  profile  taken  over  that 
station,  then  the  delays  in  Table  2  (rows  4,  5,  and  6),  appropriately  modified 
for  elevation  angle,  indicate  the  range-difference  errors  that  would  be  produced 
by  the  first  15,  000  feet  of  the  atmosphere.  That  these  differences  are  large  in 
relation  to  MISTRAM  expectations  points  out  the  seriousness  of  the  meteorologi¬ 
cal  effects  upon  MISTRAM. 

FORMATION  LEVEL  RUNS 

The  three  aircraft  made  a  total  of  seven  level  runs  in  a  tight  V-formation 
at  2000  feet  from  various  directions  into  MISTRAM  Central  during  the  two  days 
of  tests.  Two  runs  from  16  June  and  one  run  from  17  June  have  been  converted 
to  common  scales  and  plotted  in  Figs.  7,  8,  and  9. 
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Vertical  Propagation  Delays  of  a  Radio  Wave  to  a  15,  000-Foot  Altitude 
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Formation  Level  Runs,  Three  Aircraft,  2000  Feet 
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Formation  Level  Runs,  Three  Aircraft,  2000  Feet 
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Formation  Level  Runs,  Three  Aircraft,  2000  Feet 


These  plots  were  obtained  by  using  the  Gerber  Model  GDDRS-3B  Digital  Data 
Reduction  System,  modified  to  drive  a  plane  x  -  y  plotter.  The  802  record  com¬ 
pressed  80  N-units  into  about  1.  5  inches  of  curvilinear  deflection  while  the  794 
record  had  a  scale  of  25  N-units  per  2.  5  inches.  The  133  record,  made  with  a 
100  pen  recorder,  was  difficult  to  read  in  fine  detail,  and  this  fact  is  evidenced 
in  the  generally  smoother  character  of  the  133  plots  in  Figs.  7,8,  and  9.  The 
accurate  reduction  of  these  records  to  a  common  scale  is  difficult. 

The  records  from  the  three  aircraft  appear  to  agree  quite  well  in  gross 
detail  on  all  three  runs.  No  conclusions  should  be  drawn  from  the  fine  detail, 
since  the  data  reduction  process  depends  heavily  upon  the  interpretation  and 
skill  of  the  GDDRS  operator. 

These  records  also  give  a  crude  indication  of  the  lack  of  homogeneity  of 
refractivity  in  the  volume  occupied  by  the  three  aircraft.  The  run  in  Fig.  7 
was  about  9.  3  nautical  miles  (n.  m.  )  in  length.  Thus,  the  vertical  lines  in  Fig.  7 
represent  about  2500  feet  of  aircraft  travel.  If  we  assume  a  250-foot  separation 
between  802  and  794,  the  two  outside  aircraft  in  the  V-formation,  then  the  varia¬ 
tions  of  refractive  index  shown  in  Fig.  7,  for  distances  corresponding  to  one-tenth 
the  distance  between  vertical  lines,  might  be  present  between  794  and  802  in 
the  2000-foot  altitude  ranges  during  their  descents  and  ascents. 

It  was  intended  to  compare  this  data  with  the  wiresonde  data  obtained  at 
MISTRAM  Central  during  the  time  the  wiresonde  was  flown  at  a  constant  altitude. 

If  it  is  assumed  that,  on  some  time  scale,  temporal  variations  of  refractivity 
are  spatial  variations  drifting  with  the  wind,  this  process  should  allow  the  sepa¬ 
ration  of  temporal  and  spatial  variations  which  are  combined  in  the  aircraft  data. 

It  should  also  produce  an  absolute  calibration  point  for  the  aircraft  refractometers 
when  they  fly  by  the  wiresonde  balloon. 
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Since  the  winds  at  2000  feet  were  light  and  variable  on  both  operating 
days,  the  first  goal  was  not  attained.  The  wiresonde  value  of  refractive  index 
at  2000  feet  determined  the  vertical  scales  of  Fig.  7.  The  comparable  wire¬ 
sonde  value  for  Figs.  8  and  9  was  not  obtained  because  of  instrumentation 
difficulties.  Accordingly,  on  these  figures,  the  vertical  scales  were  determined 
from  the  794  data. 
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SCALE  AND  TYPE  OF  ATMOSPHERIC  REFRACTIVE  INDEX  ANOMALIES 


R.  M.  Cunningham* 


We  outlined  the  Air  Force  Cambridge  Research  Laboratory's  part  in  ESD's 
refractive  index  project  last  November.  **  D.  L.  Ringwalf's  paper  describes  the 
problems  involved  with  the  accurate  measurement  of  refractive  index  profiles 
and  bias  corrections.  I  will  illustrate  the  scale  and  type  of  atmospheric  refrac¬ 
tive  index  anomalies  we  will  be  talking  about,  from  those  that  appear  as  bias  in 
the  present  interferometer  systems  to  those  that  produce  noise.  J.  H.  Meyer 
and  R.  K.  Crane,  in  their  presentations,  will  show  how  this  type  of  data,  along 
with  aircraft  and  ground  measurements,  can  be  quantitized  in  terms  of  numbers 
familiar  to  those  handling  data  from  long  baseline  radio  interferometer  systems. 

The  atmospheric  anomalies  I  will  be  talking  about  are  those  made  visible 
by  clouds  (Fig.  1).  A  cloud  of  the  cumulus  family  often  develops  in  regions  of 
relatively  dry  air.  We  see  the  common  result  in  terms  of  the  refractive  index. 
This  type  of  trace  is  commonly  obtained  in  flight  through  the  upper-half  of  a 
growing  cumulus  cloud.  A  map  showing  the  distribution  of  cumulus  clouds, 
therefore,  is,  to  a  large  extent,  a  map  of  this  scale  of  refractive  index  anomaly. 

I  will  show  a  number  of  figures  pertaining  to  7  August  1963.  In  Fig.  2, 
we  see  the  peninsula  of  Florida  from  the  Tyros  satellite.  The  clouds  shown, 
which  are  largely  cumulus  cloud  groups ,  can  be  considered  as  areas  of  high 
refractive  index  or  high  bias  to  an  interferometer  system.  Note  the  sharp 
differences  along  the  coast  lines, as  well  as  the  effects  of  lakes,  and  of  the 
gulf  stream.  Note  also  the  occurrence  of  cloud  patches  and  cloud  lines. 

*Cloud  Physics  Branch,  Cambridge  Research  Laboratory,  Bedford,  Massachusetts 

**See  ESD-TDR-64-109,  Vol.  I,  p.  1-15. 
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Fig.  1.  Atmospheric  Anomalies  Through  a  Cloud  Section 


Fig.  2.  Florida  Peninsula  Cloud  Structure  As  Seen  From  A  Tyros  Satellite 
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Interferometer  systems  with  baselines  of  10  miles  and  up  are  obviously  affected 
by  such  atmospheric  structures.  This  scale  of  cloud  distribution  and  its  clima¬ 
tology  should  be  considered  before  locating  any  future  large  baseline  system. 

Descending  from  the  satellite  level  of  100  miles  to  the  U-2  level,  we  can 
examine  the  same  sea  breeze  clouds  over  MISTRAM  in  more  detail.  In  Fig.  3, 
we  can  see  the  MISTRAM  short  legs;  the  left  cross  denotes  the  MISTRAM  Central 
(MC)  site,  and  the  right  cross  the  100  K.  West  (KW)  site.  Various  cloud  sizes 
and  cloud  groups  are  visible.  The  clouds  inland  are  moving  with  the  wind,  while 
the  semi-permanent  cloud  "front"  remains  almost  stationary  at  the  sea  breeze 
convergence  line  shown  by  the  surface  wind  pattern  of  Fig.  4. 


Fig.  3.  MISTRAM  Site  Sea-Breeze  Clouds  as  Seen  in  a  U-2  Aircraft  (7  August  1963; 
10:31  EST) 
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Typical  MISTRAM  Site  Surface  Wind  Pattern  (7  August  1963;  10:30  EST) 


A  more  familiar  type  of  view  of  this  cloud  line  is  shown  in  Figs.  5  and  6. 
MISTRAM  Central  and  the  10  KW  and  10  K.  South  (KS)  are  just  visible  in  the  upper 
center  of  the  picture  of  Fig.  5.  The  vertical  dark  band  covering  part  of  the 
clouds  over  cape  Kennedy  is  a  propeller  blade  of  the  aircraft.  Note  that  Cape 
Kennedy  and  Merritt  Island  are  just  large  enough  to  produce  fairly  large  single 
cloud  masses.  Figure  6  is  a  ground  observers'  view  of  the  same  cloud  line. 

This  picture  was  taken  looking  north  from  the  airport  visible  in  the  lower  center 
of  Fig.  5.  Figure  7,  which  is  one  scale  size  smaller  than  the  mesocale,  shows 
a  close-up  of  one  of  the  sea  breeze  clouds  pictured  previously.  The  scale  of 
refractive  index,  irregularly  illustrated  in  the  figure,  is  given  by  the  roughness 
of  the  cloud  edge.  Turbulent  eddies  in  the  cloud  show  at  the  edge  as  round 
bumps  of  the  order  of  500  feet  or  so  in  size.  This  is  the  scale  which  has  been 
treated  and  described  most  successfully  by  statistical  methods,  by  the  micro¬ 
meteorologist  and  by  those  studying  radio  propagation.  At  the  moment,  we  will 
not  include  this  scale  in  the  papers  to  follow  (J.  H.  Meyer,  R.  K.  Crane), 
although  it  may  be  included  in  a  statistical  fashion  at  a  later  date. 

We  have  looked  at  the  mesoscale  of  cloud  and  associated  refractive  index 
for  this  typical  summer  day  over  MISTRAM;  now  let  us  look  at  the  time  scale  of 
these  features.  Before  we  look,  however,  we  should  note  the  two  following 
definitions:  a  bias  error,  as  the  term  is  used  in  connection  with  MISTRAM 
errors,  can  be  considered  a  constant  error  over  the  time  MISTRAM  is  tracking 
a  missile,  say  5  to  6  minutes,  while  a  noise  might  be  considered  a  fluctuating 
error  at  any  period  less  than  5  to  6  minutes. 

Figures  8  and  9  show  two  U-2  mosaics,  photographed  19  minutes  apart. 
Figure  10  shows  a  history  of  individual  cloud  turrets  measured  from  time-lapse 
ground  photography.  The  ground  camera  was  placed  on  the  outer  beach  SE  of 
MISTRAM  Central.  The  cloud  histories  shown  on  the  top  graph  in  solid  line 
belong  to  the  cloud  group  along  the  sea  breeze  line  at  the  bottom  of  the  two  U-2 
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Fig.  5.  Florida  Coast  Cloud  Line  Observed  from  Ground 


Fig.  6.  Florida  Coast  Cloud  Line  Observed  from  Aircraft 
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Fig.  7.  F  ine  Structure  of  Sea  Breeze  Clouds 


derived  mosaics.  The  dashed  line  represents  the  early  history  of  the  developing 
group  just  below  the  midpoint  of  the  long  E-W  leg  at  MISTRAM.  The  bottom 
graph  depicts  the  history  of  the  large  cloud  group  5  miles  west  of  MISTRAM 
Central.  (A  map  scale  can  be  derived  from  the  visible  two  mile  short  baseline 
legs  of  MISTRAM.  The  main  roads  are  on  a  N-S,  W-E  grid,  so  that  MISTRAM 
legs  are  actually  close  to  WSW  and  SSE. )  It  appears,  when  one  compares  the 
two  successive  U-2  mosaics,  that  this  cloud  group  collapsed  as  a  whole.  Its 
detailed  history  shows,  however,  that  the  individual  turrets  making  up  the  group 
were  less  and  less  successful  in  reaching  the  size  of  their  previous  turrets. 

The  time  scale  and  space  scale  of  cloud  dimension  and  history  are  posi¬ 
tively  correlated.  In  this  case,  for  instance  cloud  turrets  have  dimensions  of 
about  1  mile  and  life  histories  above  the  more  general  lower  cloud  masses  of 
4  to  7  minutes,  while  major  cloud  groups  have  dimensions  of  about  5  miles  and 
life  histories  of  1/2  to  1  hour  or  more. 
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Fig.  8.  MISTRAM  Site  Sea  Breeze  Cloud  Pattern  (7  August  1963;  10:49  EST) 


Fig.  9.  MISTRAM  Site  Sea  Breeze  Cloud  Pattern  (7  August;  11:08  EST) 
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Fig.  10.  Cloud  Growth  History 

In  summary,  then:  Clouds  of  the  cumulus  variety  allow  us  to  see  various 
scales  of  refractive  index  anomaly.  This  fact  permits  the  meteorologist  skilled 
in  mesoscale  weather  phenomena  to  contribute  a  nonstatistical  partial  answer  to 
the  radio  propagation  problems  limiting  the  accuracy  of  radio  interferometer  sys¬ 
tems.  At  a  scale  of  a  mile  or  so,  the  time  rate  of  change  of  these  cloud  features 
is  slow  enough  that  a  bias  and  a  very  low  frequency  noise  correction  could  be  made 
to  the  MISTRAM  data  in  post-fli^t  analysis.  This  would  be  done  with  the  aid  of 
cloud  mosaics  (on  days  with  the  proper  cloud  types)  for  the  time  of  a  shot  and  a 
sample  of  cloud-to-clear-air  refractive  index  taken  shortly  before  or  after  the 
same  shot.  A  real-time  correction  could  also  be  estimated  from  cloud  observa¬ 
tion  and  previous  cloud-clear  air  AN  measurements.  A  cloud  climatology 
related  to  geographic  features  would  help  in  locating  future  systems  to  minimize 
large  bias  and  large  low  frequency  noise  errors  of  the  above  type. 
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DIGITAL  ATMOSPHERIC  PROFILE  GENERATION 


J.  H.  Meyer* 

Dr.  R.  M.  Cunningham  discussed,  through  the  use  of  cloud  photographs 
taken  at  different  altitudes,  some  of  the  various  atmospheric  scale  sizes  which 
exist  in  the  troposphere.  In  this  paper,  the  mesoscale  meteorology  along  the 
E-W  leg  of  the  Valkaria,  Florida,  MISTRAM  site  and  the  method  used  to 
describe  a  particular  realization  of  the  atmosphere  at  a  given  instant  of  time 
are  discussed. 

At  the  second  refractivity  conference,  held  at  The  MITRE  Corporation  in 
November  1963,  Dr.  Cunningham  discussed  some  of  the  refractive  index  data 
gathered  in  Florida  the  previous  August.  The  analysis  of  the  data  presented 
by  Dr.  Cunningham  was  in  the  form  shown  in  Fig.  1.  Since  then,  further 
analysis  of  the  data  has  been  completed.  However,  before  proceeding,  a  review 
of  the  refractive  index  data  and  cloud  cross  section  shown  in  Fig.  1  is  in  order. 

Level-run  airborne  refractive  index  measurements  versus  altitude,  taken 
at  different  altitudes  along  a  flight  path  over  the  E-W  leg  of  MISTRAM,  are 
shown.  The  flight  path  extended  eastward  several  miles  over  the  ocean  and 
approximately  one  to  two  miles  west  of  the  100  K.  West  (KW)  site.  Representa¬ 
tion  of  the  level-run  refractive  index  data  in  the  manner  shown  in  Fig.  1  allows 
one  to  associate  discontinuities  measured  at  a  particular  altitude  with  those 
measured  at  the  other  altitudes.  Each  level  run  is  plotted  versus  time  so  that 
the  distance  between  MISTRAM  Central  (MC)  and  100  KW  decreases  with  altitude. 
This  is  because  the  aircraft  was  flown  at  a  constant  indicated  airspeed.  At  a 
constant  indicated  airspeed,  true  airspeed  or  ground  speed  increases  with 


*Technical  Operations,  Inc.  ,  Burlington,  Massachusetts 
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Level-Run  C-130  Refroctrve  Index  Traces  At-  Various  Altitudes  and  Cloud  Cross  Section  Roughly 
Corresponding  to  Conditions  on  Run  Five  From  A  Synthesis  of  Data  Shown  Above  and  U-2  arxl 
C-lSO  Photographs 


altitude  and,  therefore,  the  distance  between  MC  and  100  KW  decreases  with 
altitude.  The  cloud  cross  section  shown  at  the  bottom  of  Fig.  1  was  constructed 
to  correspond  roughly  to  the  refractive  index  discontinuities  shown  on  Run  No.  5 
from  a  synthesis  of  the  level-run  data  and  photographs  taken  by  the  C-130  and 
U-2  aircraft. 


The  next  step  in  the  analysis  was  a  check  on  the  agreement  between  the 
aircraft  data  and  the  most  representative  radiosonde  data  for  the  time  period 
over  which  the  flight  took  place.  A  mean  was  found  for  a  50-second  data  sample 
over  MISTRAM  Central  for  each  level  run.  Refractive  index  discontinuities 


caused  by  clouds  were  discounted  in  computing  the  mean.  The  mean  values  for 
the  aircraft  level  runs  were  than  compared  with  (Fig.  2)  with  the  radiosonde 
values  for  the  height  at  which  the  level  run  was  made.  The  solid  curve  in  Fig.  2 


is  the  computed  radiosonde  refractive  index  N  profile,  and  the  dashed  curve  is 
the  computed  radiosonde  saturated  refractive  index  (N  )  profile.  The  solid 
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FLT.  244  9  AUGUST  1963 

VALKARIA.  FLA., RADIOSONDE,  10:50  EST. (CURVES  Na  AND  Ns ) 
VALKARIA,  FLA.,  WIRESONDE,  09:45  EST.  (CURVE  WS  ) 
AIRCRAFT  REFRACTIVE  INDEX  VALUES,  10:20  -llill  EST. 

•  MEAN  N-VALUES  OVER  MC  DURING  LEVEL  RUN 
\  *  N-VALUES  FOR  LEVEL  RUN  DURING  PENETRATION  OF  CLOUDS 

A  N-VALUES  DURING  ASCENT  OR  DESCENT  BETWEEN 
E  \  LEVEL  RUNS  (EAST  AND  WEST  OF  MC  ) 
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Fig.  2.  Comparison  of  the  AFCRL  C-130  Refractive  Index  Values  with  the  Most 
Representative  Valkaria,  Fla.,  Radiosonde  N-Profile  for  9  August  1963 
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circles  are  the  50-second  mean  aircraft  values  of  refractive  index  over 
MISTRAM  Central.  The  solid  triangles  are  the  refractive  index  aircraft  values 
sampled  during  ascent  or  descent  between  the  level  runs,  and  are  measurements 
made  east  and  west  of  the  E-W  leg  of  MISTRAM. 

There  is  good  agreement  between  the  aircraft  and  radiosonde  refractive 
index  values,  except  below  2000  feet.  At  those  altitudes,  the  aircraft  values 
have  better  agreement  with  the  wiresonde-computed  refractive  index  values. 

This  difference,  found  between  the  radiosonde,  wiresonde  and  aircraft  data, 
was  present  in  all  the  data  gathered  during  August  1963.  The  exact  cause  for 
this  disagreement  is  unknown,  and  should  be  investigated.  A  further  check  on 
the  agreement  between  the  aircraft  and  radiosonde  refractive  index  data  is 
made  by  comparison  of  the  refractive  index  values  measured  during  aircraft 
penetration  of  growing  clouds  with  the  computed  radiosonde-saturated  refrac¬ 
tive  index .  The  asterisks  in  Fig.  2  are  the  aircraft  saturated  refractive  index 
values.  One  can  easily  see  there  is  good  agreement. 

Using  the  50-second  refractive  index  means  computed  for  each  of  the  level 
runs  over  MISTRAM  Central  for  the  data  shown  in  Fig.  1,  refractive  index 
differences  (An’.s)  were  abstracted  from  the  level-run  index  data.  Combining 
the  An  data  with  surface  An  data  computed  from  meteorological  measure¬ 
ments  taken  at  sites  along  the  E-W  leg  of  MISTRAM  and  using  the  cloud  cross 
section  shown  in  Fig.  1  as  a  model,  the  An  data  was  adjusted  for  time  and 
space  variations.  Iso  An  contours  were  then  drawn  for  the  adjusted  data. 

The  AN  cross  section  and  cloud  cross  section  for  9  August  1963  are  shown 
in  Fig.  3.  The  AN  cross  section  shown  in  Fig.  3  is  a  particular  realization  of 
the  atmosphere  at  a  given  instant  of  time,  which  allows  one  to  use  ray-tracing 
techniques  to  determine  the  effects  of  such  an  atmosphere  on  the  propagation  of 
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Cloud  Cross  Section  for  9  August  1963,  Flight  244,  Along  the  E-W  Leg  of  MISTRAM  with  the 
Corresponding  N-Cro$s  Section, 


electromagnetic  waves  through  it.  Some  of  the  more  significant  features  of 
the  An  cross  section  shown  in  Fig.  3  are: 

(a)  the  shallow  moist  layer  of  air,  as.sociated  with  the  sea  breeze, 
near  the  surface  at  approximately  15  to  20  kilofeet  east  of 
MISTRAM  Central; 

(b)  the  clouds  associated  with  the  sea  breeze  front  at  approxiamately 
10  to  15  kilofeet  east  of  MISTRAM  Central; 

(c)  the  large  changes  in  AN  as  one  goes  from  ambient  air  into  a 
cloud; 

(d)  the  difference  in  the  magnitude  of  An  for  a  growing  cloud  and 
one  that  is  dissipating  or  evaporating  (some  differences  in 
magnitude  can  range  as  high  as  30  N-units);  and 

(e)  the  general  change  of  the  ambient  An  values  from  plus  to 
minus  as  one  follows  the  cross  section  from  east  to  west. 

With  respect  to  the  features  shown  in  the  An  cross  section  for  9  August 
1963  (Fig.  3),  it  is  obvious  that  the  three  tracking  sites,  MISTRAM  Central, 

10  KW  and  100  KW,  would  experience  different  propagation  conditions.  In  this 
particular  case,  MISTRAM  Central  tracking  at  very  low  angles  would  probably 
experience  the  effects  of  the  shallow  moist  layer  to  the  east. 

Airborne  refractive  index  measurements  and  meteorological  measurements 
at  the  ground  sites  along  the  E-W  leg  of  MISTRAM  taken  on  the  7^^  and  8^^  of 
August  1963  were  analyzed  in  the  same  manner.  Figures  4  and  5  show,  respec¬ 
tively,  the  level-run  refractive  index  data  versus  altitude  and  the  associated 
cloud  cross  sections  for  7  August  and  8  August.  Figures  6  and  7  show  the 
comparison  of  the  level-run,  50-second  means  and  other  aircraft  refractive 
index  values  with  this  computed  radiosonde  N-profiles,  respectively,  for 
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Level  Run  C-130  Refractive  Index  Traces  at  Various  Altitudes  and  Cloud  Cross  Section  Roughly 
Corresponding  to  Conditions  on  Runs  4  and  5  from  a  Synthesis  of  Data  Shown  Above  and  U-2 
and  C-130  Photographs 
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Level-Run  C-130  Refractive  Index  Traces  at  Various  Altitudes  and  Cloud  Cross  Section  Roughly 
Corresponding  to  Conditions  on  Runs  5  and  6  from  a  Synthesis  of  Data  Shown  Above  and  U-2 
and  C-130  Photographs 
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FLT.  242 
7  AUGUST  1963 

VALKARIA.  FLA..  RADIOSONDE. 10:35  ESI  (CURVES  N  AND  N  ) 
AIRCRAFT  REFRACTIVE  INDEX  VALUES.  09:36-10:24  EST. 

•  MEAN  N-VALUES  OVER  MC  DURING  LEVEL  RUN 
\  n  N-VALUES  FOR  LEVEL  RUN  DURING  PENETRATION  OF  CLOUDS- 
^  N-VALUES  DURING  ASCENT  OR  DESCENT  BETWEEN 
LEVEL  RUNS  (EAST  AND  WEST  OF  MC  ) 
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Fig.  6.  Comparison  of  the  AFCRL  C- 130  Refractive  Index  Values  with  the  Most 
Representative  Volkorio,  Flo.,  Radiosonde  N-Profile  for  7  August  1963 
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FLT.  243 
8  AUGUST  1963 

VALKARIA,  FLA.,  RADIOSONDE,  10;40EST.  (CURVES  N  AND  N  ) 
VALKARIA,  FLA.,  WIRESONDE,  li;25  EST.(CURVE  WS) 

^  AIRCRAFT  REFRACTIVE  INDEX  VALUES, 10;48 -12:00  EST.  • 

*\  •  MEAN  N -VALUES  OVER  MC  DURING  LEVEL  RUN 

'v  *  N-VALUES  FOR  LEVEL  RUN  DURING  PENETRATION  OF  CLOUDS 
A  N-VALUES  DURING  ASCENT  OR  DESCENT  BETWEEN 
LEVEL  RUNS  (EAST  AND  WEST  OF  MC  ) 
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Fig.  7,  Comparison  of  the  AFCRL  C-130  Refractive  Index  Values  with  the  Most 
Representative  Valkaria ,  Fla.,  Radiosonde  N-Proflle  for  8  August  1963 
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7  August  and  8  August  1963.  Again,  there  is  good  agreement  between  the  air¬ 
craft  and  radiosonde  refractive  index  values,  except  below  the  altitude  of 
2000  feet  on  7  August  and  3000  feet  on  8  August. 

The  cloud  and  AN  cross  sections  for  7  August  and  8  August  1963  are 
shown  in  Figs.  8  and  9,  respectively.  The  significant  features  of  these  two 
AN  cross  sections  are  similar  to  those  mentioned  for  the  9  August  AN  cross 
section;  however,  on  the  7^^  and  8^^  of  August,  the  cross  sections  show  a 
difference  in  the  location  of  the  clouds  associated  with  the  sea  breeze  front. 

Notice  should  be  given  to  the  fact  that  the  shallow  moist  layer  next  to  the 
surface  is  present  in  both  the  AN  cross  sections;  however,  it  extends  inland  a 
iittle  further  than  in  the  case  of  9  August. 

On  7  August,  which  is  the  day  Dr.  Cunningham  described,  the  line  of 
building  cumulus  clouds  along  the  sea  breeze  front,  and  AN  pattern  on  the 
cross  section  show  some  of  the  effects  of  the  topographical  features  along  the 
E-W  leg  of  the  MISTRAM  system.  The  rather  disorganized  AN  pattern  associated 
with  the  clouds,  (Fig.  8)  between  30  and  85  kLlofeet  west  of  MISTRAM  Central  re¬ 
flects  the  effects  of  the  swamp-type  terrain  over  this  particular  part  of  the 
MISTRAM  system.  Convective  heating,  which  is  producing  the  well  defined 
An  pattern  associated  with  the  clouds  over  the  100  KW  site,  is  not  evident  over 
the  above  mentioned  area.  By  considering  the  propagation  of  electromagnetic 
energy  through  the  atmospheres  shown  in  Figs.  8  and  9,  one  can  easily 
note  thatThe  propagation  conditions  are  different  for  any  one  of  the  tracking  sites 
located  along  the  E-W  leg  of  the  MISTRAM  system. 

The  An  data  for  the  three  days  were  quantized,  using  a  500-foot  grid, 
to  determine  what  effects  the  three  AN  cross  sections  presented  would  have 
on  the  tracking  capability  of  the  MISTRAM  system.  R.  K.  Crane,  of  The  MITRE 
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Corporation  talks  about  the  ray-tracing  program  used  to  compute  the  range, 
range  difference  and  range-rate  errors  caused  by  the  particular  realized 
atmospheres  constructed  for  7,  8,  and  9  August  (see  page  87). 
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Cloud  Cross  Section  for  7  August  1963,  Flight  242,  Along  the  E-W  Leg  of  MISTRAM  with  the 
Corresponding  N-Cross  Section 
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Cloud  Cross  Section  for  8  August  1963,  Flight  243,  along  the  E-W  Leg  of  MISTRAM  with  the 
Corresponding  N-Cross  Section 


RAY- TRACINGS  IN  CLOUD  CROSS  SECTIONS  FOR  A  LONG 
BASELINE  INTERFEROMETER 

R.  K.  Crane* 


INTRODUCTION 

Presentations  by  R.  M.  Cunningham  and  J.  H.  Meyer  (see  pages  33  and 
43  ,  respectively)  describe  the  mesoscale  refractive  index  cross  sections  along 
the  E-W  leg  of  the  Valkaria,  Florida,  MISTRAM  site.  In  this  paper,  the  range, 
range  difference,  range  rate,  and  range  difference  rate  errors,  caused  by  the 
refractive  index  variations  described  by  the  cloud  cross  sections,  are  presented 
for  the  MISTRAM  interferometer  system. 

METHOD  OF  ANALYSIS 

The  variation  of  refractive  index  in  the  troposphere  occurs  over  many 
scale  sizes,  ranging  from  the  turbulent  to  the  synoptic  scales.  Two  ranges  of 
scale  sizes  are  considered  for  this  analysis:  those  that  contribute  to  a  refrac¬ 
tive  index  distribution  that  does  not  vary  over  the  spatial  extent  of  the  cross 
section,  and  those  that  cause  refractive  index  variations  over  distances  ranging 
from  500  feet  to  the  largest  dimension  of  the  cross  section.  The  first  of  these 
scale  size  ranges  contributes  to  the  bias  error;  the  second,  to  the  cross- 
sectional  error.  For  convenience,  these  scale  size  ranges  are  labeled  bias 
and  cross  section. 

The  refraction-induced  errors  are  computed  using  the  methods  of  geo¬ 
metric  optics.  The  scale  sizes  being  considered  imply  that  the  changes  of 
refractive  index,  over  the  dimensions  of  a  wave  length,  are  small,  and  that  the 

*The  MITRE  Corporation,  Bedford,  Massachusetts 
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ray-tracing  equations  are  valid.  The  ray  equations  used  are  simplifications 
of  Hamilton's  Canonical  Equations,  valid  for  a  variation  of  refractive  index 
with  height  above  a  spherical  earth.  The  position  of  a  ray  in  space  is  computed 
using  the  spherically  symmetric  ray  equations.  The  modification  of  the  ray 
position  by  the  variation  of  refractive  index,  as  described  by  the  cross  section, 
is  assumed  to  be  negligible.  The  only  effect  of  the  cross  section  is  to  cause 
small  changes  in  the  phase  velocity  along  the  ray. 

For  ease  of  computation,  the  bias  variation  of  refractive  index,  or  bias 
profile,  is  considered  to  be  composed  of  a  number  of  layers,  each  layer  having 
a  linear  change  in  m,  where  m  is  defined  as  the  product  of  the  index  of  refrac¬ 
tion  and  the  distance  from  the  center  of  the  coordinate  system.  With  this 
assumption,  the  ray  equations  may  be  integrated  directly  in  terms  of  trigono¬ 
metric  functions. 

This  approximate  description  of  the  bias  profile  describes  the  change  of 
refractive  index  with  height  in  almost  the  same  manner  as  the  familiar  linear 
change  in  refractive  index  description,  the  maximum  difference  between  the 
two  approximate  descriptions  being  o.lN  for  the  profiles  considered. 

Using  the  linear  m  description,  the  ray  equations  are: 


6  =  +  A  61 ; 


m 


nr 


=  (l  +  Nx  10“^)(Ro  +  h)  ; 


h 


^Subscript  i  indicates  value  at  the  lowest  height  in  the  layer. 
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The  terms  B,  h,  Rq  and  a  are  defined  in  Fig.  1. 


Fig.  1.  Geometry  for  Ray  Tracing 
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Using  these  equations,  the  central  angle  for  the  point  of  intersection  of 
the  ray  with  an  arc  of  constant  height  or  the  point  of  intersection  of  the  ray  with 
the  radius  vector  for  a  given  central  angle  may  be  computed.  The  cross  section 
description  of  refractive  index  variation  consists  of  a  table  of  AN  values,  one 
value  for  each  intersection  of  an  arc  of  constant  height,  and  a  radius  vector  as 
shown  in  Fig.  2.  The  value  AN  specifies  the  difference  between  the  measured 
refractive  index  at  a  point  in  space  and  the  value  of  refractive  index  for  that 
point  specified  by  the  bias  profile.  The  electrical  range  along  the  ray,  as 
determined  by  a  transit  time  measured,  is  given  by  Rg  =  CA^  , 
where 

ds 

V 

long  ray 


Fig.  2.  Geometry  for  Ray  in  the  Cross  Section 
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/  bias  \  cross  section 


where 


c  =  velocity  of  light  in  free  space, 

n  =  velocity  of  light  in  the  medium,  and 

s  =  arc  length  along  ray. 

The  range  error  produced  by  the  cross  section,  AR^^,  is  found  by  determining 
the  line  integral  of  AN  for  the  ray  in  the  cross  section.  The  position  of  the 
ray  is  assumed  to  be  determined  solely  by  the  bias  profile.  The  incremental 
arc  length.  As  between  each  ray  grid  line  intersection  can  be  computed  using 
the  integral  for  As  for  a  linear  m  profile  for  each  layer. 

1  +  r 

As  =  I  CSC  ad  (^1  +  r) 
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where  the  variables  are  as  defined  above. 

The  cross  sectional  range  error  is  calculated  by  assuming  that  the 
equation  may  be  approximated  as: 


^  \\N  As 


where  AN  is  taken  as  one-half  the  sum  of  the  values  of  AN  for  each  end  of 
the  incremental  segment  of  arc ,  AS  . 

The  value  of  AN  for  each  grid  line  ray  intersection  is  calculated  by  linearly 
interpolating  between  values  of  AN  for  adjacent  grid  lines. 

The  bias  range  error  ARq  is  calculated  by  subtracting  the  true  range  to 
the  target  from  the  electrical  range.  The  equation  used  to  compute  Rg 
is  given  by: 
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where  A0.  =  -0.  and  the  remaining  variables  are  as  defined  above 

The  range  error  for  a  given  ray  is  computed  as  a  function  of  the  initial 
elevation  angle  and  the  ray  end-point  heights.  The  initial  elevation  angles  for 
two  rays  that  intersect  at  the  true  target  pyosition  must  be  calculated,  see  Fig. 

3  in  order  to  determine  the  range  difference  error  for  an  interferometer 
measurement.  The  method  of  calculation  used  was  to  assume  a  value  of  initial 
elevation  angle  calculated  6  for  the  end-pbint  height,  and  compare  the 
results  with  the  given  target  position.  The  initial  elevation  angle  was  then 
corrected  using  the  difference  between  the  result  and  the  given  value,  as: 

^0  (^c  -  ^r)  «0A 

Ad  =  - 

S 

This  process  was  iterated  until  the  difference  between  6^  and  was 

_9 

less  than  5.  0  x  10  milliradians.  The  computations  were  performed  using  an 
IBM  7030  Stretch  computer,  and  the  resultant  accuracy  of  the  refraction  effects 
for  a  given  linear  profile  is  shown  in  Table  1.  These  accuracies  relate  only  to 
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computations  based  on  the  bias  approximation  to  the  refractive  index  variations 
in  space.  They  do  not  include  the  uncertainties  associated  with  using  the  bias 
profile  to  model  the  atmosphere. 


Table  1 

Computation  Accuracy  for  a  Given 
Linear  m  Profile 
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CALCULATED  RANGE  ERRORS 


The  refractive  index  cross-sections,  prepared  by  J.  Meyer,  describe  the 
N-variations  over  a  horizontal  distance  of  150,  000  feet  and  a  vertical  distance 
of  32,  000  feet.  The  deviation  of  refractive  index  from  the  bias  profile  value  is 
assumed  to  be  zero  outside  of  the  given  cross  section.  The  cross-sectional 
errors  were  calculated  for  elevation  angles  between  10  and  35  degrees  for  a 
target  at  32,  808  feet  (10  kilometers).  These  errors  apply  for  all  heights  above 
32,808  feet  for  the  indicated  elevation  angles,  because  the  cross  section  has 
zero  values  above  this  height.  The  cross-sectional  range  errors  plotted,  for 
three  different  antenna  locations  in  the  cross  section,  are  shown  in  the  figures 
on  pages  68  and  72.  These  ray  path  starting  points  correspond  to  the 
location  of  the  three  MISTRAM  sites  along  the  E-W  leg.  The  interferometer 
calculations  were  performed  using  the  assumption  that  the  target  trajectory 
lies  in  the  plane  of  the  three  sites  or  in  the  plane  for  which  the  cross  section 
is  valid. 

7  August  1963 

The  bias  profile,  rawinsonde  profile,  and  aircraft  measurements  of 
refractive  index  for  11:  00  EST  for  7  August  1963  are  shown  in  Fig.  4.  The 
bias  profile  was  generated,  using  a  best  fit  approximation,  to  the  data  points 
below  20, 000  feet,  the  rawinsonde  data  from  20,  000  to  40,  000  feet,  and  the 
lACO  standard  atmosphere  (dry  term  only)  above  40,  000  feet.  The  cross 
section  representative  of  this  particular  time  is  pictured  in  Fig.  6.  The  cross 
sectional  range  errors  for  each  of  the  MISTRAM  sites  are  given  in  Fig.  5. 

One  percent  of  the  bias  range  error  for  a  target  at  32,  808  feet  is  plotted  for 
reference. 


67 


OD 


Fig.  4.  Refractive  Index  Profile, 
7  August  1963,  Valkaria, 
Florida 
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Fig.  5.  Cross-Sectional  Range  Error,  7  August  1963,  MISTRAM  Central,  Valkaria, 
Florida 
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AN  and  Cloud  Cross  Sections,  7  August  1963,  Valkaria,  Florida 


The  rays  from  both  the  MISTRAM  Central  and  the  10  K.  West  (KW)  sites 
lie  in  the  clear  (cloud-free)  region  over  the  range  of  elevation  angles  shown.  At 
elevation  angles  between  10  to  20  degrees,  the  10  KW  ray  is  sweeping  up  into 
region  of  space  with  AN  =  +10  ;  this  region  contributes  to  the  peak  in  the  range 
error  curve  at  an  elevation  angle  of  about  10  degrees.  The  100  KW  ray  passes 
through  the  sea  breeze -induced  clouds  at  about  25,  000  feet  west  of  MISTRAM 
Central.  As  the  ray  sweeps  up,  it  gets  above  the  cloud  with  the  high  AN-values, 
and  the  cross-sectional  range  error  drops. 

8  August  1963 

The  bias  profile,  AN  cross  section  and  cross  sectional  errors  are  shown 
in  Figs.  7,  8,  and  9.  The  MISTRAM  Central  rays,  for  elevation  angles  less 
than  18  degrees,  and  the  10  KW  rays,  for  elevation  angles  less  than  13  degrees, 
pass  under  the  sea  breeze-induced  clouds.  The  cross-sectional  range  errors 
increase  for  higher  elevation  angles  as  the  ray  passes  through  the  clouds. 
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Fig.  7.  Refractive  Index  Profile, 
8  August  1963,  Valkaria, 
Florida 
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AN  and  Cloud  Cross  Sections,  8  August  1963,  Valkaria,  Florida 


Fig.  9.  Cross-Sectional  Range  Error,  8  August  1963,  MISTRAM  Central, 
Volkorlo,  Florida 
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The  100  KW  rays  pass  through  the  clouds  about  70,  000  feet  west  of  the  Central 
site  for  elevation  angles  near  25  degrees.  The  result  is  the  rapid  change  in 
cross-sectional  range  error  for  increasing  elevation  angle. 

9  August  1963 

The  bias  profile,  AN  cross  section,  and  cross-sectional  range  errors 
are  given  in  Figs.  10,  11,  and  12.  The  10  KW  rays  pass  through  the  sea  breeze 
clouds  for  elevation  angles  greater  than  10  degrees.  For  elevation  angles  above 
20  degrees,  the  ray  passes  through  a  region  of  low  AN  in  the  center  of  the 
clouds.  The  sea  breeze  clouds  are  closer  to  the  beach  than  they  were  on  either 
the  seventh  or  eighth  of  August.  This  results  in  an  increase  in  the  MISTRAM 
Central  cross-sectional  errors  at  lower  elevation  angles  than  on  the  other 
two  days. 


Fig.  10.  Refractive  Index  Profile, 
9  August,  1963,  Valkaria 
Florida 
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II.AN  and  Cloud  Cross  Sections,  9  August  1963,  Valkaria,  Florida 


Fig.  12.  Cross-Sectional  Range  Error,  9  August  1963,  MISTRAM  Central 
Valkaria,  Florida 
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Interferometer  Errors  for  7  August  1963 


The  preceding  data  gives  the  cross-sectional  range  error  as  a  function  of 
elevation  angle  for  a  single  ray.  For  an  interferometer  measurement,  two  rays 
are  combined,  each  at  a  different  elevation  angle.  The  correct  elevation  angles 
for  rays  intersecting  at  both  100, 000  and  999,  999  feet  were  computed  to  deter¬ 
mine  the  errors  for  interferometer  measurements.  The  results  are  given  in 
Figs.  13  to  22.  The  elevation  angle  error  is  defined  in  Fig.  3.  The  baseline 
comparison  value, 

1 

Sp  =  P,  -  —  P,  , 

^  10 

is  given  in  Figs.  22;  this  Figure  also  compares  the  range  difference  errors 
(Sp)  for  each  of  the  three  days  for  which  data  was  processed. 

On  7  August,  it  can  be  seen  that  the  cross-sectional  range  error  meets 
specifications  for  angles  above  10  degrees,  while  the  range  difference  error 
for  the  short  baseline  does  not.  The  comparison  of  errors  between  the  two 
baselines  does  not  conclusively  show  that  a  long  baseline  system  is  better  than 
a  short  one.  The  actual  errors  strongly  depend  upon  target  position.  The 
range  difference  rate  error  alone  does  not  meet  specifications  for  the  long 
baseline  at  angles  less  than  13  degrees.  The  short  baseline  range  rate  error 
meets  specifications  over  the  range  of  values  of  elevation  angle  that  have  been 
plotted.  These  rate  errors  are  strongly  affected  by  the  target  trajectory  and 
velocity;  the  assumed  rate  of  change  of  initial  elevation  angles  at  the  100  KW 
site  was  1  milliradian/second. 

An  actual  target  trajectory  would  have  to  be  used  to  compute  the  range, 
range  difference,  range  rate,  and  range  rate  difference  errors  for  comparison 
with  Cape  Kennedy  data.  The  results  discussed  in  this  report  pertain  to  a  simple 
constant-height  trajectory. 
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It  should  be  noted  that  the  range  difference  rate  errors  for  actual  missile 
trajectories  may  differ  from  those  indicated  above.  The  purpose  of  using  a 
simple  trajectory  was  to  enable  us  to  present  an  example  of  how  the  profile  can 
be  used  to  study  a  problem  area  of  an  operational  system.  Now  that  the  tech¬ 
nique  has  been  developed  and  the  necessary  computer  program  written,  it  is 
certainly  a  simple  matter  to  continue  using  this  approach  for  other  similar 
investigations  on  the  MISTRAM  system. 
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Fig,  13.  MISTRAM  System  Profile  and  Cross  Section,  7  August  1963;  Target  Height: 
999,999  Feet 


78 


Fig.  14.  MISTRAM  System  Profile  and  Cross  Section,  7  August  1963;  Target  Height: 
1 00,000  Feet 
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Fig.  15.  Range  Rate  Error  and  Range  Difference  Rate  Errors  Due  to  Cross  Section, 
7  August  1963,  MISTRAM  System;  Target  Height:  999,999  Feet 
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Fig.  16.  Range  Rate  Error  and  Range  Difference  Rate  Errors  Due  to  Cross  Section, 
7  August  1963,  MISTRAM  System;  Target  Height:  100,000  Feet 
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Fig.  17.  Range  Difference  Rate  Error,  7  August  1963;  Target  Height:  999,999 
Feet,  Baseline:  100,000  Feet 
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Fig.  18.  Range  Difference  Rate  Error,  7  August  1963;  Target  Height:  100,000 
Feet,  Baseline:  100,000  Feet 
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Fig.  19.  Range  Difference  Rate  Error,  7  August  1963;  Target  Height:  100,000  Feet, 
Base  li ne:  10,000  Feet 
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Fig.  20.  Elevation  Angle  Errors,  Profile  and  Cross  Section,  7  August  1963;  Target 
Height;  999,999  Feet 
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Fig.  21 .  Elevation  Angle  Errors,  Profile  and  Cross  Section,  7  August  1963;  Target 
Height:  100,000  Feet 
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Fig.  22.  Range  Difference  Errors  as  a  Function  of  Elevation  Angle 
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ESDMAUI,  HAWAH,  MEASUREMENTS 
M.  C.  Thompson,  Jr.  * 


INTRODUCTION 

The  installation  of  the  9.  4-gc  transmitter  and  receiver  for  the  X-band,  fixed 
end-point  phase  measurements  at  Maui,  Hawaii,  was  started  in  November  1963; 
the  first  data  run  was  performed  in  January  1964,  and  the  field  measurements 
were  scheduled  for  completion  in  mid- July  1964.  Current  plans  are  to  return 
the  equipment  to  Boulder,  Colorado. 

The  following  outline  summarizes  the  various  test  conditions  for  which  data 
was  obtained.  The  location  and  the  profile  of  the  propagation  path  on  Maui  are 
given  in  Vol.  1  of  the  previous  Proceedings.  **  The  location  of  the  instrumenta¬ 
tion  on  the  airstrip  is  shown  in  Fig.  1. 


*  National  Bureau  of  Standards,  Boulder,  Colorado 

♦♦See  ESD-TDR-64-103,  Vol.  1,  p.  129,  Figs.  l-8(a)  and  l-8(b). 
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MEASUREMENT  OUTLINE 


Antenna  Aperture  Smoothing 

Recorded  range  variations  simultaneously,  over  paths  using  3-  and  0.  5- 
meter  antennas,  mounted  side-by-side,  during  January  6,  7  and  8,  and  14,  15 
and  16.  (See  Fig.  1-9  (c)  Page  1-130  for  ESD-TDR-64-103. ) 

Comparison  of  Range  Difference  Between  a  Horizontal  In-Line  Baseline  and 

an  Equivalent  Vertical  Baseline 

Recorded  ranges  and  range  differences  simultaneously,  using  a  30- 
meter  vertical  baseline  and  a  corresponding  horizontal  in-line  baseline, 
during  February  11,  12  and  13.  (See  Fig.  l-9d  of  above  reference. ) 

Short  Horizontal  In-Line  Baselines 

Recorded  range  differences  simultaneously,  using  baselines  of  40-, 

100-  and  250-  meters,  on  March  4,  5,  and  6  (see  Fig.  2). 
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Fig.  2.  Horizontal  Short  In-Line 
Antenna  Arrangement 


Long  Horizontal  In-Line  Baselines 


Recorded  range  differences  using  baselines  of  250-,  400-,  855-  and 
1255-meters,  from  March  18  to  24,  and  from  April  10  to  17  (see  Fig.  3). 


Fig.  3.  Horizontal  Long  In-Line 
Antenna  Arrangement 


Short  Azimuth  Baselines 

Recorded  range  differences  simultaneously,  using  baselines  of  4.  5-, 
31.  5-  and  201.  5-  meters,  normal  to  the  propagation  path,  from  May  7  to  13 
(see  Fig.  4). 

Long  Azimuth  Baselines 

Recorded  range  differences  simultaneously,  using  baselines  of  201.  5-, 
598.  5-,  800-  and  1700-meters,  from  May  21  to  28  (see  Fig.  5). 
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Fig.  4.  Antenna  Arrangement  for  Short 
Azimuth  Baseline 


O 

Fig.  5.  Antenna  Arrangement  for  Long 
Azimuth  Baseline 
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Comparison  of  Surface  and  Airborne  Index  Measurements  for  Range  Corrections 


Recorded  one  range  and  surface  index  at  six  points  distributed  around  the 
lower  antenna.  Recorded  index  along  vicinity  of  path,  using  airborne  refrac- 
tometer,  from  June  10  to  17. 

FUTURE  PLANS 

The  data  from  these  experiments  will  be  analyzed  in  terms  of  the  spectra 
of  the  various  range  and  range  difference  fluctuations.  Both  these  spectra  and 
the  variables  will  be  examined  for  correlation  with  other  parameters,  such  as 
the  index,  for  possible  application  of  the  index  for  prediction  or  correction  of 
range  errors. 
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FIELD  TEST  RESULTS  OF  THE  LINE  INTEGRAL  REFRACTOMETER 


Harold  M.  Richardson* 


INTRODUCTION 

The  molecular  resonances  of  water  vapor  and  oxygen  molecules  have  been 
the  subject  of  theoretical  and  experimental  investigations  because  of  the  practical 
importance  of  the  effect  with  regard  to  microwave  attenuation.  The  theory  pre¬ 
dicts  that  there  will  be  a  slight  dispersion,  or  change  of  refractive  index  with 
frequency,  associated  with  these  resonances.  The  dispersion,  which  can  be 
calculated  from  the  known  physical  constants  of  the  molecules,  is  so  small 
that  it  is  very  difficult  to  measure  with  current  laboratory  techniques. 

These  dispersive  regions  cause  a  differential  in  the  propagation  transit 
time  of  two  signals  at  different  frequencies  as  the  waves  traverse  an  atmospheric 
path.  With  a  suitable  choice  of  frequencies,  a  measurement  of  transit  time 
differentials  provides  information  on  the  total  tropospheric  refractivity,  inte¬ 
grated  over  the  path.  This  measurement  could  provide  a  direct  correction  to 
precision  microwave  distance-measuring  equipment  to  compensate  for  the 
slowing  effects  of  the  clear  air  troposphere,  as  well  as  a  means  of  gaining  some 
insight  into  meteorological  and  propagation  parameters.  ** 

PRINCIPLE  OF  OPERATION 

When  the  scale  of  the  propagation  path  inhomogeneities  permits  the 
approximation  of  geometric  optics,  the  transit  time  (t)  may  be  expressed  as 


*The  MITRE  Corporation,  Bedford,  Massachusetts 

**Background  information  and  prior  testing  results  on  the  MITRE  Line  Integral 
Refractometer  (LIR)  may  be  found  in  ESD-TDR-64-103,  Vol.  I,  pp.  1-77  to  1-106. 
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where 

c  is  the  vacuum  velocity  of  light, 
ds  is  an  elemental  arc  length, 

S  is  the  curved  ray  path  over  which  the  integration  is  performed. 

The  integrated  refractivity  over  the  path  (/g  N  ds^  may  be  recognized  as 
a  measure  of  the  slowing  effect  caused  by  the  presence  of  the  material  atmosphere 
on  the  electromagnetic  waves.  This  quantity  is  the  first-order  correction  to 
inferences  of  distance  from  transit  time  measurements.  Ray  bending  produced 
by  refractivity  gradients  causes  the  arc  length  S  to  differ  from  the  straight-line 
distance.  However,  this  difference  is  a  second-order  effect,  provided  the  angle 
between  the  ray  and  the  refractivity  gradient  does  not  approach  90  degrees ,  and 
can  be  neglected  for  atmospheric  paths  at  elevation  angles  greater  than  a  few 
degrees. 


The  differential  transit  times,  phases,  and  refractivities  for  two  coherent 
frequencies  ,  (2^  >  transmitted  over  the  same  path,  are  related  as  follows: 


where 

•Al  and  1A2  are  the  phase  angles  corresponding  to  the  transit  times,  and 
Dj  2  is  the  dispersivity. 
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This  equation  shows  that  a  special  kind  of  phase  measurement  is  proportional 
to  the  line  integral  of  differential  refractivity  (dispersivity).  The  differential 
phase  is 

[(^2/^1)  'P\  ~  '/'2^  h  • 


FIELD  TESTS 

A  transmitter  and  two  identical  receivers  for  a  measurement  of  water 
vapor  have  been  constructed  and  mounted  in  truck  vans  for  field  tests.  This 
apparatus  (or  portions  thereof)  have  been  operated  over  the  seven  paths  listed 
in  Table  1. 

The  Bald  Knob  to  Mt.  Rowe,  N.  H.  ,  path,  of  0.  2-degree  elevation,  12.  8- 
mile  length,  was  used  for  calibration-testing  of  phase  versus  humidity  during 
October  and  November,  1963.  The  presence  of  strongly  variable  fluctuations* 
in  the  phase  readings  from  this  horizontal  path  led  to  an  acceleration  of  the 
slanted  path  testing  program.  The  Southern  Appalachians  provided  the  next 
four  paths,  elevated  2.  6,  1.  9,  3.  2,  and  4.  9  degrees,  and  with  lengths  as 
indicated  in  the  Table.  The  Toll  House  Gap,  N.  C.  ,  to  Johnson  City,  Tenn.  , 
path  was  slightly  vignetted  in  the  lower  side  of  the  main  lobe  of  the  15.  6-gc 
pattern,  and  was  used  primarily  for  checkout  prior  to  operation  on  the  Toll 
House  Gap  to  Pinnacle  Peak,  Tenn.  ,  path.  The  paths  from  Camp  Creek  Bald 
to  Tusculum  Airport,  Tenn.  ,  and  from  Camp  Creek  Bald  to  Moncier  cow 
pasture,  also  in  Tennessee,  were  used  in  a  beam-switching  experiment  prior 


♦See  ESD-TDR-64-103,  Vol.  I,  Figs.  1-12  and  1-13. 
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to  utilization  of  the  latter  path  for  long-term  (24-hour)  test.  The  two  paths 
were  used  primarily  for  checkout  operations  and  antenna  phase  pattern 
measurements. 


Table  1 
Test  Paths 


Paths 

Elevation 
Angle  (deg) 

Length 

(miles) 

Date 

Remarks 

Bald  Knob  to  Mt  Rowe, 

N.  H. 

0.  2 

12.  8 

Oct.  to 

Nov.  '63 

Phase- 

Humidity 

Comparison 

Toll  House  Gap,  N.  C. 
to  Johnson  City,  Teim. 

2.  6 

18.  3 

Apr.  '64 

Slanted  Path 
(Vignetted) 

Toll  House  Gap,  N.  C. 
to  Pinnacle  Peak,  Tenn. 

1.  9 

15.  3 
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The  path  profiles  are  illustrated  in  Fig.  1.  The  dashed  lines  represent 
the  position  of  the  first  antenna  null  which  clears  all  obstructions  to  the  mid¬ 
point  (from  either  end),  except  for  the  Toll  House  Gap  to  Johnson  City  path. 

The  stability  of  the  apparatus  taken  into  the  field,  as  illustrated  in  Fig,  2, 
is  within  0.  1-degree,  1-second  time  and  1-decibel  changes  in  the  15.  6-  and 
31.  2-gc  carriers. 

TEST  RESULTS 

The  elevated  path  testing  results  indicated  the  presence  of  strongly  variable 
fluctuations.  The  peak-to-peak  fluctuations  (Fig.  3)  were  determined  from  the 
first  2.  8  minutes  of  a  40-minute  test  record,  ultimately  used  to  determine  the 
long-term  power  spectra.  The  data  were  obtained  over  a  6.  5-mile,  4.  9-degree 
elevated  path. 

The  fluctuations,  as  the  figure  shows,  become  small  around  sunset  on 
5  May,  and  then  increase  for  the  period  from  midnight  of  5  May  to  about  sun¬ 
rise  of  6  May.  The  fluctuations  near  sunset  of  6  May  were  again  relatively 
small  before  becoming  alternately  noisy  and  quiet  from  near  midnight  6  May 
to  almost  noon  of  7  May.  These  fluctuations  have  been  observed  to  change 
from  approximately  2  degrees  peak-to-peak  to  10  degrees  peak-to-peak  in 
approximately  5  minutes. 

On  5,  6  and  7  May,  all  the  data  was  obtained  under  atmospheric  conditions 
of  warm,  clear,  stable,  high-pressure  air,  with  calm  to  very  light  surface 
winds.  Variable  direction  15  mile Aour -winds  were  present  at  the  transmitter 
site  on  Camp  Creek  Bald  (elevation  4440  feet).  During  the  test  period,  a 
marked  dip  and  subsequent  return  to  normal  in  water  vapor  content  occurred 
at  the  5000-foot  level,  as  indicated  by  Nashville,  Tenn.  ,  and  Greensboro,  N.  C,  , 
rawinsonde,  but  no  change  was  discernible  ,at  either  van  site. 


99 


Fig.  1.  Refroctometer  Test  Site  Profiles 
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Fig.  2.  Apparatus  Stability 


Power  Spectra 

Ten  power  spectra,  as  determined  from  40-minute  runs,  have  been  com¬ 
puted  for  the  Camp  Creek  Bald  to  the  Moncier  cow  pasture  path;  in  addition, 
two  long-term  power  spectra  have  been  derived  by  piecing  the  serial,  individual 
x'uns  with  some  interpolations.  Five  spectra  have  been  obtained  for  the  Toll 
House  Gap  to  Pinnacle  Peak  path,  two  spectra  for  the  Toll  House  Gap  to  Johnson 
City  path,  and  three  spectra  of  the  apparatus  with  100-foot  paths. 
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Fig.  3.  Variable  Signal  Fluctuatians  for  Slanted  Paths 

Some  typical,  early  SF>ectra  are  shown  in  Fig.  4  for  the  6.  5-mile,  4.  9- 
degree  elevated  path.  They  are  characterized  by  a  relatively  flat  region,  a 
break  point  between  10  ^  and  10^  cycles/second,  and  a  rapid  fall  off  at  the  higher 
frequencies.  The  power  densities  of  the  spectra  vary  considerably  in  the  flat 
region  and  to  a  lesser  degree  in  the  slope  after  the  break  point.  The  characteris¬ 
tic  shape  of  the  spectra  was  similar  to  that  obtained  during  the  Initial  testing 
over  Lake  Winnipesaukee,  N.  H. 

The  11  April,  100 -foot  curve  in  Fig.  4  indicates  the  SF>ectrum  characteristic 
of  the  apparatus.  Autocorrelation  functions  have  been  obtained  for  all  the  spectra, 
with  typical  results,  as  shown  in  Fig.  5,  for  no  atmosphere  (100  feet)  and  with 
atmosphere  (6.  5  miles)  as  indicated. 
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g.  4.  Characteristic  Spectra 


Fig.  5.  Autocorrelation  Curve  of  Signals 
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Figure  6  shows  two  interesting  spectra;  the  low-frequency  power  densities 
are  approximately  equivalent,  but  the  break  point  is  significantly  different.  The 
wind  at  the  receiver  was  essentially  nil  at  both  times  of  data  collection,  but  the 
wind  at  the  transmitter  was  approximately  20  miles/hour  for  the  6  May  spectrum 
and  about  10  miles/hour  for  7  May.  The  path  average  wind  is  unknown,  but 
these  spectra  and  wind  observations  tend  to  indicate  wind  effects. 


Fig.  6.  Comparative  Spectra 


Baseline  Experiments 

A  very  limited  set  of  data  was  taken  over  three  small  baselines,  8.  75,  16.  3 

and  184  feet.  The  two  power  spectra  for  the  8.  75-foot  baseline  are  shown  in 

-2  5  -1 

Fig.  7.  Inspection  shows  common  features  in  the  spectra  from  10  to  10 
cycles/second.  In  spectroscopic  terms,  there  appear  to  be  some  common 
"absorption  lines"  in  the  spectra.  Similar  spectra  have  been  computed  for  the  other 
baselines,  the  cospectra,  coherency  matrix,  etc.  ,  are  currently  being  calculated. 
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Baseline  Spectra  (8.75-Foot  Separation) 


The  analog  records  of  rectified  carrier  amplitude  monitors  for  the  184- 
foot  separation  are  shown  in  Fig.  8.  The  monitors  in  each  receiver  show  the 
usual  degree  of  correlation,  with  little  correlation  apparent  from  receiver  to 
receiver;  Fig.  9  illustrates  significantly  more  correlation  from  receiver  to 
receiver.  Close  examination  of  the  details  of  the  recordings,  illustrated  in 
Fig.  10,  reveal  a  relatively  high  degree  of  correlation  of  all  monitors  but  with 
a  delay  in  time,  for  Receiver  I  relative  to  Receiver  II,  of  approximately 
1  second.  This  delay  agrees  well  with  the  observed  cross  wind  at  the  receiver 
and  the  "frozen-in"  structure  model. 


11  MAY  1964 
6.5-MlLE  PATH 
184-FT  SEPARATION 


RECEIVER  II  " 

31.2  GC-—^ 

RECEIVER  II  ' 

15.6  GC - -  VV 


RECEIVER  I 
15.6  GC-^ 


RECEIVER  b 
31.2  GC-^ 


Idb 
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"•"•Sit. 


Fig.  8.  Carrier  Amplitude 
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U  MAY  1964 

6,5-MiLf  PAIH 

84T.  9- IN.  SEPARATION 

RECEIVER  [l 

31.2GC- 

\^/ 

' 

RECEIVER  1 

‘  ■■■  'i- 

31.2GC-  ,  ' 

A  . 

RECEIVER  1 

15.6GC- 

RECEIVER  11/  — : — 

15.6  GC — ^  f 

■■■--  V  -V  H 

Fig.  9.  Comparison  of  Receiver  Signal 


29  APRIL  1964 
15.3-Mll£  PATH 
16-FT.  4- IN.  SEPARATION 


RECEIVER  II 
15. 6  GC - -- 

RECEIVER  II 
31.2  GC - 


RECEIVER  I 
31. 2  GC'^ 


REaiVER  I  / 
15.6  GC — ^ 


Fig.  10.  Comparison  of  Carrier  Amplitude  Signals 
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Beam-Switching 


The  main  purpose  of  the  expedition  to  Tennessee  was  to  determine  the 
influence  of  elevation  angle  on  the  variable  fluctuations,  over  and  above  atmos¬ 
pheric  changes.  The  continued  presence  of  the  fluctuations  precluded  path-to- 
path  comparison,  except  for  near-simultaneous  or  extensive  observations.  The 
near-simultaneous  observation  was  made  in  the  transmitter  beam-switching 
(manual)  experiment,  with  Receiver  I  located  in  the  cow  pasture  (6.  5-mile  path, 
4.  9-degree  depression)  and  Receiver  II  located  at  Tusculum  airport  (10-mile 
path,  3.  2-degree  depression).  The  results  of  the  experiment  conducted  over  a 
40-minute  period,  are  shown  in  Fig.  11.  The  received  signals  experienced 
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Fig.  11.  Peak-To-Peak  Signal  Variation 
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phase  fluctuations  of  9.  4  degrees,  peak-to-peak,  in  Receiver  II  during  the  first 
three  minutes  of  the  run;  Receiver  I  measured  4.  2-degree  phase  fluctuations 
from  5  to  9  minutes  from  start;  Receiver  II  then  measured  11.  6  degree  fluctua¬ 
tions  from  13  to  17  minutes,  etc.  After  scaling  the  short-path  values  for 
comparison  with  those  of  the  longer  path,  the  effect  of  elevation  was  inversely 
proportional  to  angle,  i.  e.  ,  /^2)  =  (01/^2)  ,  or  12  degree/8  degree 

_  -4.9 

phase  *  elevation  angle.  The  data  is  too  limited  to  devise  a  general 

low  which  describes  the  dependence  of  elevation  angle  on  phase  fluctuations,  but 
it  does  indicate  that  elevation  angle  is  a  parameter  to  be  considered. 

The  following  equations  are  the  basis  for  the  scaling  ratio  and  is  given  for 

3 

measurements  of  N  =  60  N  xmits  for  10  gm/m  water  vapor  and  Nj  -  N2  is 
determined  by  the  second  equation  for  the  same  density: 
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The  microwave  portion  of  the  LIR  is  shown  in  Fig.  12  and  the  digital  phase 
shifter  portion  of  a  complete  LIR  including  the  oxygen  phase  shifter  is  shown  in  Fig. 
13.  The  oxygen  microwave  hardware  for  the  LIR  is  currently  under  development. 


flEFRACTOMCTE/< 

P£C£/VE/t 


Fig.  12.  Microwave  Electronics 
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Fig.  13.  Digital  Phase  Shifter 
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Smoothing  Times 


The  early  testing*  demonstrated  that  the  LIR  responded  to  atmospheric 
changes  when  the  fluctuations  were  smoothened  by  an  averaging  process.  Con¬ 
siderations  of  the  magnitude  of  the  fluctuations  and  the  stability  of  the  atmosphere 
lead  to  the  conclusion  that  the  short-time  fluctuations  are  caused  by  some  propa¬ 
gation  mechanism  not  representative  of  density  or  refractive  variations.  This 
is  graphically  shown  in  Fig.  14,  where  the  long-term  and  short-term  power 
spectra  are  represented  along  with  the  power  spectra  developed  by  the  National 
Bureau  of  Standards  (NBS).  The  two  strai^t  lines  of  slope  -8/3  bracket  the 


Frg.  14.  Measurement  Spectra  Results 


variable  levels  observed  by  NBS,  and  are  considered  to  represent  the  atmospheric 

fluctuations  (scaled  to  6.  5-mile  path).  The  LIR  long-term  spectra  compare  well 

-4  -3 

with  this  spectra  in  the  10  to  10  cycles/ second  region  before  breaking 

away  into  a  plateau  region  whose  level  is  highly  variable  (as  indicated  by  two 

*See  ESD-TDR-64-103,  Vol.  I,  Figs.  1-15  and  1-16. 
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typical  short-term  spectra).  See  Fig.  14.  Thus,  the  difference  between  the 
plateau  and  the  (-8/3)  regeion  represents  "free"  or  unwanted  signal  over  and 
above  what  may  be  considered  true  atmospheric  data. 

Since  power  spectra  do  not  show  the  direct  current  component  of  the 
atmosphere,  which  is  of  primary  interest  to  us,  it  is  desirable  to  compare  the 
variable  component  to  the  average  level.  This  data  was  taken  under  conditions 

3 

of  approximately  10  gm/m  of  water  vapor,  or  60  N -units  of  refractivity  in  the 
path.  A  1  percent  measurement  of  this  quantity  would  yield  0.  6  N-unit  accuracy. 
For  the  sake  of  simplicity,  the  plateau  region  may  be  approximated  as  a  constant 
power  per  unit  bandwidth,  up  to  a  complete  cutoff  frequency  representative  of  the 
break  point.  A  locus  of  such  plateau  values  and  break  point  frequencies  is  given 
by  the  line  of  slope  -2  for  a  0.  6  N  unit  fluctuation. 

It  can  be  seen  that  the  plateau  region  can  exceed  the  0.  6  N -fluctuation  at 

-3  -1 

frequencies  as  low  as  10  cycles  per  second  and  as  hi^  as  10  cycles  per 

second.  (One  spectra  had  a  plateau  lower  than  the  lowest  shown  in  Fig.  14). 

The  LIR  spectra  have  been  scaled  upward,  from  phase  measurements  of 
dispersivity,  to  represent  the  path  refractivity  by  the  ratio  of  refractivity  to 
dispersivity,  1580.  Without  this  scaling,  the  data  occur  in  the  vicinity  of  the 
Norton*  and  Barrows  explanation  of  the  spectra.  **  The  length  of  the  arrow 
(Fig.  14)  is  representative  of  the  scaling  factor.  While  the  agreement  is  good, 
the  effect  of  elevation  angle  on  the  fluctuations  and  the  observation  of  extremely 
large  amplitude  fluctuations  in  the  presence  of  visible  striations  indicate  that 
mechanisms  in  addition  to  "wave  optics"  must  be  considered. 


*See  ESD-TDR-64-103,  Vol.  I,  pp  1-143  to  1-153. 

**Ibid.  ,  pp.  1-155  to  1-193. 
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STATUS  OF  FEASIBILITY  TESTING 


The  status  of  feasibility  testing  of  the  LIR  may  be  stated  as  follows: 

(1)  smoothing  of  the  data  is  necessary  at  low-elevation  angles  (5  degrees); 

(2)  range  from  10  to  1000  seconds  for  (1  percent  precision),  100-second 
smoothing  is  typical; 

(3)  attainable  accuracy  has  not  been  demonstrated  experimentally 
because  of  the  inherent  difficulty  of  measurement  over  10-mile 
atmospheric  paths; 

(4)  instrument  precision  of  0. 13-degree  root-mean- square  corresponds 

3 

to  a  1  percent  measurement  of  nominal  humidity  (10  gm/m  )  over 
a  6-mlle  path;  and 

(5)  a  nominal  humidity  measurement  of  1  percent  corresponds  to  a 

0.  6  N-unit  measurement  of  a  0.  2  percent  refractivity  measurement. 

Comparison  with  other  Techniques 

A  comparison  with  other  methods  of  line-of-sight  radio  path  correction 
can  be  made  as  follows: 

(1)  end-point  measurements  require  long-term  averaging  time  of  approxi¬ 
mately  3  hours,  and  highly  averaged  spatial  models.  This  method 
yields  residuals  of  1-3  percent. 

(2)  Meteorological  profile  and  ray-tracing  methods  may  yield  3x1.6 
percent  correction  accuracy  with  an  estimated  time -constant  of 
about  1  hour. 

(3)  The  LIR  determines  an  average  over  a  path,  requires  10- to  1000- 
seconds  smoothing,  with  a  potential  accuracy  of  0.  6  N  -unit  or  0.  2- 
percent  accuracy. 
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OPTICAL  DEPOLARIZATION  AND  SCINTILLATION  MEASUREMENTS 

OVER  A  TERRESTRIAL  PATH 


Jurgen  R.  Meyer-Arendt* 

The  experiments  described  concern  scintillation  and  polarization  phenomena 
observed  over  a  25-km,  long,  single  terrestrial  path  and  over  a  double-path  of 
twice  that  length.  Both  types  of  phenomena  were  studied  using  a  scintillation- 
polarization  monitor  that  can  be  operated  in  two  ways  (Fig,  1),  These  are 


Fig.  1.  Single-Path  Scintillation — 
Polarization  Monitor 


♦National  Bureau  of  Standards,  Boulder,  Colorado 


referred  to  as  Single  Path  arrangement  and  Return  Path  arrangement.  In  the 
single-path,  light  source  (B)  and  receiver  (C)  are  located  at  opposite  ends  of  the 
test  path.  In  the  return  path  arrangement,  light  source  (A)  and  receiver  (C)  are 
both  placed  at  the  same  terminal.  The  light  is  returned  by  a  mirror  or  other  type 
of  reflector  (B).  We  have  used  both  a  plane,  first-surface  mirror  (practical  only 
for  short-path  lengths),  a  comer  cube  reflector  of  35  cm  cross-sectional 
diameter,  and  a  special  array  of  comer  reflectors  which  form  part  of  a  geodimeter. 

A  schematic  diagram,  and  a  photograph  of  the  instrument  set  up  for  return 
path  experiments,  are  shown  in  Figs.  2  and  3.  The  telescope,  denoted  (1),  has 
a  7.  5-cm.  aperture  achromatic  objective  (d)  of  114-cm  focal  length.  A  prismatic 
member  (2),  located  near  the  focal  plane,  has  a  divergent  lens  cemented  to  its 
front  face  and  a  plane  mirror  cemented  to  its  rear  face.  A  hole,  about  6  mm. 
in  diameter,  is  drilled  through  the  prism  and  the  mirror.  The  light  then  passes 
throu^  a  polarizing  filter  (3)  with  a  central  hole,  through  an  iris  diaphragm  (e), 
a  rotating  polarizing  filter  or  a  Glan-Thompson  prism  (5),  and  a  ground  glass 
(6),  and  is  incident  on  the  light-sensitive  surface  of  a  photomultiplier  tube  (7). 

The  light  comes  from  a  high-pressure  mercury  arc  (a)  and  passes  through  a 
collimating  lens  (b)  on  to  the  target  (c).  Part  of  the  light  can  reach  the  photo¬ 
tube  directly  via  the  window  (4).  A  large-diameter  polarizing  filter  is  placed 
in  front  of  the  collimating  lens  (b);  the  other  part  of  the  light  falling  directly, 
without  traversing  the  path,  on  the  photo  surface  is  polarized  by  the  outer  portion 
of  the  filter  (3). 

The  experiments  were  carried  out  from  9  to  19  June  1964  on  Maui,  Hawaii, 

using  the  same  path  that  had  been  used  for  radio  refraction  measurements  by 

[l]  * 

Norton  and  his  coworkers  ,  and  which  is  currently  being  used  again  by 


♦Numbers  in  brackets  designate  References  on  page  122. 
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Polarization  Monitor,  Schematic 


Fig.  3.  Polarization  Monitor,  Control  Panel 

Two  Representations  of  Double-Path  Scintillation 
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M.  C.  Thompson,  Jr.  ,  and  his  group.  The  island  of  Maui  offers  both  a  semi- 
tropical  humid  climate,  (easily  accessible  by  car,  on  the  top  of  Mt.  Haleakala 
(3055  m. ))  an  arid  climate  where  subfreezing  temperatures  sometime  prevail. 
Between  these  two  regions  there  are  heavy  air  turbulence  and  hazy-to-dense 
clouds  moving  into  and  out  of  the  path,  smoke  from  burning  sugar  cane  fields, 
and  almost-constant  strong  winds  (the  trade  wind,  blowing  mostly  in  a  northerly 
direction),  normal  to  path  of  observation.  The  length  of  the  path  (one  way) 
is  24.  88  km. 

The  results  of  our  experiments  can  be  briefly  described  as  follows: 

(a)  The  single  path  recordings  show  strong  scintillation,  i.  e.  ,  intensity 
fluctuations,  with  occasional  high  spikes  or  flashes  of  light  of  short 
duration. 

(b)  Single-path  polarization  records  show  a  well  discernible  sinusoidal 
curve  on  which  scintillation  fluctuations  are  superimposed. 

(c)  In  all  return  path  experiments,  the  noise  level  is  higher  because 
much  light  is  back  scattered  in  the  path.  The  reflectors  act  as  a 
(stationary)  secondary  light  source,  producing  essentially  the  same, 
and  sometimes  even  less  pronounced,  scintillation,  an  effect  that 
might  be  the  result  of  the  larger  cross  section  of  the  reflectors  as 
compared  to  the  primary  light  source. 

(d)  In  the  return  path  polarization  experiments,  the  influence  of  back 
scatter  is  rather  distracting.  It  can  easily  be  shown  that  the  sinus¬ 
oidal  curve  ordinarily  obtained  is,  in  no  small  part,  produced  by  the 
polarized  light  returned  by  scatter  from  the  initial  light  beam.  This 
is  evident  from  the  fact  that  moving  clouds  which  completely  block 
out  the  return  signal  only  moderately  reduce  the  sinusoidal  signal 
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INTENSITY 


recorded  (see  the  last  one-third  in  the  upper  tracing  in  Fig.  4),  In 
order  to  minimize  this  effect,  an  iris  diaphragm  was  placed  in  the 
prime  focus  of  the  telescope  objective  and  carefully  adjusted  so  as 
to  admit  light  forming  the  reflector  image  only.  The  result  is  shown 
in  the  lower  portion  of  Fig.  4. 


TIME 

Fig.  4.  Return  Path  Polarization 

In  summary,  we  come  to  the  following  conclusions: 

(a)  Based  on  the  experiments  performed,  there  is  no  evidence  of 

depolarization  produced  by  atmospheric  turbulence.  Attenuation, 

caused  by  haze,  limits  the  signal  transmission  before  any  depolariza- 

[2] 

tion  becomes  evident.  Polarization  phenomena  caused  by  aerosols 

[3] 

or  by  background  skylighr  are  excluded  in  this  study. 
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(b)  Scintillation  is  a  major  limitation  to  optical  telecommunication.  Its 
Influence  increases  with  the  distance  traversed,  the  characteristics 
of  the  atmosphere  traversed  (turbulence,  etc. ),  and  with  the  degree 

of  colllmation  of  the  light  source  (higher  collimation  is  less  favorable). 

(c)  Any  receiver  for  optical  information  transmittal  based  on  polarization 
must  have  provision  for  compensation  of  scintillation  fluctuations. 

This  can  be  accomplished,  for  instance,  by  a  beam-splitter,  two¬ 
photocell  arrangement.  (An  interesting  proposal  for  a  photoelectric 
method  for  simultaneously  measuring  all  four  Stokes  parameters  is 
found  in  Reference  4. ) 


Finally,  it  should  be  mentioned  that  the  results  presented  in  this  paper 
are  preliminary  in  nature,  and  that  they  will  be  further  analyzed.  It  appears 
probable  that  some  kind  of  turbulence  might  arise  from  thermal  expansion  pro¬ 
duced  by  a  very  strong  laser  beam.  This  possibility  is  currently  being  explored 
at  the  National  Bureau  of  Standards,  Boulder  Laboratories,  in  a  series  of 
experiments  in  which  a  laser  pulse  is  directed  into  one  arm  of  a  Mach-Zehnder 
interferometer  (Fig.  5).  Local  perturbations  would  then  show  as  distortions 
of  the  interference  pattern  produced  in  the  interferometer  by  light  from  another, 
auxiliary  light  source. 
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Mach-Zehnder  Interferometer 
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FADING  ON  MICROWAVE  LINE-OF-SIGHT  PATHS 


H.  T.  Dougherty* 


INTRODUCTION 

The  received  signals  for  line-of-sight  propagation  paths  tend  to  vary  with 
time  because  of  the  refractive  index  gradient  variations  along  the  propagation 
ray  path.  Although  this  is  experienced  at  lower  frequencies,  the  effect  can 
become  particularly  severe  at  microwave  frequencies;  fading  can  result,  which 
reduces  the  reliability  of  the  system  below  an  acceptable  value.  This  fading  can 
be  considered  under  two  general  categories. 

1.  Multipath  or  phase  interference  fading. 

This  effect,  the  result  of  refractive  index  gradient  variations, 
is  generated  by  the  interference  between  the  direct  wave  and  (a)  the 
specular  component  of  the  ground- reflected  wave;  (b)  nonspecular 
components  of  the  ground-reflected  wave;  (c)  partial  reflections 
from  elevated  layers  or  atmospheric  sheets;  or  (d)  additional  direct 
wave  paths  produced  by  horizontally  distributed  changes  in  refractive 
index.  These  are  illustrated  in  Fig.  1. 

2.  Fading  caused  by  attenuation  or  defocuslng 

This  results  from  partial  isolation  of  the  transmitting  and 
receiving  antennas  because  of  (a)  intrusion  of  the  earth's  surface 
into  the  propagation  path  (diffraction  fading);  partially  reflecting 


*National  Bureau  of  Standards,  Boulder,  Colorado 
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elevated  layers  interpositioned  between  the  terminal  antenna  eleva¬ 
tions;  or  (c)  ducting  formations  containing  only  one  of  the  terminal 
antennas.  These  are  illustrated  in  Fig.  2.  A  comprehensive  descrip¬ 
tion  of  the  foregoing  and  additional  fading  mechanisms  is  given  by 

[l]* 

Bechmann  and  Spizzichino 
BEADING  CHARACTERISTICS 

When  fading  is  observed  on  microwave  line-of-sight  paths,  it  is  useful 
to  characterize  the  fading.  This  may  be  done  in  terms  of  the  duration  and 
depth  of  the  fades,  their  diurnal  variation,  and  any  selectivity  with  respect  to 


♦Numbers  in  brackets  designate  References  on  page  135. 
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space  or  frequency.  The  characterization  may  then  permit  identification  of  the 
fading  mechanism  (s)  and  suitable  techniques  to  reduce  their  effects. 

The  depths  of  the  fades  encoimtered  in  multipath  phenomena  can  be  quite 
severe,  depending,  of  course,  upon  the  effective  reflection  coefficients  or  the 
relative  amplitudes  of  the  component  waves.  The  interference  between  the 
direct  wave  and  the  reflected  wave  from  groimd  or  elevated  layer  can  produce 
fades  persisting  from  minutes  to  hours.  During  such  fades,  the  nonspecular 
component,  normally  small  for  grazing  angles,  can  cause  additional  interference 
resulting  in  deeper,  more  rapid  fades  whose  durations  are  of  the  order  of  sec¬ 
onds.  The  multipath  fades  tend  to  be  less  frequent  during  the  afternoon  hours 
when  conditions  more  closely  approach  the  "standard  atmosphere"  for  which  most 
paths  are  designed. 
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These  fades  are  selective  in  space  and  frequency  so  that  their  effects  can 

[2,  3,  4,  5] 


be  reduced  substantially  by  space  or  frequency  diversity. 


Other  tech¬ 


niques  are  also  available.  For  example,  the  Japanese  have  employed  anti- 

r  6 1 

reflective  wave  antennas^  ^  to  reduce  the  contribution  of  ground-reflected 
waves.  Furthermore,  for  propagation  over  rough  terrain,  the  proper  choice 
of  antenna  sites  can  also  reduce  or  eliminate  the  ground-reflected  wave. 


Attenuation  fading  caused  by  diffraction  of  the  earth  can  persist  for  several 
hours  and  to  depths  of  up  to  20  or  30  db.  The  depth  of  fading  is  space  and  fre¬ 
quency  selective,  and  tends  to  be  less  frequent  during  the  afternoon  hours.  This 
may  be  avoided  by  antenna  heights  which  provide  clearances  of  several  Fresnel 
zones  for  the  standard  atmosphere.  Deep  fades,  up  to  16  db,  lasting  for  hours 
or  days  (fadeouts),  characterize  the  attenuation  produced  by  ducts  and  elevated 
layers.  These  are  selective  in  space  but  not  in  frequency,  and  tend  to  be  less 
prevalent  during  the  daylight  hours.  These  fadeouts  can,  therefore,  be  en¬ 
countered,  to  a  degree,  by  space  diversity. 


MISTRAM  I  OBSERVATIONS 


The  multipath  phenomenon  has  been  observed  as  a  source  of  fading  on 
the  MISTRAM  I  System  of  the  Atlantic  Missile  Range  (AMR).  This  has  been 
particularly  troublesome  for  that  part  of  the  communication  link  subsystem  which 
consists  of  microwave  relay  circuits  transmitting  all  signals,  including  data  and 
timing  synchronization  between  stations  in  the  System.  Signal  recordings  at 
approximately  7  gc,  obtained  from  17  September  1963  through  20  April  1964, 
were  provided  by  the  Pan  American  Technical  Staff  of  the  Air  Force  Eastern 
Test  Range  (AFETR)  for  two  of  the  overland  links. 
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The  following  is  the  result  of  the  analysis  by  the  National  Bureau  of 
Standards  (NBS).  Both  paths  were  similar,  30.  5  km,  in  length  with  effective 
antenna  heights  of  41, 1  and  47.  7  meters.  A  portion  of  the  recorded  signal 
exhibiting  deep  fades  is  presented  in  Fig.  3.  The  numerous  fades  of  the  first 
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Fig.  3  Deep  Fades  on  MISTRAM  System 


hours  or  so  shown  are  typical  of  about  1.  4  percent  of  the  observation  period. 

The  signal  during  the  next  hour  shown  is  representative  of  about  5,  6  percent 
of  the  period,  and  the  last  hour  is  ts^jical  of  the  other  93  percent  of  the  time. 

Even  during  the  hours  when  fading  was  experienced  the  signal  generally  returned 
to,  or  above,  the  free-space  value  (a  basic  transmission  loss  of  approximately 
139  db)  —  the  behavior  that  is  characteristic  of  a  direct  plus  ground-reflected 
wave  type  of  phase  interference.  The  occurrence  of  the  fading  was  not  predictable. 
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since  it  exhibited  no  definite  diurnal  pattern;  that  is,  there  were  many  days  for 
which  no  fading  occurred.  When  fading  did  occur,  it  occasionally  extended  into 
daytime  hours. 

These  particular  links  employed  frequency  diversity  reception  for  which 
Af  was  250  me.  A  comparison  of  the  recorded  signals  for  the  two  frequencies 
is  shown  for  a  period  of  several  hours  in  Fig,  4,  The  strong  similarity  in  the 
signal  variation  pattern  is  as  expected  for  multipath  fading.  The  "lag"  between 
the  two  patterns  is  generally  not  constant  because  it  is  a  function  not  only  of  the 
frequency  spacing  and  the  gradient  of  refractive  index  involved,  but  also  of  the 
time  rate  of  change  of  the  gradients  encountered.  For  certain  strong  negative 
gradients,  the  lag  will  go  to  zero,  making  the  frequency  diversity  ineffective. 

This  effect  was  illustrated  with  a  film*  following  the  oral  presentation; 

[5] 

the  film  showed  the  results  of  an  experimental  study'"  ^  of  a  48.  3-km  over-water, 
microwave  line-of-sight  propagation  path.  By  employing  a  frequency  sweeping 
technique  for  both  the  transmitter  and  the  receiver  local  oscillator,  Lewin 
obtained,  for  two  vertically  spaced  receiving  antennas,  an  oscilloscopic  compari¬ 
son  of  the  two  multipath  lobing  patterns.  The  transmitter  frequency,  4.  3  gc  was 
swept  ±  300  mcs.  The  local  oscillator  (IF  =  60  me,  BW  =  10  me,)  was  swept 
at  a  much  slower  rate.  The  height  of  the  transmitting  antenna  was  198  meters 
and  the  heights  of  the  receiving  antennas  were  169.  8  and  164  meters.  The  film, 
representing  six  days  of  observations  sampled  at  one-minute  intervals,  showed 
the  relative  displacement  of  the  two  patterns  as  the  atmospheric  conditions 
changed. 


*Film  based  on  L.  Lewin' s  experimental  results. 
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Fig.  4  Comparison  of  Signals  at  Two  Frequencies 
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Figure  5  indicates  a  few  of  the  situations  shown  in  the  film.  The  plots  of 
and  versus  f  represent  the  envelopes  of  the  oscilloscope  traces  for 
the  upper  and  lower  antennas,  respectively.  The  oscilloscope  sweep  was 
synchronized  with  the  sweep  frequency  of  the  transmitter.  The  plot  at  the  top 
of  the  slide  is  representative  of  the  condition  which  was  generally  observed  for 
optimum  spacing  between  antennas.  The  plot  at  the  center  of  the  figure  repre¬ 
sents  the  observed  comparison  for  an  excessive  antenna  spacing.  The  plot  at 
the  bottom  shows  the  effect  of  attenuation  of  the  field  at  the  lower  antenna 

possibly  due  to  partial  isolation  of  the  lower  antenna.  Sample  illustrations  have 

[5] 

been  given  by  Lewin. 


Fig.  5  Envelopes  of  Signals  from  Spaced  Antennas 
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DIVERSITY  RECEPTION 


The  signal  recordings  mentioned  above,  and  those  illustrated  in  Fig.  4, 
indicate  that  a  substantial  improvement  would  be  expected  for  a  diversity  recep¬ 
tion  system  which  provides  for  wide  variations  in  the  gradients  of  refractive 
[5] 

index.  Lewin  ^  describes  a  procedure  for  determining  the  optimum  antenna  or 
frequency  spacing  for  such  a  situation.  Briefly,  the  procedure  consists  of 
determining: 

(a)  the  terminal  antenna  heights  required  to  provide,  for  a  given  path 
length  and  operating  frequency,  a  specified  protection  (a  20-db  fade 
margin,  for  example)  against  diffraction  fading  for  those  extreme 
positive  values  of  refractive  index  gradients  that  might  be  encountered, 
and 

(b)  the  diversity  spacing  required  to  provide  the  specified  protection 
against  multipath  fading  for  those  extreme  negative  refractive  index 
gradients  that  might  be  encountered. 

The  MISTRAM  I  paths,  described  previously,  can  be  used  to  illustrate 
the  results  of  applying  such  a  procedure.  The  present  terminal  heights  are 
adequate  to  provide  a  20-db  fade  margin  for  positive  gradients  up  to  approxi¬ 
mately  +200  N  units/km. 

An  examination  of  the  wiresonde  data  available  at  various  times  during 
the  period  of  the  MISTRAM  I  signal  recordings  indicates  refractive  index  gradi¬ 
ents  generally  of  the  order  of  -75  N  units/km,  which  corresponds  to  an  effective 
earth's  radius  factor,  k,  of  1.9,  The  average  gradient  for  the  first  100  meters 
or  so  above  the  surface  ranged  from  -206  to  +75  N  units/km.  This  corre¬ 
sponds  to  a  range  of  k-values  from  approximately  -3.  1  to  0.  68.  To  insure  a 
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20-db  fade  margin  for  gradients  of  -314  N  units/km  or  greater,  a  diversity 
spacing  of  approximately  11  meters,  or  1150  me,  would  be  required. 

Either  the  antenna  or  the  frequency  spacing  would  provide  the  required 
20  db  fade  margin  for  diffraction  fading  or  multipath  (direct  plus  ground-reflected 
wave)  fading  for  these  particular  paths.  In  general,  however,  the  possibility 
of  attenuation  fading  because  of  elevated  layers  or  ducting  conditions  would 
suggest  that  the  antenna  spacing  provides  a  preferable  form  of  diversity  recep¬ 
tion.  In  addition,  the  use  of  spaced  antennas  permits  addition  of  the  two  received 

[5] 

signals  at  RF  by  means  of  a  phase  combiner.  This  technique  has  been 

[7l  [4] 

employed  successfully  both  in  England  and  Japan. 

CONCLUSIONS 

The  experimental  data  shows  that  frequency  diversity  spacing  of  250  me 
is  inadequate  for  the  paths  described.  The  author  expects  that  diversity  recep¬ 
tion  by  means  of  vertically  spaced  antennas  will  provide  a  substantial  improve¬ 
ment  in  system  performance  for  the  fading  mechanisms  commonly  encountered 
over  the  microwave  line-of-sight  paths  for  the  flat,  maritime  terrain  of 
MISTRAM  I.  Although  the  effects  of  elevated  layers  and  ducting  conditions,  if 
encountered,  may  be  partially  offset  by  the  diversity  reception  with  spaced 
antennas,  their  most  severe  effects  would  require  increased  antenna  gains, 
more  sensitive  receiving  systems,  or  increased  transmitter  power. 
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SOME  COMMENTS  ON  MICROWAVE  LINK  CHARACTERISTICS  AND  PROBLEMS 


B.  M.  Hadfield* 


GENERAL 

All  methods  of  radio  communication  between  two  points  on  the  earth's 
surface  suffer  from  defects,  but  the  variety  and  degree  diminish  as  the  trans¬ 
mitted  frequency  is  raised.  This  is  true  primarily  because  reasonably  sized 
directive  antennas  become  feasible,  so  that  the  transmitter  energy  can  be 
concentrated  into  a  narrow  beam  directed  at  the  receiving  point.  Hence, 
microwave  links  operate  roughly  above  1000  me,  are  strictly  line-of-sight.  and 
require  reasonably  high  masts  to  elevate  the  antennas  above  the  earth's  curva¬ 
ture  and  any  intervening  hills. 

Microwave  links  owe  their  existence,  as  highly  competitive  alternates 
for  overhead  wires  and  underground  cables,  to  commercial  exploitation  on  the 
basis  of  equal  or  greater  reliability,  hi^er  channel  capacity,  and  lower  main¬ 
tenance  costs  with  comparable  initial  costs.  By  comparison  with  other  forms 
of  radio  link,  microwave  methods  have  been  refined  until  there  is  literally  no 
comparison  in  terms  of  reliability.  For  example,  99.  9  percent  reliability  is 
accepted  as  conventional  practice  in  this  field.  The  equipment  needs  to  attain 
this  performance  are  well  known,  and  are  obtainable  from  many  manufacturers 
on  a  production  basis. 

A  narrow  beam  line-of-sight  link  provides,  in  effect,  an  invisible  cable 
between  the  terminals,  whose  power  loss  between  isotropic  antennas  is  given  by 
22  +  20  d/L  db  (where  d  is  distance  and  L  is  wavelength),  or  6  db/octave 

*The  MITRE  Corporation,  Bedford,  Massachusetts 


137 


of  distance  or  frequency.  External  noise  has  a  relatively  minor  effect  on 
design,  because  of  the  narrow  beam  and  the  high  frequency.  Hence,  the  only 
basic  design  factors  are 

(1)  the  generation  of  sufficient  power  to  overcome  the  overall 
attenuation  and  provide  a  sufficient  signal-to-noise  (S/N) 
ratio  at  the  receiver,  and 

(2)  the  effects  generated  because  the  beam  is  not  ideally  narrow 
and  has  to  traverse  a  medium  which  causes  it  to  deviate 
from  a  straight  line  between  the  terminals. 

While  the  following  sections  will  stress  the  problems  associated  with 
setting  up  and  maintaining  a  microwave  link,  it  must  be  emphasized  that  many 
of  these  problems  exist  in  other  methods  of  point-to-point  communication  to 
an  even  greater  degree,  together  with  other  difficulties  more  or  less  incapable 
of  being  overcome  except  by  a  usually  expensive  tradeoff.  Microwave  links  are 
not  expensive,  and  some  additional  figures  will  be  given  later  as  to  their 
extreme  reliability,  when  properly  installed. 

LIMITING  FACTORS  DEPENDENT  ON  LINK  DISTANCE 

Free -Space  Loss  (Between  Isotropic  Antennas) 

The  6  db/octave  of  distance  effect  is  small  compared  to  the  "basic"  loss 
at,  say,  1  mile.  This  "basic"  loss,  often  termed  the  launching  loss,  is  100  to 
120  db  at  1  nautical  mile  (n.  m. )  and  for  frequencies  between  1000  and  10,  000 
me.  Between  5  and  60  miles  and  1000  to  10,  000  me,  the  loss  varies  between 
114  and  154  db.  This  is  not  a  highly  significant  change  if  about  80  db  of  the 
"basic"  114  db  loss  can  be  neutralized  by  antenna  gains,  since  the  receiver 
sensitivity  can  be  less  than  0.  1  microwatt  (2700  uv  on  75  ohms).  At  10,000 
me,  antennas  of  4-foot  diameter  will  give  a  total  gain  of  80  db,  and  at  3000  me. 
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the  diameter  would  be  13  feet.  Distances  of  up  to  50  miles  have  been  achieved, 
on  a  99  percent  reliability  basis,  with  not  more  than  1  watt  of  generated  power. 

Mast  (Antenna)  Height 

The  link  distance  is  limited  primarily  by  how  high  a  mast  can  be  tolerated, 
or  engineered,  to  withstand  the  elements  with  minimum  twist  or  sway.  If  a 
maximum  of  150  feet  is  decided  on,  then  the  link  length  can  be  18  n.  m.  at 
1000  me,  23  n.  m.  at  3000  me,  and  26  n.  m.  at  10,  000  me,  assuming  that  the 
intervening  terrain  is  smoothly  curved. 

The  required  antenna  height  is  the  sum  of  a  component  proportional  to  the 
2 

square  of  distance  (d  /4.  5,  d  measured  in  nautical  miles),  and  a  component 
proportional  to  distance  and  inversely  proportional  to  transmitter  frequency. 

The  latter  component  becomes  negligibly  small  at  distances  of  about  50  miles. 

Antenna  height  for  line-of-sight  over  obstacles,  more  than  any  other 
factor,  limits  the  use  of  microwave  beamed  links  to  around  30  to  50  miles  for 
permanent  installations  where  some  choice  can  be  made  in  the  relative  disposi¬ 
tion  of  the  sites. 

Microwave  Link  Fading,  Signal  Strength  Margin  and  Design  Criteria 

Normal  propagation  losses,  at  6  db/octave  of  distance,  occur  mainly  in 
the  daytime,  when  the  lower  atmosphere  is  thoroughly  mixed  by  rising  currents 
and  winds.  Alternatively,  normal  propagation  losses  occur  when  the  humidity 
content  of  the  lower  atmosphere  is  low,  as  in  winter  months. 

Abnormal  propagation  losses  are  caused  either  by  abnormal  refraction 
index  gradients,  reflection  from  some  midpath  surface,  or  absorption  by  water 
vapor,  rain  or  snow.  Whereas  absorption  gives  a  steady  decrease  in  signal 
strength  for  the  duration  of  the  cause,  refraction  index  changes  and  reflections 
give  rise  to  variable  decreases  known  as  fading. 


139 


In  regards  to  absorption  losses ,  it  is  sufficient  to  say  that  they  are  propor¬ 
tional  to  distance,  to  the  density  of  the  water  content,  and  to  the  square  of 
frequency.  At  around  the  usual  microwave  link  distance  of  20  to  30  miles, 
absorption  can  be  a  serious  matter  for  frequencies  above  3000  me.  In  any  event, 
absorption  losses  should  always  be  allowed  for,  since  any  relatively  steady 
depression  of  signal  strength  only  tends  to  enhance  the  fading  losses. 

Signal  variations  caused  by  reflections  are  primarily  the  result  of  overall 
bending  or  twisting  of  the  beam  path  between  the  antennas,  because  of  refraction 
index  changes,  and  only  occasionally  the  result  of  a  dense  atmospheric  layer 
above  the  line-of-sight.  The  effects  of  beam-edge  reflection  from  the  ground  in 
the  midpath,  are  normally  avoided  by  stipulating  an  additional  height  to  the 
antennas  over  and  above  that  for  insuring  a  line-of-sight.  The  additional  height 
is  proportional  to  distance  and  inversely  proportional  to  frequency,  i.  e.  ,  about 
60  feet  at  20  miles  for  3000  me. 

For  the  main  cause  of  fading,  namely  refractive  index  gradient  changes, 
the  variable  loss  in  signal  strength  is  produced  by  the  partial  cancellation  of 
energy  between  beam  ray  paths  that  have  travelled  different  distances.  This  is 
the  well  known  multipath  effect.  Almost  all  fading  levels  of  more  than  15  db 
below  normal  daytime  signal  level  are  caused  by  this  interference  between 
direct  and  Indirect  ray  path  lengths  within  the  normal  beamwidth. 

Even  under  normal  conditions,  the  beam  is  bent  downward  by  the  normal 
decrease  in  atmosphere  dielectric  constant  with  height.  This  is  allowed  for  in 
design  by  assuming  that  the  earth's  radius  is  multiplied  by  4/3,  thus  giving 
basic  antenna  heights  somewhat  less  than  would  be  required  for  a  tangential 
line-of-sight  path. 

Changes  in  humidity,  pressure  and  temperature  cause  both  local  and 
widespread  changes  in  dielectric  constant,  and  even  inversions  of  the  normal 
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gradient.  Thus,  the  received  signal  may  come  successively  from  various  parts 
of  the  beam  lobe  structure,  or  from  two  or  more  different  ray  paths.  In  addition, 
the  overall  bending  of  the  beam  may  occasionally  cause  it  to  hit  a  highly  reflective 
midpath  surface,  such  as  water  or  marshy  land.  With  the  lower  portion  of  the 
beam  being  reflected  and  the  upper  still  going  direct,  conditions  are  similar  to 
VHF  space  wave  propagation  in  which  the  propagation  loss  is  12  db/octave  of 
distance. 

Thus ,  if  the  fading  is  attributed  to  multiple  ray  paths  generated  by  changes 
in  the  vertical  refractive  index  gradients ,  then  vertical  space  diversity  of  an¬ 
tennas  should  help  to  overcome  fading.  But  if  the  fading  is  the  result  of  a  mid¬ 
path  reflection  or  to  mere  bending  at  the  beam,  space  diversity  is  unlikely  to 
help.  Frequency  diversity  will  also  help  in  those  cases  produced  by  multiple  ray 
paths.  There  is  always  a  very  fine  structure  associated  with  fading,  and  where 
this  is  important,  frequency  diversity  is  the  only  help.  Fine  structure  fading  can 
produce  short  fades  of  more  than  45  db  below  the  normal  level. 

When  the  normal  reduction  of  atmosphere  dielectric  constant  with  height 
becomes  inverted,  then  the  beam  will  be  deflected  up  off  the  receiving  antenna, 
and  only  a  weak  signal  will  be  received  by  diffraction  from  the  intervening 
ground.  This  is  particularly  true  if  the  design  has  not  allowed  enough  initial 
beam  clearance  from  the  intervening  terrain,  with  a  late-evening  radiation- 
type  ground  fog,  and  in  late  summer. 

Alternatively,  it  is  possible  to  get  abnormally  high  signal  strengths, 
.especially  when  the  atmosphere  is  still  and  the  earth  is  hotter  than  the  air, 
because  of  the  formation  of  a  waveguide  atmosphere  duct  which  prevents  the 
beam  from  spreading  normally. 
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In  general,  there  is  always  a  diurnal  change  of  signal  level,  and  use  of 
too  narrow  a  beam  or  insufficient  path  clearance  very  much  enhances  multipath 
fading.  Multipath  fading  level  variations  of  between  +  15  and  -  30  db  over  30 
miles  are  liable  to  occur. 

A  specimen  design  allowance  for  fading  at  2000  me  would  be  small  or 
zero  for  short  distances  up  to  10  miles,  since  beam  deviations  clearly  increase 
with  distance.  From  10  to  20  miles,  the  design  allowance  increases  to  20  db. 
Above  25  miles,  however,  it  is  found  that  proportionally  larger  allowances 
are  required;  a  margin  to  cope  with  signals  30  to  35  db  below  normal  is  usually 
quite  adequate. 

The  fading  rate  tends  to  increase  with  frequency,  but  the  integrated 
duration  of  low  signal  levels  is  about  the  same  at  all  frequencies  between  1000 
and  10,  000  me. 

One  of  the  most  reliable  microwave  systems,  the  Bell  TD-2  between  the 
East  and  West  coasts,  uses  20  to  30  mile  links  (sited  on  mountains),  and  is 
designed  to  accommodate  fading  from  +5  to  -25  db,  with  average  rates  of  fade 
of  10  db/second. 

The  long-term  reliability  of  a  single  30-mile  link,  without  diversity,  is 
about  90  percent  for  a  10-db  fading  margin,  99  percent  for  a  20-db  margin, 
and  99.  9  percent  for  a  30-db  margin,  at  frequencies  of  about  2000  me. 

In  general,  abnormal  fading,  caused  by  abnormal  climatic  conditions  or 
to  a  difficult  route,  is  usually  compensated  for  by  diversity  transmission  and 
reception,  using  spaced  antennas,  frequencies  or  different  polarizations. 
Diversity  reception,  of  one  term  or  another,  is  found  to  be  imperative  for  over¬ 
water  paths,  if  high  reliability  is  desired. 
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Attenuation  Due  to  Atmospheric  Absorption  or  Scattering 


This  problem  is  peculiar  to  microwave  links,  since  it  is  significant  for 

frequencies  above  1000  me.  Its  magnitude  is  also  a  function  of  distance. 

Snow  or  heavy  rainfall  are  the  only  conditions  which  produce  absorption  at 

lOR 

frequencies  up  to  10,000  me,  and  a  reasonable  formula  is  ^.y^avelength)^ 

db/mile,  where  R  is  the  rainfall  in  inches/hour,  and  the  wavelength  is  in 
centimeters.  This  attenuation  is  only  significant  for  frequencies  above  3000 
me,  amounting  to  20  db  at  10,000  me,  for  20  n.  mi.  ,  and  1  inchAour  of 
rainfall.  This  is  of  the  same  order  as  multipath  fading  at  20  miles  and 
10,  000  me. 

Attenuation  caused  by  this  condition  is  a  serious  matter  because  it  can 
result  in  complete  outage  for  the  duration  of  the  storm.  It  is  primarily  caused, 
however,  by  the  use  of  frequencies  above  5000  me  and  by  the  lack  of  a  system 
design  margin.  Users  of  frequencies  above  5000  me  rely  on  the  normal  design 
fading  margin  to  compensate  for  this  effect,  since  it  is  true  that  normal  fading 
is  not  present  during  rainfall  because  of  the  turbulent  condition  of  the  atmosphere. 
Also,  the  probability  of  a  prolonged  storm  of  1  inch/hour  over  the  whole  of 
a  microwave  system  is  not  very  great,  according  to  present  operating  experience 
in  this  country.  However,  in  other  parts  of  the  world,  for  example,  the  monsoon 
regions,  outages  attributable  to  these  conditions  would  be  very  probable. 

In  general,  then,  although  rainstorm  and  snowstorm  attenuation  is 
possible  on  microwave  links,  in  practice  it  need  not  be  a  limiting  factor, 
provided  frequencies  below  5000  me  are  used  with  Unit  lengths  less  than 
30  miles,  and  provided  the  design  is  engineered  with  an  adequate  margin 
for  fading  all  descriptions. 
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OTHER  DESIGN  FACTORS  NOT  PRIMARILY  DEPENDENT  ON  DISTANCE 
Route  Survey  and  Terminal  Siting 

In  common  with  all  point-to-point  communication  systems,  but  more 
particularly  on  account  of  the  line-of-sight  requirement,  beamed  microwave 
links  have  to  be  readily  adaptable  to  the  terrain  between  the  terminals ,  and 
In  llic  available  terminal  sites  which  are  economically  situated  with  respect 
to  tho  user. 

Much  has  been  written  on  the  subject  of  how  to  survey  such  links,  most 
of  which  reduces  to  a  requirement  that  the  designer  should  have  a  complete 
picture  both  in  plan  and  cross  section  of  the  terrain  between  the  desired 
terminals,  during  both  winter  and  summer  (due  to  the  changing  foliage).  By 
normal  "hand  and  foot"  methods,  this  sort  of  survey  can  be  an  exhausting 
process,  requiring  perhaps  several  visits  if  the  early  assumptions  based  on 
survey  maps  prove  inaccurate,  as  is  often  the  case. 

However,  aerial  photographic  surveys  in  depth,  plus  aerial  contour 
mapping,  plus  trial  propagation  tests  (with  proper  regard  for  the  prevailing 
atmospheric  conditions)  between  helicopters  equipped  with  standard  trans¬ 
mitters  and  receivers  hovering  over  the  proposed  sites  at  the  proposed  mast 
heights,  all  reduce  the  time  element  for  the  geographic  system  design  to  a 
matter  of  days  rather  than  weeks. 

In  general,  then,  the  main  requirement  in  designing  a  microwave  system 
(such  as  operating  frequency,  power  generated,  receiver  sensitivity,  channel 
usage,  antenna  sizes  and  mast  heights)  is  to  have  acquired  sufficient  experi¬ 
ence  of  operating  microwave  links  in  the  designated  area.  There  is  plenty  of 
such  experience  for  links  anywhere  in  the  United  States,  and  a  fair  amount 
for  other  latitudes  and  climates. 
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Mast  Rigidity,  and  the  Use  of  Reflectx)rs 


Because  the  mast  heights  for  links  of  20  miles  or  more  are  greater  than 
loo  feet  for  a  smooth  curved  earth,  and  because  thay  have  to  support  antennas 
(or  reflectors)  weighing  some  400  to  1000  pounds,  with  beamwidths  of  4  degrees 
or  less,  mast  rigidity  is  an  important  factor  in  overall  long-period  system  per¬ 
formance.  A  twist  of  more  than  2  degrees  will  clearly  give  received-signal 
level  variations  comparable  to  normal  multipath  fading. 

Normal  atmospheric  loading  conditions  are  1/2-inch  of  ice  and  85 
m.  p.  h.  winds,  for  installations  in  the  United  States,  and  usually  the  mechani¬ 
cal  design  permits  2-1/2  degrees  angular  beam  deflection  for  antennas 
mounted  at  the  mast  head.  For  designs  using  a  ground-based  antenna  beamed 
upward  to  a  reflector  at  the  mast  head,  the  usual  permissible  mast  deflection 
is  1-1/2  degrees,  since  the  resulting  beam  deflection  is  twice  this  amount. 

A  typical  permanent  steel  mast  construction  uses  a  triangular  trussed 
beam  weighing  620  pounds/ 20-foot  section,  and  can  be  erected  to  a  hel^t  of 
340  feet.  An  aluminum  mast  is  similar  in  design,  weighs  50  pounds/lO-foot 
section,  and  can  be  errected  to  200  feet. 

Portable  or  mobile  masts  are  now  available,  which  are  erectable  on  a 
small  baseplate,  by  means  of  an  integral  winch.  A  typical  example  is  a 
100-foot  mast  telescoping  into  a  12-foot  package  which  can  be  erected  by  4 
men  on  a  3-foot  square  baseplate,  in  2  hours,  including  guying.  The  total 
weight  of  the  mast  plus  reflector  is  about  800  pounds.  Its  rigidity  is  such  as 
to  limit  beam  deflections  to  -0.  75  degree  in  elevation  and  1.  5  degrees  in 
azimuth,  under  the  normal  ice  and  wind  loading  conditions. 

It  is  sometimes  possible,  and  necessary,  to  use  passive  reflectors 
situated  along  the  link  and  sited  so  as  to  get  around  or  over  some  obstruction 
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which  would  otherwise  demand  prohibiting  mast  heights.  However,  such 
reflectors  are  really  permanent,  since  to  intercept  a  reasonable  portion  of  the 
beams,  they  have  to  be  large  (20  feet  x  24  feet,  say),  very  rigid  (not  more 
than  0.  14  degree  twist),  and  very  flat  (to  0.  5  inch).  They  also  reduce  the  beam- 
width  to  about  0.  5  degree,  and  although  they  thereby  contribute  gain,  the  net 
result  is  an  enhancement  of  multipath  fading. 

Operating  Reliability  and  Maintenance  Costs  of  Commercial  Equipment 

and  Systems 

The  American  Gas  and  Electric  Company  has  a  5 -link  system,  one  link 
of  which  uses  a  passive  reflector  on  a  hilltop,  which  has  been  in  operation  for 
2  years.  They  have  kept  records  of  outages  for  any  cause,  including  regular 
maintenance  and  some  initial  troubles  while  learning  the  maintenance  techniques. 
Their  measured  system  reliability  is  99.  89  percent,  and  they  expect  to  improve 
this  to  99.  95  percent. 

Southern  California  Edison  Company  operate  40,  000  miles  of  private  wire 
lines  and  7000  miles  of  microwave  circuits.  The  maintenance  costs  of  their 
private  wire  lines  is  $0.  68/mile/month,  while  the  microwave  system  varies 
between  $0.  32  and  $0.  51/mile/month.  They  confirm  the  high  reliability  of 
microwave  systems ,  and  say  that  the  best  possible  wire  line  reliability  has 
been  found  to  be  98.  2  percent. 

Cleveland  Electric  Illumination  Company  has  an  extensive  microwave 
system  for  transfer -trip  relaying,  which  is  considered  to  require  the  most 
reliable  type  of  communications  channel.  Their  preliminary  tests  over  5 
months,  including  summer  fading,  showed  a  reliability  of  99.  98  percent  over 
a  distance  of  19  miles  during  which  no  false  signals  at  all  were  received. 

From  the  above,  and  other  evidence  too  long  to  be  incorporated  here, 
there  can  be  no  doubt  that  properly  engineered  beamed  microwave  links  have 
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reliabilities  as  good  as,  and  probably  better  than,  private  wire  and  cable  links, 
with  much  lower  maintenance  costs  and  probably  lower  initial  or  amortization 
costs. 

OBSERVATIONS  ON  THE  MISTRAM  FADING  PROBLEM 

While  we  have  not  had  enough  time  to  look  at  all  the  recorded  signal  level 
charts  in  detail,  it  would  seem  that  there  are  at  least  two  types  of  fading;  one 
which  is  of  short  period  with  deep  excursions,  and  the  other  which  is  slower 
and  superposed  on  a  general  reduction  in  signal  level. 

Frequency  diversity  should  compensate  for  the  short-period  deep  fades, 
but  the  general  depression  of  received  level  is  indicative  of  midpath  trouble, 
and  is  more  likely  to  be  overcome  by  increasing  the  antenna  height  than  by 
vertical  diversity  with  the  present  masts. 

It  is  pertinent  to  inquire  whether  any  propagation  tests  were  run  over  a 
period  of  a  year  between  the  baseline  sites  and  against  various  antenna  heights. 
The  whole  subject  of  microwave  link  siting  is  very  involved  and  very  important 
to  successful  operation.  The  designer  should  have  a  complete  picture  in  plan 
and  cross  section  of  the  terrain,  during  summer  and  winter,  and  tests  using  a 
pair  of  helicopters  are  quite  common. 

No  mention  has  been  made  in  the  presentations,  or  the  discussions,  of 
a  design  fading  allowance  in  the  MISTRAM  baseline  links  (30  db  is  a  common 
allowance  for  fading).  Furthermore,  in  a  well-designed  microwave  link, 
even  when  all  such  allowances  have  been  made,  an  operating  margin  of  20  db 
is  usually  still  available.  This  surplus  margin,  more  than  anything  else,  is 
what  accounts  for  the  very  hi^  reliability  of  conventional,  commercial 
microwave  links.  It  would  be  interesting  to  have  available  a  typical  breakdown 
of  the  transmitter,  propagation  and  reception  losses  and  margins,  for  the 
MISTRAM  baseline  microwave  links. 
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FINDINGS  AND  RECOMMENDATIONS  OF  THE  NATIONAL  ACADEMY  OF 

SCIENCE  "PANEL  ON  TRACKING  DATA  ANALYSIS"  IN  THE  AREA  ’ 

OF  TROPOSPHERIC  PROPAGATION  AND  RE  FRACTION 

Dr.  A.  T.  Waterman,  Jr.  * 

The  Panel  dealt  with  a  variety  of  problems  at  the  Eastern  Test  Range,  (ETR), 
of  which  propagation  was  one.  My  role  was  only  in  connection  with  the  propagation 
aspects. 

Everyone  else  who  has  spoken  here  has  said  something  about  what  he  has 
actually  done.  The  remarks  that  I  have  to  make  refer  just  to  giving  advice  to 
people  on  how  to  do  things.  I  might  mention  also  (I  think  it's  entirely  appropriate 
to  do  so)  a  report  that  was  put  out  a  couple  of  years  ago  by  the  National  Academy 
of  Sciences  Ad  Hoc  Panel  on  Electromagnetic  Propagation  of  which  Dr.  Straiten 
was  a  member.  Most  of  the  technical  investigations  of  propagation  problems 
associated  with  missile  guidance  and  tracking  were  comprehensively  summarized 
in  that  report.  As  a  result,  the  actual  amount  of  work  of  that  nature  remaining 
to  be  covered  by  our  Panel  was  relatively  small. 

There  are  two  basic  questions,  I  think,  with  regard  to  propagation  and  its 
relationship  to  the  role  of  tracking  accuracy.  One  is  the  question  of  what  should 
be  done  to  improve  certain  errors;  the  other  is  the  question  of  the  ultimate  limit 
of  accuracy  set  by  the  atmosphere.  Now,  it  seems  to  me  hi^ly  important  to 
recognize  that  these  two  questions  are  of  a  different  nature.  The  first  question 
is  one  which  can  be  answered  by  engineering  techniques:  by  taking  existing 
knowledge  and  turning  it,  through  one  approach  or  another,  to  best  advantage  in 
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working  out  a  methodology  for  estimating  errors.  The  second  question,  the 
ultimate  limit  of  accuracy,  is  more  a  scientific  than  an  engineering  question  — 
the  sort  of  thing  one  should  be  quite  cautious  in  answering,  since  new  and  unfore¬ 
seen  developments  could  lead  to  a  substantial  improvement  in  the  future. 

I  think  it  is  also  worth  noting  (it  is  perfectly  obvious  to  everyone  here)  that 
there  have  been  a  couple  of  schools  of  thought  on  how  best  to  approach  this  basic 
propagation  problem.  One  places  emphasis  on  the  spectral  description  of  the 
atmosphere  and,  to  a  large  extent,  uses  statistical  techniques.  The  other  is 
more  deterministic.  It  thinks  of  the  atmosphere  in  terms  of  the  physical  processes 
taking  place.  It  is  faced  with  more  serious  measurement  problems,  but,  perhaps 
in  the  long  run,  might  hold  hope  for  greater  advances  because  of  a  more  profound 
eventual  understanding  of  the  processes  at  work.  The  viewpoint  taken  by  the 
Panel,  or  at  least  by  me  in  my  role  as  a  member,  has  been  that  these  two 
approaches  are  not  necessarily  in  conflict,  but,  rather,  that  they  can  be  used  in 
a  complementary  fashion. 

I'll  try  and  outline  briefly  the  recommendations  that  relate  most  nearly  to 
the  propagation  problems  involved.  These  fall  into  three  general  categories. 

One  is  meteorological;  the  other  is  basically  radio;  the  third  is  organizational. 

Meteorologically,  there  were  three  recommendations;  they  are  fairly 
obvious  sorts  of  things,  and  are,  in  fact,  being  implemented  to  a  certain  extent. 

The  first  has  to  do  with  vertical  soundings  in  the  Cape  Kennedy  region.  The 
recommendation  is  that  these  soundings  be  used  to  develop  insight  into  what  a 
sampling  error  is  in  a  sounding.  In  this  regard,  it  was  quite  gratifying  to  see  the 
records  which  D.  L.  Ringwald  showed.  I  think  there  is  also,  by  implication,  a 
need  for  a  multiplicity  of  soundings  and,  therefore,  of  any  improvement  in 
cheaper  methods  of  obtaining  accurate  soundings.  In  this  regard,  there  is  a 
possibility  of  further  development  of  the  concept  of  an  expendable  type  of 
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refractometer  as  distinct  from  one  in  which  the  entire  instrument  has  to  be  flown 
in  a  plane  and,  therefore,  is  expensive  enough  to  have  to  be  preserved. 

Second,  still  in  the  meteorological  vein,  there  is  a  recommendation  that  a 
surface  network  of  refractive  index  measurement  sites  be  set  up  in  the  region  of 
the  two  MISTRAM  stations,  with  spacings  ranging  from  the  order  of  antenna  sizes 
to  distances  comparable  with  the  largest  baseline  anticipated  in  guidance  systems, 
and  over  time  intervals  of  a  length  comparable  with  that  required  for  the  air  to 
drift  across  the  instrumentation  range.  (This  recommendation  relates  to  the 
approach  adopted  by  the  National  Bureau  of  Standards). 

The  third  recommendation  with  regard  to  meteorologic  measurements  is 
just  that;  with  a  large  amount  of  instrumentation,  the  opportunity  of  obtaining 
as  comprehensive  a  measure  of  atmospheric  characteristics  as  possible  should 
be  fully  exploited  on  some  occasions. 

Now,  with  regard  to  the  recommendations  relative  to  radio  measurements, 
there  were  several  things  that  were  considered.  Some  of  them  were  discarded 
for  reasons  that  were  not  wholly  related  to  propagational  or  meteorological 
problems ,  but  to  other  practical  considerations .  Some  of  the  things  considered 
were  (a)  the  establishment  of  fixed-path  measurements  in  the  range  area,  (b) 
a  continuation  of  measurements  obtained  in  the  past  in  which  location  by  optical 
means  of  a  beacon  in  a  plane  is  compared  with  the  operation  of  the  MISTRAM 
system  with  a  transponer  in  the  same  plane,  and  (c)  a  satellite  specially  out¬ 
fitted  with  a  MISTRAM  transponder ,  which  could  then  serve  as  a  means  of 
calibrating  the  system.  For  one  reason  or  another,  none  of  these  received  the 
full  support  of  the  Panel.  The  concept  of  fixed-path  measurements  has  limits 
in  the  range  area.  Most  fixed  paths  would  be  relatively  flat,  and  at  low  angles, 
and  in  situations  which  would  not  be  pertinent  as  far  as  missile  tracking  is  con¬ 
cerned.  There  is  also  the  problem  of  surface  reflections  which  would 
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complicate  the  situation  considerably.  The  concept  of  the  plane  flights  was 
argued  as  not  being  an  effective  means  of  calibrating  the  system,  apparently 
because  of  problems  pertaining  to  ballistic  cameras:  the  accuracy  and  cost 
of  reducing  data.  The  idea  of  using  a  satellite  equipped  with  a  MISTRAM  trans¬ 
ponder  was  argued  only  to  the  extent  that  special  satellite  should  be  used.  On 
the  other  hand,  the  concept  of  putting  a  transponder  on  a  satellite  that  did  have 
some  extra  room  on  it  received  favorable  response,  and  constituted  one  of  the 
recommendations.  The  most  important  recommendation  in  connection  with 
calibrating  the  propagation  errors,  however,  was  that  measurements  obtained 
during  a  missile  launch  should  be  continued  during  the  free-flight  portion  of 
the  missile  trajectory.  The  orientation  of  the  missile  would  be  stabilized, 
and  this  free-flight  portion  could  then  be  used  as  a  more  accurate  means  of 
calibrating  propagational  errors. 

Finally,  with  regard  to  the  organizational  recommendations  of  the  Panel, 
it  was  felt  that  the  amount  of  effort  and  manpower  that  was  being  placed  on  the 
problems  related  to  propagation  was  small  relative  to  its  overall  importance. 

The  propagational  problem  itself  will  be  one  that  continues  to  be  emphasized  and, 
therefore,  its  importance  might  be  reflected  in  an  organizational  set-up  in  which 
there  is  greater  centralization  of  control,  authority  and  responsibility  for 
refractive  and  propagational  measurements ,  as  well  as  for  the  analysis  of 
propagational  data. 

Now,  in  this  connection  there  are  some  advantages  that  have  not  been 
utilized,  and  the  discussion  in  the  recommendations  includes  mention  of  these. 
Work  done  in  other  places  and  the  methods  developed  by  other  groups,  such  as 
the  Cambridge  Research  Labs. ,  the  National  Bureau  of  Standards  and  The  MITRE 
Corporation,  could  be,  perhaps,  more  fully  utilized  in  the  reduction  of  actual 
missile  flight  data.  One  possibility  would  be  to  furnish  the  data  from  some  of 
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the  flights  directly  to  these  concerns.  They  could  perform  their  own  analysis 
of  the  data;  and  the  range  would  then  be  in  a  position  to  compare  the  relative 
merits  of  each  type  of  analysis. 

Finally,  I  think  it  is  important  in  this  field  to  keep  abreast  of  all  new 
developments.  Certainly,  equipment  such  as  the  line  integral  refractometer 
should  be  exploited  quite  fully  in  order  to  ascertain  its  ultimate  merits. 
Undoubtedly,  other  ideas  will  arise.  It  is  important  to  remain  constantly  alert 
to  such  new  suggestions. 
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SPACED  REFRACTOMETER  EXPERIMENT 


C.  F.  Martin* 

Before  discussing  what  the  spaced-refractometer  experiment  is  all  about 
and  what  we  hope  to  get  out  of  it,  I  would  like  to  briefly  review  the  refractometer 
background,  particularly  in  regard  to  calibration  problems.  Besides,  it  should 
also  be  of  some  interest  in  its  own  right  and  will,  I  hope,  perhaps  prompt  some 
helpful  suggestions  from  those  of  you  who  have  been  involved  with  similar  problems. 

REFRACTOMETER  CALIBRATION 

As  most  of  you  probably  know,  the  AFETR  now  has  in  operation  two  precision 
electronic  missile  trajectory  measurement  systems  which  are  commonly  referred 
to  as  MISTRAM  I  and  MISTRAM  H.  MISTRAM  I,  located  at  Valkaria,  Florida, 
about  35  miles  south  of  Cape  Kennedy,  consists  of  a  central  range -measuring 
station,  a  North-South  baseline  and  an  East-West  baseline,  each  with  range- 
difference  measuring  stations  at  10  K  feet  and  100  K  feet  from  the  central  station. 
MISTRAM  II,  located  on  the  island  of  Eleuthera,  about  300  miles  southeast  of 
Cape  Kennedy,  is  a  similar  system  except  that  it  lacks  the  short  baseline  stations 
and  the  baselines  intersect  at  about  120-degree  angle.  Specifications  for  these 
systems  call  for  a  range  measurement  accuracy  of  0.  4  foot  and  a  range  difference 
accuracy  of  5  ppm.  Although  such  specifications  can  be  met  by  the  system 
electronics,  the  current  state  of  knowledge  of  (tropospheric)  propagation  effects  is 
such  that  such  accuracies  cannot  now  be  met. 


*  Pan  American  World  Airway  s/GMRD,  Patrick  AFB,  Florida 
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In  designing  the  MISTRAM  systems,  it  was  recognized  that  propagation 
effects  would  likely  be  the  limiting  factor  in  final  system  accuracy,  so  provision 
was  made  for  including  as  an  integral  part  of  the  system  at  each  measurement 
station  the  finest  refractometer  that  money  could  buy.  Originally,  it  was  planned 
to  record  each  refractometer  reading  20  times  per  second  just  as  the  range  and 
range  differences  were  measured  20  times  per  second.  For  various  reasons, 
including  the  state  of  calibration  of  the  refrac  tome  ter  s  and  the  fact  that  the  data 
channel  was  needed  for  other  purposes,  this  rate  of  refractive  index  recording 
has  long  since  been  discontinued.  Nevertheless,  we  still  have  the  5  Valkaria 
refractometers  and  the  3  Eleuthera  refractometers  for  use  with  the  MISTRAM 
systems  in  whatever  manner  will  give  us  the  maximum  benefit  in  accounting  for 
atmospheric  effects  on  MISTRAM  data. 

I  am  not  quite  sure  of  the  reasons  for  the  selection  of  ground-level  refrac¬ 
tometers  by  General  Electric  as  the  basis  for  MISTRAM  refraction  corrections. 
No  doubt  they  were  influenced  to  a  considerable  extent  by  the  work  of  Bean  and 
Thayer  of  NBS,  even  though  they  reportedly  were  somewhat  less  than  pleased  by 
the  refraction  correction  accuracies  estimated  by  Bean  and  Thayer.  Perhaps  it 
was  also  felt  that  with  $100,  000  worth  of  microwave  refractometers  at  the 
MISTRAM  I  site,  highly  accurate  refraction  corrections  would  be  automatically 
insured.  If  this  were  the  case  they  were  certainly  overly  optimistic.  Almost 
equally  as  certain,  however,  is  the  fact  that  the  refractometers  can  be  put  to 
greater  use  than  they  have  been  until  now. 

Several  months  ago.  Dr.  A.  J.  Mallinckrodt,  who  is  a  consultant  on  radio 
propagation  problems,  to  Pan  American  Airways/GMRD,  proposed  an  experiment 
which  would  utilize  the  MISTRAM  refractometers  for  obtaining  information  on  the 
space-time  structure  of  the  refractive  index  near  groimd  level  in  the  MISTRAM 
area.  The  performance  of  this  experiment  would  require  that  at  least  6  refrac¬ 
tometers  be  in  operating  condition,  that  all  have  almost  identical  performance, 
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but  not  necessarily  that  the  zero  point  of  each  refractometer  be  determined.  In 
other  words,  the  experiment  would  require  only  relative  refractive  index  meas¬ 
urements.  But  since  considerable  effort  was  necessary  in  order  to  verify 
satisfactory  relative  operation,  and  since  absolute  operation  is  certainly  desirable, 
inasmuch  as  it  is  essential  for  use  in  any  refraction  correction  procedure  for 
MISTRAM  tracking  data,  it  was  decided  that  the  first  order  of  business  would  be 
to  put  all  eight  MISTRAM  refractometers  into  operational  condition  and  then  install 
them  on  site  after  the  conclusion  of  the  spaced-refractometer  experiment.  As  a 
matter  of  fact,  the  MISTRAM  I  refractometers  would  be  installed  on  site  for  the 
experiment  itself. 

The  type  of  refractometer  that  we  have  been  talking  about  was  designed  at 
the  National  Bureau  of  Standards  and  is  commonly  referred  to  as  the  Vetter  re¬ 
fractometer.  The  ones  that  we  have  were  built  by  Colorado  Research  Corporation 
which,  I  understand,  is  now  defunct;  but  presumably  this  has  no  connection  with 
our  calibration  difficulties.  The  design  of  the  Vetter  refractometer  system  is 
such  that  its  accuracy  of  operation  is  relatively  insensitive  to  the  electronics 
involved  and  we  need  concern  ourselves  only  with  the  two  resonant  cavities  which 
are  used. 

One  cavity,  referred  to  as  the  reference  cavity,  is  sealed  and  has  a  servo- 
controlled  probe  for  varying  the  resonance  frequency  of  the  cavity.  The  resonance 
frequency  of  this  cavity  must  be  calibrated  as  a  function  of  probe  position,  since 
the  only  readout  of  the  instrument  is  really  nothing  more  than  probe  position. 

The  other  cavity  is  called  the  sampling  cavity  and  is  the  one  which  is  ex¬ 
posed  to  the  atmosphere  whose  refractive  index  is  to  be  measured.  If  we  can 
establish  the  resonance  frequency  of  the  sampling  cavity  at  a  time  when  we  know 
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the  true  refractive  index  of  the  atmosphere  in  the  cavity,  then  we  can  obtain  the 
refractive  index  with  some  other  atmosphere  by  measuring  the  cavity's  resonance 
frequency  and  using  the  relation 
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The  resonance  frequency  of  the  sampling  cavity  we  measure  by  the  use  of 
the  reference  cavity;  the  servo-loop  positions  the  probe  of  the  reference  cavity 
such  that  the  two  cavities  have  the  same  resonance  frequency,  and  the  probe- 
position  calibration  then  gives  the  sampling-cavity  resonance  frequency.  The 
atmosphere  of  a  known  refractive  index  is  most  accurately  obtained  by  placing 
the  sampling  cavity  in  a  bell  jar  and  evacuating  it.  Since  the  refractive  index 
of  even  a  poor  vacuum  deviates  from  unity  by  an  insignificant  amount  for  our 
purposes,  we  have  an  absolute  calibration  of  the  instrument. 

This  all  sounds  rather  simple  and  straightforward,  but  it  really  doesn't 
turn  out  to  be  that  way.  That  this  might  be  so  we  might  suspect  by  considering 
what  we  are  actually  measuring  and  what  affects  our  measurements.  The  reso¬ 
nance  frequencies  are  dependent  only  upon  geometry  and  we  are  interested  in 
obtaining  an  overall  accuracy  of  better  than  1  N  unit  over  a  period  of  days. 

Over  short  periods  of  time  we  expect  to  detect  changes  of  0.  1  N  unit.  When  we 
realize  that  1  N  unit  corresponds  to  a  change  in  cavity  dimensions  of  1  ppm  and 
that  the  sampling  cavity  must  be  open  to  the  atmosphere  along  with  its  fluctuating 
temperature  and,  especially  in  Florida,  its  numerous  varieties  of  flying  and 
crawling  pests,  we  can  appreciate  the  possibility  of  difficulties  arising.  As  a 
matter  of  fact,  even  before  going  out  into  an  uncontrolled  environment  we  have 
encountered  some  problems. 

The  calibration  of  the  reference  cavities,  or  at  least  finding  a  workable 
routine  for  it,  can  be  considered  successfully  accomplished  after  considerable 
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effort.  Since  each  reference  cavity  is  a  little  bit  different  from  any  other,  and 
the  cavity  resonance  frequency  is  not  quite  a  linear  function  of  probe  position  in 
any  case,  each  one  must  be  calibrated  individually  and  its  own  calibration  curve 
of  probe  position  versus  resonance  frequency  drawn  up.  In  a  typical  case,  the 
error  incurred  by  taking  the  refractometer  readout  to  be  true  N  -unit  reading 
amounts  to  about  5  N  units.  Without  this  correction  we  might  as  well  not  use 
the  refractometer  at  all,  since  standard  meteorological  measurements  made  with 
a  little  care  give  better  results. 

The  calibration  curve  is  obtained  by  replacing  the  sampling  cavity  by  a 
sealed  calibration  cavity.  This  cavity  also  has  a  probe  which  is  positioned  by 
a  vernier  dial  on  the  end  of  the  cavity.  Its  calibration  of  resonance  frequency 
versus  vernier  reading  was  performed  by  the  Bureau  of  Standards  Calibration 
Services. 

Some  difficulty  has  been  encountered  with  a  rather  long-term  drift  in  one 
or  two  of  the  reference  cavities.  This  I  understand  is  most  likely  due  to  the  type 
of  seal  material  used  and  should  stabilize  after  a  warmup  of  several  hours. 

An  apparent  instability  of  the  sampling  cavities  which  has  been  encountered 
in  attempts  to  measure  their  temperature  coefficients  is  causing  more  concern 
at  the  moment.  I  use  the  word  apparent  because  I  don't  actually  believe  that  they 
are  as  real  as  indicated.  In  going  through  a  temperature  cycle,  a  hysteresis 
loop  is  expected  due  to  the  lag  in  tiie  temperature  compensating  device  of  the 
cavity,  but  the  loop  should  close  instead  of  being  3  to  4  N  units  away  from  the 
initial  reading  as  has  been  foxmd  to  be  the  case. 

In  order  to  insure  both  adequate  ventilation  of  the  sampling  cavity  and  to 
help  prevent  the  intrusion  of  undesired  pests,  a  forced  ventilation  shelter  (with 
screens)  is  being  designed  for  use  at  the  refractometer  operating  sites.  The 
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shelter  is  designed  so  that  it  may  be  readily  sealed  for  air  evacuation.  The 
evacuation  process  is  necessary  in  order  to  check  the  refractometer  zero  set, 
but  the  frequency  at  which  it  will  be  required  is  not  yet  known. 

SPACED  REFRACTOMETER  EXPERIMENT 

The  motivation  behind  the  spaced-refracto meter  experiment  lies  in  the 
current  paucity  of  information  on  refractive  index  fluctuations  over  long  periods 
of  time  at  widely  separated  spatial  points.  Considering  the  at  least  potential 
availability  of  accurate  refractometers  at  the  already  rather  widely  separated 
MISTRAM  stations,  it  is  proposed  to  use  these  instruments  —  together  with  one 
of  the  Eleuthera  refractometers  as  an  extension  of  one  of  the  baselines  —  to 
record  the  near  ground  level  index  of  refraction  at  a  rate  of  about  6/minute  for 
a  period  of  1  week.  The  current  plan  is  to  locate  the  extra  refractometer  at  the 
Azusa  Mark  n  site  on  Cape  Kennedy,  giving  a  maximum  spatial  separation 
between  refractometers  of  about  55  miles.  The  Azusa  location  was  selected,  in 
preference  to  earlier  suggestions  of  Tampa  or  Lakeland,  Florida,  partly  because 
of  its  accessibility  and  partly  because  1  found  it  very  difficult  to  believe  that 
atmospheric  conditions  at  either  of  these  towns  had  any  relation  to  conditions  high 
out  over  the  Atlantic  ocean. 

This  experiment  is  scheduled  to  be  carried  out  as  soon  as  the  refractometers 
are  calibrated.  At  the  present  time,  it's  impossible  to  say  when  this  will  be, 
but  a  crude  guess  would  be  sometime  in  September. 

We  might  now  briefly  consider  how  we  are  going  to  treat  the  data  and  how 
it  will  be  useful.  Primarily,  we  are  interested  in  obtaining  a  space-time  corre¬ 
lation  function  for  refractive  index  fluctuations.  As  a  matter  of  definition,  this 
function  is  given  by 


C  (p,  r)  =  <  S  N  (r,  t)  S  N  (r  +  p,  t  +  r)  > 
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where  the  <>  indicates  a  time  and  spatial  average.  How  we  handle  the  data  depends 
partly  upon  what  we  expect.  Measurements  made  at  a  single  point  have  supposedly 
indicated  that  the  time  correlation  function  is  approximately  exponential, 


K  the  atmospheric  inhomogeneities  really  were  frozen,  then  they  should  drift  by 
the  measurement  station  which  would  then  really  be  measuring  spatial  inhomogene¬ 
ities.  The  resulting  form  for  the  space-time  correlation  function  is  then 
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C(p,T)  =  e 

if  V  is  the  drift  velocity  of  the  inhomogeneities.  Even  on  physical  grounds, 
there  are  several  things  wrong  with  this  reasoning.  First,  one  should  not  expect 
perfect,  or  even  near  perfect,  correlation  of  space  and  time  fluctuations  except 
over  very  short  periods  of  time.  When  the  wind  direction  changes  many  tim^^ 
during  the  data  sampling  period,  one  could  hardly  expect  the  data  from  spaced 
refractometers  to  yield  a  well  defined  wind  velocity.  In  the  Cape  Kennedy  area, 
the  wind  direction  changes  from  a  predominantly  southeasterly  direction  during 
the  day  to  a  predominantly  northwesterly  direction  in  the  evening,  so  a  7 -day 
test  would  cover  at  least  14  major  changes  in  wind  direction. 

A  further  source  of  difficulty  is  that  even  a  simple  time  correlation  function 
will  not  necessarily  have  the  above  simple  form.  An  analysis  of  some  refractive- 
index  data  obtained  at  the  MISTRAM  sites  for  Mr.  Ringwalt  shows  that  the  function 
goes  to  zero  at  about  18  hours,  and  then  to  a  negative  value  of  -^0.  35  at  around  28 
hours.  After  this  period,  it  again  approaches  zero,  although  rather  slowly. 
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Such  behavior  is  somewhat  disconcerting  and  not  at  all  what  we  would  like. 
The  reason  we  suspect  to  be  a  rather  large  daily  fluctuation.  Some  of  Mr. 
Ringwalt's  data  is  now  being  analyzed  after  first  taking  out  the  daily  mean.  Hope¬ 
fully  this  procedure  will  yield  correlation  functions  which  die  off  more  rapidly 
and  have  better  behavior.  If  so,  we  will  then  feel  that,  for  the  one  week  data 
period,  we  can  adequately  represent  the  refractive  index  fluctuations  of  the 
atmosphere  as  the  sum  of  a  sinusoid  of  a  1-day  period  plus  a  random  component. 

The  use  to  which  a  space-time  correlation  function  can  be  put  will,  we 
hope,  lie  in  the  area  of  improving  the  accuracy  of  the  refraction  corrections  used 
in  MISTRAM  data  reduction.  In  addition,  an  assessment  can  be  made  of  the 
residual  errors  in  the  refraction  correction.  This  possibility  arises  in  the  follow¬ 
ing  way.  Suppose  that  a  mean  range  refraction  correction  for  a  certain  period  — 
perhaps  several  hours  or  even  several  days  —  is  known  but  that  it  is  desired  to 
improve  upon  this  accuracy.  Such  a  mean  correction  could  be  based  upon  ground 
level  measured  values  of  index  refraction,  but  this  is  not  necessarily  so. 

Now  let  us  suppose  that  the  improvement  in  the  range  correction  which  we 
shall  make  will  be  proportional  to  some  time-filtered  function  of  ,  the  devia¬ 
tions  of  the  surface  value  of  index  of  refraction  about  its  mean  during  this  period. 
Thus,  we  propose  to  make  the  additional  range  refraction  correction 


with  G(t'-t)some,  at  present  unkown,  filter  function.  It  seems  logical  that 
this  function  should  be  chosen  in  such  a  manner  that  the  resulting  range  error 
will  be  a  minimum,  at  least  in  a  statistical  sense.  With  the  introduction  of 
the  space-time  correlation  function  defined  above,  this  may  indeed  be  done. 
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If  we  demand  that  G  be  the  function  which  minimizes  the  spectrum  of  the 
refraction  correction  errors  (assuming  the  mean  correction  to  be  without  error), 
then  we  are  led  to  a  G  which  has  the  Fourier  transform 


Y  (w) 
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where  $  (s,  w)  is  the  Fourier  transform  of  C(s,  r),  (w)  is  the  Fourier  trans¬ 

form  of  c.  (0,  r),  and  the  integration  is  along  the  ray  path  from  the  tracking  station 
to  the  space  vehicle.  The  spectrum  of  the  remaining  range  error  is  then 
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It  is  of  interest  to  note  that  all  the  spectra  which  enter  the  above  equations 
are  cross  spectra  of  refractive  index  fluctuations  between  different  points  along 
the  actual  ray  path  to  the  tracked  vehicle.  Thus,  the  use  of  spectra  obtained 
from  ground  level  measurements  implies  certain  assumptions  about  the  isotropy 
of  the  space-time  correlation  function.  However,  it  is  not  presently  feasible  to 
make  long-period  refractive  index  measurements  along  an  approximate  ray  path. 
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THE  FABRY-PEROT  REFRACTOMETER  AND 
THE  PRECISION  NRL  LYMAN-ALPHA  HUMIDIOMETER 


F.  C.  MacDonald,  D.  L.  Randall,  and  D.  Stillwell* 

INTRODUCTION 

The  measurement  of  microcells  of  atmospheric  moisture,  which  are 
believed  to  be  responsible  for  transhorizon  tropospheric  radio  propagation  at 
VHF  and  UHF  by  the  scattering  mechanism,  has  been  a  challenge  to  the  sci¬ 
entist  because  of  inadequate  meteoroligcal  instriunentation.  An  aircraft  is  the 
most  convenient  platform  for  mounting  these  instriunents,  but,  unless  the  time 
response  of  the  sensing  instriunents  is  very  rapid,  the  speed  of  the  aircraft 
permits  only  large-scale  measurements, 

DISCUSSION 

To  make  fine-grained  meteorological  measurements  from  an  aircraft, 

NRL  is  developing  two  basic  instriunents;  The  Fabry-Perot  Ref racto meter; 
and  the  Precision  NRL  Lyman- Alpha  Humidiometer, 

THE  FABRY-PEROT  REFRACTOMETER*-^^** 

The  Fabry-Perot  refractometer  is  similar  to  the  conventional  Pound 
stabilized  refractometer  except  that  it  uses  a  much  faster  sampling  cavity. 

This  new  resonator  (Figs.  1  and  2)  consists  of  two  equal  radii  spherical  caps 
hollowed  from  1-inch-thick  aluminum  plates.  These  concave  reflecting  surfaces 


*U,  S,  Naval  Research  Laboratory,  Washington,  D.  C, 

**Numbers  in  brackets  denote  References  on  page  193, 
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Fig.  1.  Fabry-Perot  Refractometer 
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Fig.  2.  Open-Sided  Spherical  Cavity 
Resonator 


of  7  inches  radii  of  curvature  are  spaced  5.  5  inches  apart  by  three  alloy  rods 
(Weinschel  LA  685).  Unlike  the  cylindrical  (TEOl)  cavity,  this  cavity  has  no 
perforated  end-plates  to  slow  down  the  ventilation  rate  and  to  introduce  errors 
in  refractive  index  resulting  from  d3Tiamic  temperature  and  pressure  changes 
in  the  measured  sample.  The  selection  of  spherical  surfaces  over  parabolic  or 
plane  surfaces  was  made  because  of  the  relative  insensitivity  of  Q  and  frequency 
to  lateral  displacements  of  the  plate  axes.  This  configuration  has  the  advan¬ 
tage  that  any  condensation  on  the  spherical  surfaces  will  be  rapidly  removed 
by  the  air  flow  across  the  surface.  In  addition,  corrosion  of  the  surface  has  a 
minimal  effect  on  cavity  Q  and  manifests  itself  primarily  in  a  long-term 
frequency  drift  of  the  resonator.  For  this  reason  it  is  not  necessary  to  gold  or 
silver  plate  the  reflecting  surfaces  to  obtain  high  Q  as  in  the  case  of  the  cylindrical 
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cavities.  The  loaded  Q  of  this  aluminum  cavity  is  about  10,  000,  which  is  compar¬ 
able  to  the  conventional  cylindrical  cavity.  Accurate  data  is  not  yet  available 
to  state  the  affect  of  temperature  on  the  performance  of  the  cavity,  though 
preliminary  measurements  indicate  compensation  to  within  ±  2  N-units  over 
the  range  -20°C  to  +40°C.  This  variation  will  not  adversely  affect  the  use  of 
the  instrument  for  the  study  of  fast  refractive-index  changes  measured  from 
an  aircraft. 

THE  PRECISION  NRL  LYMAN-ALPHA  HUMIDIOMETER*-^^ 

The  Precision  NRL  Lyman-Alpha  Humidiometer  is  an  improvement  over 
the  basic  instrument  described  in  the  International  Symposium  on  Humidity  and 
Moisture,  held  in  Washington,  D.  C.  ,  May  20-23,  1963,  and  in  "The  Report  of 
NRL  Progress,  "  June  1963,  in  that  an  alternating-current  recording  channel 
has  been  added  for  recording  fast  fluctuations  between  5  to  100  cps. 

The  block  diagram  in  Fig.  3  shows  the  principal  parts  of  the  basic  instru¬ 
ment.  When  energized  by  the  high-voltage  power  supply,  the  hydrogen  lamp 
emits  Lyman-alpha  radiation  at  1216  A  which  passes  through  the  lithium- 
fluoride  window  in  the  lamp,  across  a  measuring  path,  and  through  a  lithium- 
fluoride  window  of  the  detector  tube  which  is  filled  with  nitric-oxide  gas.  By 
photoionization  processes,  the  radiation  from  the  hydrogen  lamp  causes  a 
current  to  flow  in  the  detector  tube  which  is  connected  to  an  electrometer,  a 
metering  circuit,  and  a  recording  system.  This  current  is  a  logarithmic  func¬ 
tion  of  the  water  vapor  density  in  the  measuring  path.  The  absorption  coefficients 

of  the  different  atmospheric  gases  at  Lyman-alpha  have  been  mentioned  in  previous 

[2] 

papers  and  will  not  be  repeated,  except  to  mention  that  there  is  a  window 
which  permits  selective  absorption  of  water  vapor. 
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Fig.  3.  NRL  Lyman-Alpha  Humidiometer  (Block  Diagram) 


Figure  4  illustrates  the  basic  humidiometer  components  shown  in  Fig,  3.  The 
hydrogen  lamp  is  mounted  opposite  the  nitric-oxide  ion  chamber  and  p>otted 
electrometer  with  a  plastic  spacing  block  between  them.  For  the  humidiometer 
used  in  the  measurements  we  shall  later  discuss,  the  measuring  path  spacer 
block  is  3/4  inch  and  the  air  intake  diameter  is  5/8-inch, 

Figure  5  is  a  view  of  the  probe  in  which  this  unit  is  mounted.  The  probe 
is  then  attached  to  the  side  of  the  aircraft  redome.  From  the  outside  of  the 
radome  surface  to  the  air  intake  is  12  inches  for  this  particular  probe.  The 
connections  on  the  base  of  the  probe  are  for  the  hydrogen  lamp  power,  the 
power  to  the  electrometer,  and  the  output  of  the  electrometer. 
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Components  of  Humidiometer 


air  Witake  to 

MEASURING  PATH 


CABLE  CONNECTORS  FOR: 

(1)  HYDROGEN  LAMP 

(2)  ELECTROMETER  AMPLIFIER 


Fig.  5.  Humidiometer  Probe  Assembly 


Figure  6  shows  two  of  these  probes  mounted  on  the  top  radome  of  a  Navy 
Lockheed  Constellation  aircraft.  The  distance  between  the  air-intake  holes 
of  the  two  humidiometers  is  6  feet.  The  air  intake  of  the  top  probe  is  18  inches 
from  the  outside  of  the  radome.  This  places  the  two  air  intake  holes  in  a 
vertical  line.  Figure  7  is  a  forward  view  of  the  two  humidiometer  probes. 

Figure  8  is  a  schematic  diagram  of  the  humidiometer.  Tube  V5  is  the 
hydrogen  lamp,  and  tube  V4  is  the  nitric-oxide  ion  chamber.  The  output  of  the 

9 

ion  chamber  is  dropped  across  a  4.  7  x  10  ohm  resistor  in  this  circuit.  In 
order  to  achieve  faster  time  response  in  the  electronics,  this  resistor  has  been 

9 

changed  to  1.  0  x  10  ohms.  The  cathode  follower  electrometer  output  goes  to 
a  logarithmic  amplifier  and  then  to  the  recorders.  Tube  V7  is  used  to  control 
the  current  flowing  through  V5.  For  frequency  calibration  of  the  humidiometer, 
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Schematic  Diagram  of  HumicJiometer 


the  secondary  of  an  audio  frequency  transformer  is  connected  in  series  in  the 
cathode  circuit,  and  a  signal  generator  is  used  to  modulate  the  hydrogen  Lyman- 
alpha  beam  received  by  the  nitric-oxide  ion  chamber. 

Figure  9  shows  the  basic  humidiometer  components  plus  the  recording  sys¬ 
tem.  The  DC  recording  system  is  a  presentation  on  a  photoelectric  Visicorder. 
The  AC  recording  system  consists  of  an  audio-signal  generator,  a  variable-gain 
amplifier,  a  chopper,  and  a  tape  recorder.  The  magnetic  tape  is  played  into  a 
Packard-Bell  250  computer. 

DC  Recording  System 

The  DC  recording  system  makes  it  possible  to  get  absolute  values  for 
plotting  soundings  as  is  customarily  done  for  refractive  index,  temperature, 
and  humidity.  Figure  10  presents  a  portion  of  mission  12.  Here  you  see  plots 
of  temperature,  dew-point  temperature,  total-refractive  index,  and  moist- 
refractive  index.  In  this  case,  the  dew-point-temperature  curve  is  used  to 
compute  the  moist-refractive-index  values.  Notice  that  the  differences  in  the 
moist-refractive  index  are  slightly  greater  than  for  the  total  refractive  index 
measured  with  the  refractometer.  This  is  probably  due  to  the  fact  that  in  a 
layer  of  increasing  temperature,  the  temperature  term  of  the  refractive-index 
equation  subtracts  from  the  moisture  term.  The  DC  recording  channel  will 
observe  changes  in  refractive  index  and  moisture  to  the  order  of  10  cps. 

AC  Recording  System 

As  shown  in  the  block  diagram  (Fig.  9),  a  capacitor  blocks  the  DC  from 
the  AC  channel.  This  cut  off  frequency  is  about  1/2  cps.  The  AC  channel  is 
designed  to  cover  a  range  from  1/2  to  100  cps.  This  makes  it  possible  to  sense 
refractive-index  fluctuations  heretofore  unmeasurable  by  our  system.  The  moist 
component  of  refractive  index  is  the  dominant  term  in  the  refractive-index 
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NRL  Lyman-Alpha  Humidiometer  and  its  Recording  Syst* 


PRESSURE  HEIGHT  PRESSURE  HEIGHT 

(THOUSANDS  OF  FEET)  (THOUSANDS  OF  FEET) 


MICROMETEOROLOGY 
MISSION  12-MAY  27,  1964 


N-p  =  TOTAL  REFRACTIVE  INDEX 
N^  =  MOIST  REFRACTIVE  INDEX  TERM 


Fig.  10.  Sample  DC  Record  from  Sounding 
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equation,  having  a  weight  of  4. 14  for  a  1-millibar  change  in  water  vapor,  while 
the  temperature  term  has  a  weight  of  -0.  89  for  a  1-degree  change  in  temperature 
over  the  same  height  interval.  Of  course,  this  presupposes  an  atmosphere  of 
1000  millibars  total  pressure,  a  temperature  of  300°K,  and  a  water  vapor  pres¬ 
sure  of  10  millibars. 

Figure  11  depicts  a  typical  humidiometer  calibration  curve  made  with  a 

9 

humidiometer  having  a  measuring  path  of  3/4  inch  and  a  drop  resistor  of  1  x  10 
ohms.  The  slope  of  this  curve  between  2  to  10  millibars  is  13.  9  microamperes 


mbs  — »■ 

Fig.  n.  Typical  Humidiometer  Calibration  Curve 


per  millibar.  Above  10  millibars  there  is  a  knee  in  the  curve,  and  at  14  milli¬ 
bars,  the  slope  has  lessened  to  5  microamperes  per  millibar.  At  7  millibars, 

_3 

in  the  straight  portion  of  the  curve,  the  noise  level  is  0.9  x  10  millibars,  or 
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3.8x10  ^  moist  N-units  at  300°K.  At  14  millibars  the  noise  level  is  2.  5  x  10 
millibars  or  10,  3x10^  moist  N-units  at  300°K.  These  are  the  basic  sensitivity 
values  for  this  tj^pe  humidiometer.  The  sensitivity  can  be  increased  by  increas¬ 
ing  the  value  of  the  drop  resistor,  but  this  increases  the  time  constant  so  that 
high  frequencies  cannot  be  measured. 

Figure  11  also  presents  the  method  of  calibrating  the  water-vapor  variation 
indicated  by  the  system.  The  Lyman-alpha  source  is  modulated  by  a  10-cps 
sine  wave  of  fixed  amplitude  in  microamperes.  The  tape  recorder  shows  this 
variation  in  bits  (1  volt  equals  62  bits)  which  are  directly  proportional  to  the 
amplitude  of  the  sine  wave.  The  calibration  curve  is  then  used  to  convert  this 
output  to  a  fluctuation  in  millibars  of  water-vapor  pressure. 

With  the  basic  sensitivities  quoted,  it  is  possible  to  amplify  the  incoming 
signals  so  that  they  will  always  fall  within  the  0-  to  2-volt  input  range  of  the 
tape  recorder.  Four  gain  steps  are  provided  so  that  weak  signals  can  be  ampli¬ 
fied  to  fit  the  capabilities  of  the  input  to  the  tape  machine. 

Figure  12,  curves  No.  1,  No.  2,  No.  3,  and  No.  4,  displays  the  modulating 
frequency  at  different  gain  settings.  To  keep  the  output  of  the  amplifier  within 
the  2-volt  range  of  the  tape  recorder  it  is  necessary  to  attenuate  the  input  of 
the  hydrogen  lamp.  In  curve  No.  1  of  Fig.  12,  the  amplifier  is  set  for  gain 
position  1.  Notice  that  there  is  very  little  noise  on  the  envelope  of  the  sine 
wave.  For  this  setting,  0.  3  millibars  equals  1.  2  moist  N-units,  or  62  bits. 

In  curve  No.  2  of  Fig.  12,  there  is  a  slight  increase  in  noise  superimposed  on 
the  10-cps  modulating  frequency  with  the  amplifier  set  in  gain  position  2.  For 
this  position,  0,  15  millibars  equals  0.  62  moist  N-units,  equals  62  bits.  Curves 
No.  3  and  No.  4  of  Fig.  12  describe  the  same  type  of  information.  In  curve 
No.  4  of  Fig.  12,  it  is  very  obvious  what  would  happen  if  the  noise  equaled  the 
amplitude  of  the  incident  signal. 
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Fig.  12.  Gain  Calibration  Curves 


Figure  13  presents  the  frequency  response  of  the  system  in  terms  of 
the  amplitude  at  10  cps,  which  is  100  percent. 


Fig.  13.  Variation  of  Gain  with  Frequency 


Some  Experimental  Results 

Now  that  humidiometer  calibration  procedures  have  been  worked  out  for 
both  the  DC  and  AC  channels  some  measurements  will  be  considered. 

On  Mission  16,  made  July  7,  1964,  the  aircraft  was  slowly  descending 
through  the  top  of  a  haze  layer.  The  zone  of  greatest  change  was  between  8550 
to  7880  feet,  a  layer  670  feet  thick.  The  refractive  index  gradient  as  measured 
by  the  refractometer  was  52  N-units/lOOO  feet,  and  the  temperature  gradient 
was  +  2.  7°  C/1000  feet.  The  oscilloscope  presentations  for  the  two  humidiometers 
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were  fluctuating  in  unison  as  they  normally  do,  but  then  a  layer  was  encountered 
at  8250  feet  in  which  the  oscilloscope  presentations  of  the  probe  fluctuations  were 
out  of  phase. 

The  right-hand  portion  of  Fig,  14  is  a  sketch  of  what  was  seen  on  the 
oscilloscope.  Part  of  the  time  the  two  presentations  were  180  degrees  of  arc 
out  of  phase.  This  condition  lasted  intermittently  for  about  a  minute,  and  then. 
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Fig.  14.  Out-of-Phase  Variation  of  Humidiometer  Probe  Sensors 


as  the  aircraft  lost  altitude,  the  humidiometer  fluctuations  were  in  unison  again. 
When  the  aircraft  ascended,  the  layer  was  found,  and  the  fluctuations  were  still 
out  of  phase.  This  lasted  for  nearly  5  minutes. 
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The  left-hand  portion  of  Fig.  14  is  a  possible  explanation  of  this  phenomenon. 
Suppose  that  there  was  a  filament-like  layer  of  moisture  below  which  the  moisture 
gradient  was  increasing,  and  above  which  the  moisture  was  decreasing.  Suppose 
also  that  this  layer  had  travelling  waves  of  a  sinusoidal  form.  This  would  then 
account  for  the  out-of-phase  variation  of  the  two  humidiometers.  This  would 
mean  that  the  filament  must  be  very  thin  because  the  probes  are  only  6  feet  apart. 

To  get  more  information  about  this  phenomenon,  a  portion  of  this  record 
was  printed  out  by  the  PB-250  computer  and  some  actual  values  measured.  This 
is  shown  in  Fig,  15.  The  assumption  was  made  that  the  mean  water  vapor  density 

of  the  fluctuations  for  either  probe  was  the  same.  Then  by  using  the  gain  and 
frequency  calibrations  described  in  Fig.  12,  moist  values  of  refractive  index 


Fig.  15.  Dimensions  of  a  Microcell  of  Moisture 
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were  assigned  to  the  deflections.  Consider  the  half  cycle  between  5  and  9 
hundredths  of  a  second.  For  humidiometer  1  there  was  an  excursion  of  1.  2  N- 
units,  and  for  humidiometer  2  there  was  an  excusion  of  0.  3  N-units.  Now  if 
the  lower  probe  was  in  a  region  where  the  refractive  index  gradient  was  increas¬ 
ing,  and  the  upper  probe  in  a  region  where  the  refractive  index  gradient  was 
decreasing  we  have  a  minimum  refractive  index  gradient  of  0.  9  N-units/6  feet 
or  150  N-units/1000  feet.  By  knowing  the  horizontal  speed  of  the  aircraft  we 
can  compute  a  minimum  horizontal  gradient  of  90  N-units/1000  feet.  Of  course, 
one  microcell  of  moisture  in  itself  does  not  mean  much,  but  when  there  is  a 
thin  layer  covering  many  square  miles  with  variations  of  refractive  index  as 
noted  in  this  cell,  then  the  problem  becomes  more  meaningful. 

Figure  16  is  a  power  spectra  of  a  6-second  sample  made  233  seconds 
before  the  measurements  shown  in  Fig.  15.  The  magnitudes  of  the  cells  are 


Fig.  16.  Power  Spectra  of  Moisture 
Fluctuations 
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read  from  the  abscissa.  Consider  the  maximum  for  k  equal  to  .  08.  This  states 
that  the  size  of  the  cells  at  this  particular  maximum  were  78.  54  feet.  Consider 
the  maximum  at  k  equal  to  0.  8.  This  gives  a  cell  size  of  7.  8  feet. 


[2] 

The  ordinate  in  Fig.  14  is  (delta  bits)  .  In  Fig.  12,  the  gain  calibration 
curves  show  the  equivalent  of  bits  in  moist  refractive  index  units.  This  6-second 
sample  was  made  on  gain  position  2.  Using  this  information,  the  delta  N  in 

_3 

moist  refractive  index  units  was  about  6.  1  x  10  for  the  78-foot  cell,  and  1.  3  x 

-4 

10  for  the  7.  8-foot  cell.  These  are  horizontal  gradients  made  from  one 
humidiometer.  The  sample  from  which  these  illustrations  are  computed  was  6 
seconds  in  length  and  covered  a  horizontal  distance  of  about  1500  feet. 


CONCLUSION 


The  results  cited  above  are  a  beginning  of  things  to  come.  The  Fabry- 
Perot  Refractometer  when  calibrated  and  tested  will  be  recorded  on  the  visi- 
corder  and  the  tape  recorder  in  the  same  manner  as  the  two  humidiometer 
probes.  The  tape  recording  technique  and  the  Packard-Bell  250  computer  make 
it  possible  to  print  out  fast  fluctuations  of  the  instrument  up  to  100  cps  which 
were  heretofore  not  discernible  from  the  visicorder  record,  and  to  make  power 
spectra  of  samples  of  3  or  6  seconds  in  length. 

So  far  the  humidiometer  is  a  fair  weather  measuring  instrument  because 
it  loses  sensitivity  for  water  vapor  pressure  above  8  millibars.  Plans  are  under 
study  for  making  an  all-weather  humidiometer. 

The  object  of  the  fast  response  instrument  development  program  is  to  make 
studies  of  microcells  in  all  types  of  weather  conditions  so  that  the  meteorologist 
will  have  some  basis  for  forecasting  radio  scattering  conditions  from  a  weather 
map  and  meteorological  soundings  of  temperature,  dew  point,  and  refractive  index. 
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INSTRUMENTATION  FOR  LOWER  TROPOSPHERIC 
SOUNDINGS  OF  REFRACTIVITY 


D.R.  Hay* 

This  paper  describes  one  aspect  of  our  tropospheric  refraction  studies  at 
the  Iftiiversity  of  Western  Ontario.  We  are  concerned  with  special  instrumenta¬ 
tion  for  sounding  the  lower  troposphere.  The  region  of  interest  lies  between 
the  ground  and  3000  feet  above,  the  so-called  layer  of  frictional  influence. 

The  motivation  for  this  work  stems  from  our  research  on  radar  angels. 

These  are  clear-air  reflections  which  are  observed  on  sensitive  microwave  radars. 
We  have  been  studying  microwave  angels  for  several  years  with  a  special 
7000-mc  radar  which  has  two  10-foot  antennas  spaced  15  feet  apart  and  pointing 
in  the  vertical  direction.  Our  attempts  to  interpret  the  anomalous  clear-air 
reflections  lead  us  to  conclude  that  they  arise  from  very  thin  strata  of  contrast- 
ing  temperature  and  humidity  injected  into  the  environment.  *■  ’  Further,  it 
appears  that  they  must  be  of  the  order  of  radar  wavelength  in  depth.  This  means 
that  we  are  looking  for  changes  in  refractivity  of  the  order  of  one  or  two  parts 
per  million,  through  a  vertical  height  interval  of  only  a  few  centimeters.  This 
type  of  structure  may  be  related  to  turbulence  and  other  well-known  weather 
features. 

Direct  evidence  of  air  micro  structure  with  this  degree  of  fineness  is  very 
meager.  Instruments  which  will  measure  temperature  and  water  vapor  pres¬ 
sure  with  sufficient  precision  and  time  or  spacial  resolution  are  generally  not 


♦University  of  Western  Ontario,  London,  Canada. 
♦♦Numbers  in  brackets  designate  References  on  page  205 . 
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commercially  available.  We  have  been  working  on  a  program  on  instrument 

[3  4  5  6] 

development  to  aid  us  in  obtaining  these  measurements  for  some  time,  J 

The  measurements  are  made  from  a  ground  station  which  is  adjacent  to 
our  microwave  radar,  and  which  permits  the  same  volume  of  air  to  be  examined 
simultaneously  by  different  techniques.  Two  types  of  instruments  are  used.  One 
is  a  smoke  generator  in  the  nose  of  a  rocket  which  deposits  an  instantaneous 
vertical  column  that  is  studied  for  wind  shear.  The  second  instrument  is  a 
combination  which  measures  air  refractivity  and  air  temperature  simultaneously, 
with  a  rapid  speed  of  response  and  a  high  degree  of  spatial  resolution.  The 
information  from  these  instruments  will  be  used  not  only  to  deduce  the  detailed 
profiles  of  temperature  and  water  vapor  pressure,  but  also  to  obtain  the  local 
Richardson  numbers  to  describe  the  regime  of  turbulence  in  this  lower  part  of 
the  troposphere. 

The  remainder  of  my  talk  will  be  on  the  details  of  these  instruments. 

Figure  1  shows  the  rocket-bome  smoke  generator  which  deposits  the  vertical 
smoke  trail.  The  rocket  vehicle  commercially  known  as  a  "Cricket,  "  is 
propelled  by  carbon  dioxide  in  both  liquid  and  gaseous  form.  The  launching 
device  is  readily  p>ortable.  It  will  be  seen  that  the  man  on  the  left  is  holding  a 
launching  tube  extension,  which  is  about  10  feet  long.  After  launch,  the  rocket 
vehicle  ascends  for  about  12  seconds  and  then  turns  over  to  begin  its  descent  under 
parachute.  Figure  2  shows  the  nose  cone  being  fitted  into  the  rocket  vehicle. 

This  smoke  generator  is  not  commercially  available;  it  is  similar  to  one 

[7] 

described  and  used  by  Gill  at  the  University  of  Michigan,  Figure  3  is  a 
close-up  of  the  smoke-depositing  nose  cone.  It  is  10  inches  in  length  and  has 
a  pressure  chamber,  in  the  upper-half,  for  gaseous  carbon  dioxide  at  70  psi. 

The  clear  lower  section  contains  titanium  tetrachloride  which  is  ejected  through 
fine  nozzles  to  produce  titanium  oxide  when  in  contact  with  moisture  in  the  air. 
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Fig.  1.  Cricket  Rocket  Being  Lowered  into  the  Launching  Tube 


Fig,  2,  Fitting  the  Nose  Cone  on  the  Rocket  Vehicle 
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Fig.  3.  Smoke-Generating  Nose  Cone 


It  is  necessary  to  perform  two  tasks  after  the  ascent  of  the  rocket,  as 
shown  by  Figs.  4  and  5.  Since  the  rocket  is  recoverable  and  reusable,  the 
impact  point  must  be  determined  by  tracking  the  rocket  during  descent.  This 
is  done  with  the  aid  of  an  optical  range  finder  at  the  launching  site  and  an  astro- 
compass  at  a  remote  location  (sere  Fig.  4).  The  second  task  involves  photographic 
recording  of  the  smoke  trail  during  its  10-second  lifetime.  Two  cameras  are 
used  for  this  purpose.  One  is  a  photo-theodolite  with  a  telephoto  lens  to  observe 
fine  details  in  the  trail  deformation.  The  other  has  a  wide-angle  lens;  here  the 
film  advances  at  1-second  intervals  to  record  the  time  sequence  in  the  position 
of  the  whole  smoke  trail.  These  cameras  are  located  normal  to  the  plane  of 
drift  of  the  trail  as  determined  from  a  preliminary  pilot  balloon  ascent.  Figure 
5  shows  the  two  cameras  at  a  field  site. 
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Fig.  4.  Optical  Ranging  Instruments 


Fig.  5.  Photographic  Equipment  for  Observing  the  Smoke  Trail 
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The  second  instrument  for  airborne  soundings  measures  air  temperature 
and  refractivity  simultaneously.  It  is  launched  on  a  manner  similar  to  a  stand¬ 
ard  radiosonde,  with  balloon  and  parachute  (see  Fig.  6).  It  will  be  seen  that  the 
trailer  contains  a  winch  for  controlling  the  height  of  ascent;  the  ascent  rate 
normally  is  about  1000  feet  per  minute.  The  response  of  the  instrument  is  ade¬ 
quate  to  follow  changes  in  temperature  or  refractivity  within  a  height  of  interval  of 
only  a  few  centimeters.  Figure  7  shows  the  mobile  ground  station  with  a  4-foot 
paraboloid  antenna  on  its  roof  for  receiving  and  recording  the  signal  from  the 
airborne  instrument. 


Fig.  6.  Launching  the  N-T  Sonde 
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Fig.  7.  Recording-Receiver  Van 


Two  views  of  the  N-T  sonde  are  found  in  Figs,  8  and  9.  This  device  weighs 
about  4  pounds  and  is  built  in  modular  form.  It  contains  two  400-mc  transmitters; 
one  carries  the  temperature  information  and  the  other  the  refractivity  modulation. 
The  section  above  the  transmitters  contains  the  temperature  sensor,  which  uses 
a  very  fine  element  of  platinum  wire.  It  can  detect  changes  in  temperature  of  a 
fraction  of  a  degree  with  a  response  time  of  one  millisecond.  At  the  top  is  the 
capacitor-type  (Hay)  refractometer.  The  sensing  element  is  a  ring-type  air 
capacitor  with  a  vertical  depth  of  1  cm,  for  the  desired  spatial  resolution.  All 
of  the  circuits  of  this  dual  instrument  have  been  transistorized,  with  the  excep¬ 
tion  of  the  400-Mc  transmitters. 

Figure  10  is  a  block  diagram  of  the  N-T  sonde.  The  upper  section  repre¬ 
sents  the  refractometer  circuit;  the  lower  section,  the  temperature  circuit,  and 
the  middle,  a  common  control  circuit.  The  air  capacitor  on  the  left  which  is  the 
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Fig.  9.  N-T  Sonde,  Showing 

Refroctometer  Sensing  Capacitor 
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Fig.  10.  Block  Diagram  of  the  N-T  Sonde 

refractivity-sensing  element,  forms  part  of  the  resonant  circuit  of  a  9-mc 
Clapp  oscillator.  Changes  in  this  capacitance,  which  are  due  to  changes  in  the 
dielectric  constant  of  the  air  passing  through  it,  result  in  variation  of  the  oscil¬ 
lator,  and  the  audio  difference  output  used  to  modulate  one  of  the  400-mc  trans¬ 
mitters.  The  temperature  sensing  element  is  platinum  resistance  wire,  having 
a  diameter  of  2.  5  microns.  It  is  part  of  the  circuit  of  a  unijunction  transistor 
multivibrator.  The  output  of  this  circuit  also  is  an  audio  signal  which  is  used  to 
modulate  the  second  400-mc  transmitter.  In  each  case,  the  sensing  element  is 
replaced  periodically  by  a  reference  element,  for  a  brief  interval,  for  calibration 
purposes.  * 


*The  calibration  includes  temperature  compensation  of  the  sensing  element  for 
the  ambient  air  temperature.  It  should  be  noted  that  the  simplicity  of  construc¬ 
tion  has  been  an  important  factor  in  the  design  of  the  N-T  sonde,  and  that  circuit 
complexity  is  similar  to  that  of  the  conventional  radiosonde.  This  suggests  its 
possible  use  as  an  expendable  instrument. 
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At  present,  we  are  studying  the  photographic  techniques  for  smoke-trail 
observation,  and  the  calibration  of  the  refractometer  sensor.  Examples  of  the 
temperature  profile  in  the  lower  troposphere  that  were  obtained  with  the  tempera¬ 
ture  sonde  are  shown  in  Fig.  11;  the  small  fluctuations  superimposed  upon  the 
slowly  varying  bacl^ound  are  apparent. 
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Profiles  of  Air  Temperature 
Obtained  with  the  T  Sonde 


A  number  of  preliminary  soundings  with  the  instruments  have  been  carried 
out  in  the  field.  The  separate  components  have  been  tested  individually,  and 

[4  6] 

reports  on  part  of  this  work  have  been  published  elsewhere. 

The  research  and  development  of  the  two  airborne  instruments  described 
in  this  paper  have  been  supported  in  part  through  a  subcontract  from  the 
Northern  Electric  Company  of  Canada,  under  Contract  DDP  9  PC  1-436, 
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